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PROCEEDINGS OF THE SPACE SHUTTLE 


ENVIRONMENTAL ASSESSMENT WORKSHOP 
ON STRATOSPHERIC EFFECTS 

9 , 

By Andrew E. Potter, Compiler 


1. SUMMARY 


Representatives of the various National Aeronautics and Space Administra- 
tion centers and supporting contractors presented working papers assessing the 
impact of Space Shuttle exhaust on the environment. Participants in the work- 
shop reviewed and discussed the working papers and made the following conclu- 
sions . 


Two physical effects on the stratospheric environment might result from 
Space Shuttle exhaust products. 

1. Alxmiinimi oxide dust particles might alter the radiation balance of the 
stratosphere and modify weather. 

2. The gaseous hydrogen chloride and nitric oxide and the aluminum oxide 
dust might increase the rate of ozone destruction and lead to a reduction of 
ozone. The average levels of ultraviolet radiation at the surface of the Earth 
would then increase, with possible adverse biological consequences. 

A preliminary evaluation indicated that the only significant effect would 
probably be ozone depletion by hydrogen chloride. Consequently, the major 
emphasis was on a more detailed investigation of hydrogen chloride and its 
effect on the ozone layer. A parallel investigation of low-chlorine propellant 
combinations based on ammonium nitrate as possible substitutes for the base- 
line propellant was pursued, should it be necessary to reduce the chlorine 
content of Space Shuttle exhaust. 

A review of ozone-reduction predictions from the various stratospheric 
models indicated that a reasonable estimate of the ozone reduction from Space 
Shuttle operation was 0.2 percent (hemi spherically averaged). The biological 
consequences of this reduction were estimated to be too small to be realistic- 
ally predicted with the current state of knowledge. 


^ASA Lyndon B. Johnson Space Center. 
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2. INTRODUCTION 


The Space Shuttle Environmental Assessment Workshop was held at the NASA 
lyndon B. Johnson Space Center (JSC) to assess the impact of Space Shuttle ex- 
haust on the environment. Working papers were presented on the following top- 
ics: propellant characteristics, ozone (O^) depletion, particulates, plume 

afterburning, and medical/biological effects. These papers are included in 
their entirety in appendixes A to H. The attendees (listed in table 2-l) re- 
viewed and discussed these topics; the results of their reviews are summarized 
in this section. 

As an aid to the reader, where necessary the original units of measiire have 
been converted to the equivalent value in the Systeme International d'Unites 
(si). The SI mits are written first, and the original units are written par- 
enthetically thereafter. 


PROPELLANT CHARACTERISTICS 


The propellant formulation currently planned for use by the Space Shuttle 
solid rocket booster is an aliminiim/ammonium perchlorate (Al/AP) combination 
that releases aluminum oxide (Al^O^) particles and gaseous hydrogen chloride 

(HCI) into the stratosphere dirring la\inch burn. ’’Frozen" chemical equilibrium 
in the rocket nozzle and afterbwning with ambient air entrained in the exhaust 
plume lead to the production of nitric oxide (NO). 


OZONE DEPLETION BY HCI 


Stratospheric models were used to predict the effects of HCI on 0^ deple- 
tion. The simplest models are one dimensional; that is, they consider only 
vertical mass transport. Transport by horizontal winds is important, since 
Space Shuttle exhaust products will be introduced into the stratosphere at only 
two laiinch points, NASA John F. Kennedy Space Center, Florida; and Vandenberg \ 

Air Force Base, California. Latitudinal winds will then spread the products 
in a band around the globe after only a few days. However, north/south treins- 
port is considerably slower. Most stratospheric 0^ depletion predictions have 

employed one-dimensional models in which the slowness of transequatorial trans- 
port was approximated by using chlorine (Cl) concentrations averaged over the 
Northern Hemisphere. The Space Shuttle 0^ depletion effects predicted by five 

different one-dimensional models were compared in detail. Predictions agreed 
within a factor of 3. 

A consensus was reached that the substantial differences in predictions 
from various models could be explained by systematic differences in assumed 
input parameters. Two types of iincertainties should be considered for the 
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input parameters. One is systematic; it involves factors such as incorrect 
initial assumptions on transport coefficients, errors in the laboratory-measured 
coefficients, or the failure to include important reactions. This type of un- 
certainty cannot be evaluated by statistical methods but rather must be deter- 
mined by stratospheric measurements that compare model predictions with obser- 
vation. These measiarements have not yet been made for the Cl chemistry asso- 
ciated with the Space Shuttle exhaust. Therefore, the systematic errors in the 
model predictions cannot be fully evaluated at this time. Some indication of 
the possible range of systematic errors might be gained from comparison of re- 
sults from different stratospheric models. These errors range over a factor 
of 3. Random errors, such as standard deviation associated with rate-constant 
determinations, are susceptible to statistical analysis; for one model, these 
random errors lead to a factor of 3 at a 50-percent confidence level. 

The range of reasonable Space Shuttle 0^ depletion predictions is from 0.1 

to 0.3 percent for 60 launches per year averaged over the Northern Hemisphere. 
This range represents a lower limit to the uncertainty in the perturbation. An 
additional probability distribution of a factor of 3 on either side can be at- 
tached to any chosen centerline value. For planning purposes, the JSC Environ- 
mental Effects Project Office chose a centerline value of 0.2 percent. This 
value fits within the random-error bounds of all model predictions considered 
and represents a reasonable compromise on the systematic errors associated with 
the choice of assumed input parameters. 

Preliminary calculations using multidimensional models have provided an 
approximate picture of the north/south distributions of total chloride and its 
depletion effects. Approximately 30 percent of the Cl eventually reaches the 
Southern Hemisphere in the steady state (rather than 0 percent as assumed for 
the one-dimensional models). Preliminary calculations also show a seasonally 
changing "corridor" effect, such that the 0^ depletion near latitude 30® N is 

2 times the Northern Hemisphere average during the summer and 1.5 times the 
average during the spring and fall. 

The 0^ reduction effects estimated for the various Space Shuttle exhaust 

products (one-dimensional models, hemispheric ally averaged products) are as 
follows . 


Propellant 

Depletion by 
HCl, percent 

Depletion by 
oxides of 
nitrogen, 
percent 

Depletion by ^120^ 
particles, percent 

Baseline (21 percent HCl) 

0.2 

<0.001 

<0.001 

High-Cl alternate 
(5.8 percent HCl) 

• 

o 

00 

<0.001 

<0.001 

Low-Cl alternate 
(3 percent HCl) 

.Olt 

<0.001 

<0.001 
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EFFECTS OF PARTICULATE EXHAUST EMISSIONS 


The Al^O^ particles in the exhaust plume could function as catalytic agents 

for 0^ decomposition or they coiild alter the radiation balance of the Earth by 

reflecting siinlight. Both of these effects have been analyzed and found to be 
negligible. The AlgO^ particles could also function as nucleation centers for 

svilfuric acid (HgSO^) aerosols which occur naturally in the stratosphere. Any 

change in the natural H 2 S 0 j^ aerosols could be climatologically significant. 

(This problem is currently being studied. ) 

PLUME AFTERBURNING EFFECTS 


"Frozen" chemical equilibriimi in the nozzle and an afterbiirning of the ex- 
haust plume produces NO, which could cause additional 0^ depletion. Conclusions 
on this subject were as follows. 

1. The oxides of nitrogen (NO^) production rate for the baseline propel- 
lant in the stratosphere is two orders of magnitude lower than the HCl produc- 
tion rate. 

2. The alternate propellant "A" NO^ production rate in the stratosphere 

is as much as 10 times greater relative to HCl than for the baseline propellant. 

In any case, whether the baseline or alternate propellants are used, the 0^ de- 
pletion effect of afterburning will be at least an order of magnitude less than 
that caused by the HCl content. It is consequently considered to be negligible 
and will not be further considered. 


MEDICAL/BIOLOGICAL EFFECTS OF ESTIMATED 0^ REDUCTION 

The estimated 0^ reduction can be converted by a series of assumptions and 

calculations into a rough estimate of the increase in biologically harmful ul- 
traviolet (BHuv) radiation from the Sm. It is generally assumed that X per- 
centage of 0^ reduction resiilts in an approximately 2X increase in BHuv ra- 
diation for small 0^ reductions. Thus, a 0.1- to 0. 3-percent reduction in 0^ 

results in a 0.2- to 0.6-percent increase in BHuv radiation. A O.U-percent in- 
crease in BHuv radiation was chosen for assessing the potential impact on the 
biosphere. 
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Based on the limited available biological data, the impact on the biosphere 
of an assumed O.U-percent increase in BHuv radiation has been assessed, and the 
following conclusions were drawn. 

1. The effects on the biosphere of a 0.1*-percent increase in BHuv radia- 
tion will not be detectable even after decades of observation; the natural 
BHuv irradiances are highly variable, and so are the responses of organisms 
and ecosystems to given doses of BHuv radiation at given dose rates. However, 
the fact that a cause-and-effect relationship of an increase in BHuv radiation 
cannot be established does not mean that an increase in BHuv radiation can be 
ruled out as a contributor to some deleterious future event ; lack of detecta- 
bility is not equatable with a lack of effect. 

2. There is a very low probability that an average O.i^-percent increase 
in BHuv radiation will have any unacceptable effect on agricultural plants or 
natural ecosystems. This conclusion is based primarily on the wide range of 
BHuv irradiances to which agricultiiral and natural ecosystems are currently 
exposed without apparent impact and on the fact that more important and per- 
vasive factors (e.g., temperature, precipitation, interspecies competition) 
generally limit the geographical distribution of vegetation types. 

3. Some increase in the number of skin cancer cases may occur among sus- 
ceptible individuals as a result of a 0.2-percent 0^ reduction; however, be- 
cause of the great variability of normal BHuv radiation and biological re- 
sponses mentioned previously, the long latent period (20 to 60 years) for the 
appearance of skin tiomors , and the many other factors that already may be 
tending to increase or decrease the number of skin cancer cases (both reported 
and vmreported), the increase will not be detectable, 

k. The number of skin cancer cases resulting from a 10-year-duration, 

0.2-percent 0^ reduction is not realistically predictable, (The widely publi- 
cized "prediction" that a 1-percent decrease in stratospheric 0^ will lead to 

approximately 6000 more cases of skin cancer each year in the United States 
is based on questionable data and refers to a lifetime exposure). The 
extant experimented, and epidemiological data, although clearly suggesting a 
contribution of solar radiation exposure to skin cancer incidence, are in- 
adequate for making quantitative predictions with reasonable limits of un- 
certainty for the Shuttle case. 

Three reasons why the data are inadequate for quantitatively predicting 
the additional nvimber of individuals developing skin cancer as a result of a 
0.2-percent 0^ reduction are as follows. 

l. The true current baseline annual incidence is not accurately known, 
and the number of cases reported does not necessarily reflect the number of 
cases that occur; moreover, the reported incidence is increasing rapidly even 
in the absence of any apparent 0^ reduction. 

2. The reported annual incidence data apply to cases, not to patients; a 
sizable niomber of new cases occur in previously affected susceptible individuals. 
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3. The analyses of the epidemiological data do not adequately consider 
other contributing factors such as carcinogenic chemicals in the environment, 
changes in lifestyle (increased outdoor living and changes in clothing and 
hair style), geographic location of domicile prior to the development of skin 
cancer, and ethnic character. 

Finally, the basic epidemiological data for skin cancer incidence are 
obviously for a full lifetime exposure to sunlight. The induction period for 
skin cancer is 20 to 60 years, averaging near hO years. The Shuttle effects 
will last about 10 years. Current predictive models for skin cancer do not 
take the exposure period into account eind therefore cannot be applied correctly 
to the Shuttle case. 
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TABLE 2-1.- ATTENDEES OF THE WORKSHOP 


Attendee 

Agency 

J. Q. Miller 

NASA MSFC^ 

J. B. Stephens 

NASA MSFC 

R. D. Hudson 

NASA GSFC^ 

R. S. Stolarski 

NASA JSC 

C. E. Baker 

NASA Lewis Research Center 

R. D. Rundel 

NASA JSC 

S. C. Liu 

University of Michigan 

D. M. Butler 

Rice University 

D. E. Robbins 

NASA JSC 

G. Varsi 

NASA JPL*^ 

A. E. Potter 

NASA JSC 

J. Stanley 

NASA JSC 

G. L. Pellett 

NASA Langley Research Center 

R. Stanley 

NASA Langley Research Center 

I. G. Poppoff 

NASA ARC'^ 

M. H. Farlow 

NASA ARC 

R. C. Whitten 

NASA ARC 

0. B. Toon 

NASA ARC 

D. Davis 

University of Maryland 


George C. Marshall Space Flight Center. 
^Goddard Space Flight Center. 

^Jet Propulsion Laboratory. 

'^Ames Research Center. 
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TABLE 2-1.- Concluded 


Attendee 

Agency 

R. Cicerone 

University of Michigan 

J. S. Chang j 

LLL ERDA® 

D. S. Nachtwey 

NASA JSC 

L. R. Greenwood 

NASA Langley Research Center 

P. Wetzel 

NASA Headquarters 

R. T. Watson 

jpl'^ 

L. Keyser 

JPL 

W. B. Moore 

JPL 

1 

P. J. Armitage 

NASA JSC 

H. S. Pergament 

AeroChem 

B. G. Cour-PaJLais 

NASA JSC 

D, J. Kessler 

NASA JSC 

R. P. Turco 

R&D Associates 

C. Riegel 

San Jose State University 

Ro Prinn 

MIT^ 

Pc B. Burband 

NASA JSC 

H. Rudolph 

NASA KSC® 

J. King 

NASA Headquarters 


Q 

Jet Propulsion Laboratory. 

0 

Lawrence Livennore Laboratory Energy Research and 
Development Administration. 

f 

Massachusetts Institute of Technology. 
s 

°John F. Kennedy Space Center. 
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3. BASELINE AND ALTERNATE PROPELLANTS: 


BALLISTICS AND EXHAUST COMPOSITION 


A study of low-Cl solid propellants was initiated in August 197^ with the 
objective of defining possible alternates to the baseline propellant. Complete 
replacement of the perchlorate propellant with a nitrate-based propellant causes 
unacceptable ballistic performance. Consequently, a composite propellant grain 
was considered in which the baseline propellant is burned out at a maximum al- 
titude of 20 kilometers (63 000 feet). Then, the burn changes to an alternate 
low-Cl formulation for passage through the stratosphere. 

Approximately 250 different possible alternate propellant formulations were 
evaluated. These formulations included various combinations of AP and non-HCl- 
producing oxidizers such as ammonium nitrate (AN), cyclotetramethylenetrinitra- 
mine (HMX), trimethylethanitrinitrate , nitrocellulose, nitroglycerine, and 
cyclotrimethylenetrinitramine. Formulations were first evaluated on a theoretical 
basis in comparison with the baseline AP-oxidized propellant. From this com- 
parison, the most promising alternate propellant formulations were selected for 
preliminary process, compatibility, hazard, burn rate, pressure exponent, and 
density impulse characterization. Based on results from these characteriza- 
tions , the best HMX- and non-HMX-containing alternate propellant formulations 
were selected for additional evaluation. A formulation and ballistic property 
comparison of these two alternate formulations ("A" and "B") with the baseline 
formulation is presented in table 3-1. To achieve acceptable ballistic prop- 
erties, both alternate formulations require the addition of AP. The theore- 
tical exhaust composition for the baseline and alternates "A" and "B" is pre- 
sented in table 3-II. 

The most promising design concept developed was a double-web star propel- 
lant grain configuration (dual-grain concept). With this concept, the baseline 
propellant is completely burned out at an altitude of approximately 20 kilo- 
meters (63 000 feet) and the alternate propellant (low HCl) is burned at alti- 
tudes above 20 kilometers (63 000 feet); thus, minimum cost and performance 
impacts result. Alternate propellant formulations "A" and "B" were separately 
used with the baseline propellant to design two booster solid rocket motors in- 
corporating the dual-grain concept. Figures 3-1, 3-2, and 3-3 contain data on 
thrust as a fimction of time and mass flow characteristics for the baseline, 
alternate "A," and alternate "B" designs, respectively. 

Results from preliminary burn rate tests conducted on the alternate "B" 
propellant indicate that AP levels as high as 20 percent by weight may be re- 
quired to obtain acceptable burn rates and combustion efficiencies. This change 
will increase the HCl percentage weight in the exhaust from 3.1 to approximately 
5.8. Other exhaust constituent percentage weights are not significantly af- 
fected by this change. Therefore, for ciirrent environmental assessments, it 
is recommended that an HCl percentage weight range of 3.1 to 5.8 be used for 
alternate propellants. Use of this concept results in a ballistic performance 
loss estimate of approximately 3 percent for the entire grain (i.e., 25 percent 
of the grain consists of propellants having approximately 10 percent less spe- 
cific impulse). 
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TABLE 3 - 1 ..- PROPELLANT FORMULATION AND EALLISTIC PROPERTY DATA 
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TABLE 3- I I.- THEORETICAL PROPELLANT EXHAUST 


COMPOSITION 
[Percentage weight] 


Constituent 

Baseline 

Alternate 

"A" 

"B" 

HCl 

20.9 

3.1 

3.1 

Molecular hydrogen (H^) 

2.1 

3.7 

2.9 

Water (H^O) 

9.1* 

6.1 

16.2 

Carbon monoxide (CO) 

2i*.l 

31.8 

19.6 

Carbon dioxide (CO^) 

3.5 

3.8 

6.5 

Molecular nitrogen (N^) 

8.7 

22.9 

23.2 

AI 2 O 3 

30.2 

28.2 

28.3 

Other 

1.1 

.1* 

.2 
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Vacutim thrust, N (lb) 


17.8 X lO'^ (4.0 X lO"^) 



Time, 

sec 

Altitude, 
m (ft) 

Mass flow rate, 
kg/sec (lb /sec) 

0 

0 

0 

10 

305 (1 000) 

5 263 (11 603) 

20 

1 372 (4 500) 

5 441 (11 996) 

30 

2 896 (9 500) 

4 735 (10 439) 

40 

4 877 (16 000) 

4 166 (9 185) 

50 

7 315 (24 000) 

4 050 (8 929) 

60 

10 363 (34 000) 

4 091 (9 020) 

70 

14 326 (47 000) 

4 198 (9 255) 

80 

18 288 (60 000) 

4 396 (9 692) 

^83. 36 

19 507 (64 000) 

4 301 (9 483) 

90 

22 860 (75 000) 

4 232 (9 329) 

100 

28 651 (94 000) 

3 833 (8 450) 

110 

34 442 (113 000) 

3 456 (7 620) 

120 

41 453 (136 000) 

505 (1 114) 

12 5 

44 196 (145 000) 

115 (253) 


^Transition to alternate propellant. 


Figure 3-1.- Baseline thrust as a function of time and mass flow 

characteristics . 
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Vacuum thrust, N (lb) 


Time, 

sec 

^ Altitude, 

m (ft) 

0 


0 

10 

305 

(1 000) 

20 

1 372 

(4 500) 

30 

3 048 

(10 000) 

40 

5 182 

(17 000) 

50 

7 620 

(25 000) 

60 

10 363 

(34 000) 

70 

14 326 

(47 000) 

80 

18 288 

(60 000) 

90 

23 774 

(78 000) 

100 

28 956 

(95 000) 

110 

35 357 (116 000) 

120 

42 672 (140 000) 

12 5 

45 720 (150 000) 


5 325 
5 478 
4 751 
4 274 
3 947 




0 

1 

740 

2 

076 

0 

475 

9 

423 

8 

702 

8 

972: 

9 

469 











Vacuum thrust, N (lb) 


17.8 X 10'^ (4.0 X lO"^) 



Time, 

sec 

Altitude, 
m (ft) 

Mass flow rate, 
kg/sec (Ib/sec) 

0 


— 

0 



10 

305 

(1 000) 

5 263 (11 603) 

20 

1 372 

(4 500) 

5 441 (11 996) 

30 

2 743 

(9 000) 

4 735 (10 439) 

40 

4 724 

(15 500) 

4 166 

(9 185) 

50 

7 163 

(23 500) 

4 050 

(8 929) 

60 

10 211 

(33 500) 

4 091 

(9 020) 

70 

14 326 

(47 000) 

4 198 

(9 2 55) 

80 

18 288 

(60 000) 

4 396 

(9 692) 

83.4 

19 507 

(64 000) 

3 472 

(7 654) 

90 

22 860 

(75 000) 

3 481 

(7 675) 

100 

28 346 

(93 000) 

3 248 

(7 160) 

110 

34 138 (112 000) 

3 039 

(6 700) 

120 

40 538 (133 000) 

2 694 

(5 940) 

130 

46 634 (153 000) 

328 

(723) 

135 

49 682 

(163 000) 

111 

(262) 


a 

Transition to alternate propellant. 


Figm-e 3-3.- Alternate ”B" thrust as a function of time and mass flow 

characteristics. 
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h, STRATOSPHERIC MODELING 


Compo\inds of Cl and nitrogen (N) catalyze the reaction of 0^ with atomic 
oxygen (O). Therefore, HCl and NO^ released from the Space Shuttle exhaust in- 
to the stratosphere can result in a net reduction of 0^ and consequently in- 
crease the amount of solar uv radiation at the Earth surface. 


The 0^ 


depletion effects of exhaust products may be estimated by use of 


theoretical stratospheric models. In these models, two major factors are con- 
sidered: (l) the transport of chemical compounds into and out of the strato- 

sphere and (2) chemical reactions in the stratosphere. Transport coefficient 
input data for the models are estimated by the analysis of vertical concentra- 
tion profiles through the stratosphere. The relative importance of various 
chemical reactions in the models is estimated from laboratory measurements and, 
in special cases, from stratospheric measiirements of the concentrations of in- 
termediate species, such as hydroxyl (OH). 


In the sections which follow, the current situation with respect to chem- 
ical rates and stratospheric measurements will be surveyed. Following this 
survey will be discussions of various models and their 0^ depletion predictions. 


STRATOSPHERIC REACTION RATES 


The complex network of chemical reactions involved in 0_ depletion is il- 

^ 1 

lustrated in figirre ^-1, wherein C is carbon, F is fluorine, D is an energy 
state, Hv is photon energy, and M is an energy-providing molecule. Not all 
the possible reaction processes are shown i — only those believed to be most im- 
portant. . Recently meas\ired rate expressions for l8 reactions are listed in 
table U-I (refs. U-1 to U-5). This table also shows the estimated precision 
(as opposed to accuracy) where such information is available. Additional data 

are contained in reference 4-6 and an unpublished article.^ 

Despite recent advances in reaction rate measiarements, not all the required 
data are available. Improved data are particularly needed for reactions in- 
volving HO^, of which the most important are two: 


OH + HO^ ^ H^O + 0^ 

(4-1) 

0 + HO^ -»■ OH + 0^ 

(4-2) 


T. Watson: Rate Constants of CIO^ of Atmospheric Interest. To be 

published in the Jo\irnal of Physical Chemistry. 
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Programs to measure these rate constants are underway in several 
laboratories . 

Other reactions for which improved measurements are needed are those in- 
volving formation and destruction of chlorine nitrate (CINO^), oxidation of 

methane (CH^), and formation of HCl. Specific reactions include the following. 


NO + HO,, N0_ + OH 

(i*-3) 

Cl + CHj^ HCl + CH^ 

(U-U) 

Cl + HO,, HCl + Op 

(l*-5) 

CH^Op + NO ^ products 

(h-6) 

CH^Op + HOp products 

(M) 

CIO + NOp + M ^ CINO^ + M 

(l*-8) 

CINO^ + 0 -► products 

(i*-9) 

CINO^ + OH products 

(U-IO) 


Photo chemi cal studies of the following processes are also needed. 


CINO^ + hv products not defined (U-ll) 
HOCl + hv ^ products not defined (h-12) 
HO^ + hv -»■ products not defined (U-13) 


where h 


CH_0 + hv 

is the Planck constant and 


CO + 

H + HCO 

V is photon frequency. 


(U-ilt) 


h-2 


fill 


fill 
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A continuing investigation should be maintained to identify any possible 
reactions of significance which may have been omitted from the present strat- 
ospheric chemical models. 


MEASUREMENT OF STRATOSPHERIC MINOR CONSTITUENTS 


Stratospheric measurements are needed to provide experimental estimates 
of eddy diffusion coefficients, to narrow the uncertainty range of importeint 
reaction rates, and to evaluate the general validity of the stratospheric 
models. In general, transport terms can be derived from measurements of the 
vertical concentration profiles of sources and sinks which have clearly de- 
fined loss rates in the stratosphere (as CH^ and 1 ^q)» although terms derived 

in this way do depend on the details of the model used. 

Important reaction rate information has been derived from measurements of 
unstable species, such as OH. Measiirements aimed at verification of the com- 
plex of reactions shown in figure U-1 have not been made, with the possible 
exception of NO/nitrogen peroxide (NO^). For this reaction, simultaneous 

measurements of associated species must be made. Thus, for the Cl chemistry 
associated with the Space Shuttle problem, simultaneous measurements of CIO, 
HCl, OH, NO, 0^, 0, uv flxix, and a transport tracer, such as CH|^ or nitrous 

oxide (NgO), are needed. 


Chlorine Compounds 

Chlorine conrpounds to be measured are HCl, CIO, Cl, halocarbons (both nat- 
ural and manmade), and particiilate chloride. Measurements of HCl have been 
done; these show increases from less than 0.1 part per billion by volume (ppbv) 
near an altitude of 15 kilometers (i+T 000 feet) to 0.5 to 1 ppbv at an altitude 
of 25 to 30 kilometers (82 000 to 98 000 feet). This altitude profile is at- 

tributed to a stratospheric soiarce (probably photodissociation of halocarbons) 
and a tropospheric sink (rainout). No satisfactory measxirements of CIO or Cl 
exist at this time, although CIO measurements are in progress by several 
investigators . 

Altitude profiles have been measured for several halocarbons (e.g., CF^Cl^ 

and CFCl^) with results that are consistent with photodissociation at altitudes 

above 20 kilometers (63 000 feet). Particulate matter suspended in the strat- 
osphere has been analyzed and found to contain negligible amounts of chloride. 
Consequently, these particles do not appear to be important sources or sinks 
for Cl. 
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Hydrogen Compounds 


Vertical prol'lJes for H^O and CHj^ have been measured by several investiga- 

’.orp; the CHj^ profiles have been used in the estimation of transport coeffi- 

cieri'.'. Different investigators have measured OH in different altitude ranges. 
The most significant measurements relative to the stratospheric problem (ref. 
U-T ) are in the range of 30 to ^3 kilometers (98 000 to lUO 000 feet). These 
results have been used to "tune" some stratospheric chemistry models by ad- 
justing uncertain rate constants to match the observed altitude profile. Re- 
cent ground-based measurements (R. Eiu-nett, private communication) are signi- 
ficant because they show the total OH concentration (vertical column content) 
to be "patchy," varying by almost an order of magnitude. 


Nitrogen Compounds 

Odd-N composites (NO, NO^, and nitric acid (HN0_)) are produced by reaction 

of 0”"D with N^O. There are n\jmerous measurements at different times and places 

(O' NO, NO 2 , HNO^, and Ng, but in situ simultaneous measurements are lacking. 

Below an altitude of approximately Uo kilometers (131 000 feet), the chemistry 
is such that odd-N compounds are split somewhat evenly among the three species, 
sc that individual measurements of one or even two of these species in this 
region do not yield accurate estimates of the total odd-N content of the strat- 
osphere to which a model can be adjusted. This fact is important to the Space 

Shuttle problem since predicted Cl perturbations are quite sensitive to the 
amount of NO present. Various measurements of NO differ by as much as an order 
of magnitude. Estimates of odd-N mixing ratios in the near-asymptotic region 
above 30 kilometers (98 000 feet) made from the existing measurements range 
from a few parts per billion by voliune to a few tens of parts per billion by 

volume. Most models contain asymptotic odd-N mixing ratios of 10 to 20 ppbv. 

More high-altitude NO measurements are necessary, in addition to in situ sim- 
ultaneous measurements of all the major odd-N constituents, to test the chemical 
reaction schemes. 

The major source of odd-N compounds in the stratosphere is N^O, Ground- 
level measurements of N 2 O (R. A. Rasmussen, private communication) yield re- 
sults in the range of 0.3 to 0.35 part per million by volume (ppmv) of N^O. 

These are higher than earlier measurements (ref. ^-8), which were 0.25 to 0.2T 
ppmv. In contrast to previous measurements, the measurements in remote areas 
show almost no random variation. This uniformity possibly indicates a longer 
atmospheric residence time. If the stratosphere is the only sink and if the 
longer residence time is real, a source/sink balance problem, which is cur- 
rently being debated, still remain's. 

New stratospheric profiles of N^O (A. L. Schmeltekopf , private commimica- 
tion) show a much more rapid decrease with altitude than did previous results. 
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This difference may be a manifestation of the variability of atmospheric dy- 
namics . 


Transport Terms 

Effective vertical eddy diffusion coefficients have been derived from the 
analysis of vertical concentration profiles for substances with known loss rates 
in the stratosphere. For example, the vertical distribution of CHj^ is deter- 
mined by the effective upward transport rate and the reaction of CHj^ with O^D, 

OH, and Cl. The loss rate by chemical reaction can be estimated and used with 
the profile data to obtain an effective vertical eddy diffusion coefficient. 

Different species (as ceurbon-ll* ) , seasons, or locations yield different 
eddy diffusion coefficients. This conclusion is illustrated in figure k-2, in 
which six different profiles proposed by different investigators are compared. 
Substantial differences are seen, and these are reflected in results from 
stratospheric models. It is evident that work is needed to reduce or to under- 
stand the uncertainties associated with vertical transport. 


MODELING OF 0^ DEPLETION BY HCl DEPOSITION FROM SPACE SHUTTLE 


Hydrogen chloride deposition by Space Shuttle booster combustion commences 
at launch and continues until burnout at U 5 kilometers (155 000 feet). The 
rate of HCl emission in metric tons per kilometer (ref. h-9) is shown in fig- 
ure I+-3. At peaJk operational activity, 60 launches per year are planned. Launch 
operations are expected to commence in 1978 ; they build up to 60 per year by 
I 98 U according to the sched\ile shown in figiore 4-U. The rate of HCl emission 
from the alternate propellants described in the previous section can be obtained 
by multiplying the baseline emission rate by appropriate factors derived from 
table 3 - 1 I. 

Five groups of investigators, each with a different one-dimensional strat- 
ospheric model, were represented at the workshop. Present were R. C. Whitten 
from the NASA Ames Research Center (ARC), R. Turco from R&D Associates, 

R. Cicerone and S. C. Liu from the University of Michigan, J. Chang from the 
Lawrence Livermore Laboratory (LLL) , and R. Rundel and R. Stolarski from JSC. 

In addition, written and verbal (telephone) information was provided by S. Wofsy 
of Harvard University. These groups were asked to run their models for Space 
Shuttle injection of Cl at the rate of 60 launches per year, with north/south 
transport rates taken crudely into acco\mt by assuming the source to be averaged 
over the Northern Hemisphere. 

Although all the models used the same amount of added Cl, they differed in 
the details of the reaction schemes; in choice of rate constants and cross sec- 
tions, boundary conditions, and transport parameters; and in details of the 
mathematical manipulations. The results are shown in the third column of 
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table U-II. A range of calcixLated 0^ depletions from 0.1 to 0.3 percent for 
the hemi spherically averaged Space Shuttle source is found. 

In the following sections of this report, more detailed comparisons of the 
various models are described. For this purpose, an attempt was made to separate 
the Cl perturbation predictions into effects caused mostly by transport (and 
boundary conditions) and those caused mostly by chemistry. The calculated 
asymptotic (i.e., above approximately 35 kilometers (llU 000 feet)) mixing 
ratio of total Cl (C1 ) is listed in the second column of table H-II and the 

X 

efficiency of the pertxorbing Cl in terms of percentage change of 0 per part 

X j 

per billion by volimie of added Cl^, is shown in the fourth colimin. The former 

depends only on transport, whereas the latter indicates to a large extent dif- 
ferences in assumed chemistry. 


Model Comparison 

Differences in the models used are discussed in the following paragraphs. 

Transport parameters .- Most of the models employed to evaluate the effect 
on the 0^ layer of HCl emitted by Space Shuttle launch vehicles are one dimen- 
sional. Transport is thus confined to the vertical dimension and must be pa- 
rameterized by an "eddy diffusion coefficient" profile. The physical basis 
for such a set of parameters is weak. In practice, one can make use of the 
concept by adjusting the eddy diffusivity profile such that it leads to pre- 
dictions of vertical distributions of tracers such as carbon-1 CHj^, and total 

odd N that are close to observed values. Vertical eddy diffusion, then, rep- 
resents the real transport with all its spatial and temporal variability av- 
eraged over the globe or hemisphere and over a long time period. 

The eddy diffusion profiles used by various groups to study Space Shuttle 
effects on 0^ are shown in figure U-2. To predict measured tracer distributions 

adequately, eddy diffusivity profiles must have the following characteristics. 

1. A rapid decrease from large tropospheric values to a deep minimim 
located between 13 and 15 kilometers (U2 000 and 000 feet) 

3 3 2 

2. A value at the minimxjm ranging from 1 x lO to 5 x 10 cm /sec 

5 2 

3. A rapid increase from the minimum to a value of >10 cm /sec at 
50 kilometers (l63 000 feet) 

Beyond these broad characteristics, the various model eddy diffusivities differ 
in that they predict somewhat different tracer distributions and are dependent 
on model structure. 

It was generally felt that uncertainties were approximately ±2 near the 
minimum and ±3 near 50 kilometers (l63 000 feet). The coefficients referred 
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to in table U-II are shown in figure U-2. The labels in figiire h-2 refer to 
those in the second column of table U-II for each model. Relatively large vol- 
umes of the coefficients result in small values for the steady-state asymptotic 
Cl^ mixing ratios from the Space Shuttle and vice versa. 

Boundary conditions .- It was generally agreed that, within the realm of 
reasonable choices, boundary conditions were not particularly crucial to the 
Space Shuttle problem. The problem of choosing an upper boundary condition for 
CHj^ flux was discussed. This value assumes importance when CHj^ is being used 

to deduce an eddy coefficient. Boundary conditions for the Cl injected by 
Space Shuttle should be a very small or zero mixing ratio at the tropopause 
and zero flux at some altitude above U 5 kilometers (155 000 feet). 


Effect of oxides of hydrogen on Cl perturbations of 0 ^«~ In this and the 


next two subsections, aspects of the chemistry of the models which determine 
the efficiency of 0^ perturbation as given in the fourth column of table ii-II 

are discussed. Some of the differences can be ascribed to the choices made in 
calculating the oxides of hydrogen (HO^). The key constituent is OH because its 

direct reaction with HCl forms Cl, which participates in the catalytic 0^ des- 


truction cycle. Computed values of the concentration of OH at 35 to itO kilo- 
meters (III+ 000 and 131 000 feet) are shown in the seventh and eighth columns 
of table 4-II, The estimates used for two crucial rate parameters, OH + HO^ 


H 2 O + 0^ and 0 + HOg OH + 0^, are given in the fifth and sixth columns. These 

rates are crucial because they exert a controlling influence on the calculated 
OH concentration and because their rates are virtually unknown. One technique 
is to adjust these rates within a limited range to fit the recent OH measure- 

ments (ref. U-IO), ~2xl0 cm at hO kilometers (131 000 feet). This "tuning” 
process was discussed together with the uncertainties caused by the apparent 
variability of OH as indicated in Burnett's colimin content measurements. It 
was agreed that the calculated OH shown in the table fairly closely reflects 
the choice of rates in the fifth and sixth colimins and that the efficiency is 
directly related to the computed OH around 40 kilometers (131 000 feet). 


Effect of No., on Cl perturbations of 0^,- An inverse relationship exists 
between NO^ concentrations and the efficiency of Cl perturbations. This re- 
lationship occurs for two reasons. 

1. The NO catalytic cycle is the principal destruction mechanism for 0_ 

X 3 

in the stratosphere. Leirger NO concentrations produce a larger 0_ destruction 

X 3 

rate; hence, the fractional effect of a given Cl^ concentration is less. 

2. The CIO reacts with both 0 and NO, but the latter leads to no net 0^ 
destruction; hence, increased NO leads to a lower efficiency of 0^ destruction. 
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Colimin 10 of table U-II shows the calculated asymptotic mixing ratio of total 
odd N (no + NO^ + HNOg) for each of the models. Two factors exert the major 

control over the differences. One is the eddy diffusion coefficient and the 

other is the choice of O^D reaction rates. The effect of the eddy diffusion 
coefficient is such that a larger coefficient tends to diffuse more N_0 to al- 

1 ^ 

titudes where it can react with 0 D to produce NO; thus, a larger source is 
provided while the produced NO is diffused more rapidly toward its sink in the 
troposphere. The net effect is calculated to be in the same direction as that 
for the Space-Shuttle-injected Cl^; that is, faster diffusion implies less odd 

N. However, because of the competing N^O diffusion, the odd-N concentration 

is not as sensitive to the eddy diffusion coefficient as is the Cl concentra- 

X 

tion. This fact is clearly illustrated in columns 2 and 10 of table U-II. 

It was generally agreed that the gross features of the differences in Space 
Shuttle perturbation calculations were explained on the basis of the differ- 
ences in the eddy coefficients and the choice of reaction rates including OH 

and O^D. The OH and the odd-N contents reflect these differences. 


Chlorine nitrate formation .- It has been suggested that the formation of 
CINO^ could be important to the stratospheric chemistry of Cl. Its formation 

by way of a three-body reaction (CIO + NO 2 + M CINO^ + M), its destruction by 

way of photodissociation, and its reaction with 0 or OH were hypothesized. If 
this suggestion is correct, the net effect is to tie up both Cl and NO in a 

temporary nonreactive compound and to reduce the expected 0^ perturbation. 


At the stratospheric workshop, none of the model calculations reported 
included possible effects caused by CINO^; no generally accepted measurement 

of its formation rate has been made. However, calculations made after the 
workshop by the JSC and University of Michigan groups indicate that if the 

31 6 

formation rate of CINO^ is approximately 10~ cm /sec, reported 0^ depletions 


are reduced by a factor of approximately 2. The ARC group calculated a factor 
of 2 to it (appendix D). 


Application of Models to the Space Shuttle Program 

The various models were applied to questions related to the Space Shuttle, 
as shown in the following paragraphs. 

Evaluation of suggested alternate propellants .- Three of the models were 
used to evaluate the change in the estimated 0^ depletion which would resiilt 

from a switch from the baseline propellants to the alternates "A" and "B" de- 
scribed earlier in this report. The improvement factor for each alternate 
calculated by the various models is shown in table i+-III. There is general 
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agreement that a switch to alternate "A" at 20 kilometers (63 000 feet) will 
decrease the predicted 0^ depletion by a factor of U to 5 if alternate "A" has 

3.1 percent HCl at the exit plane and a factor of 3 if it has 5.8 percent HCl 
at the exit plane. For alternate ”B," the improvement is slightly less, a 
factor of U for a 3.1 percentage and a factor of 2.5 for a 5.8 percentage. 
Alternate "B" is slightly worse than alternate "A" because the decrease in burn 
rate is more than compensated for by the increase in the altitude of burnout. 

A parametric study was conducted for the alternate propellants as a fxinc- 
tion of switchover altitude. These results, shown in figiire *+-5> are not sub- 
stantially different from the previous studies by JSC and ARC in which a switch- 
over to a completely clean propellant was studied. The detailed shape of this 
curve will vary with exactly which eddy coefficient is chosen. Results indicate 
that increasing the altitude of switchover to 25 kilometers (82 000 feet) would 
not greatly change the 0 ^ depletion effects. 

Corridor effect .- The latitude vsiriation of 0^ reduction was approximated 

with the one-dimensional models by assuming the exhaust products to be distri- 
buted uniformly over the Northern Hemisphere only. A more exact calculation 
of the latitude variation of 0^ reduction was made by using the Massachusetts 

Institute of Technology three-dimensional model and the ARC two-dimensional 
models , 

The largest effect is during the Northern Hemisphere summer, when the de- 
pletion near latitude 30° N is approximately a factor of 2 larger than in the 
remainder of the Northern Hemisphere. The Northern Hemisphere depletion is a 
factor of 2 to 3 larger than in the Southern Hemisphere. During the spring, 
the latitude 30° N to Northern Hemisphere average ratio is approximately 1.5, 
whereas the Northern Hemisphere average to Southern Hemisphere average is 
approximately 2 . 

If these ratios are applied to the one-dimensional hemispheric ally averaged 
calculations, it would mean that 20 to 30 percent of the 0^ depletion will occur 

in the Southern Hemisphere and that the Northern Hemisphere depletion can be 
correspondingly decreased by this amoimt. This resulting Northern Hemisphere 
value will have a factor of 'ui .5 enhancement in the simnner in a corridor aro\ind 
latitude 30° N and a factor of 'uO.5 reduction in the spring. These remarks are 
illustrated in figure U- 6 , in which the summer and winter latitude distribu- 
tions are shown. The curves labeled "summer" and "winter" have been adjusted 
so that the area under them is eqxial to the area tinder the broken line labeled 
"one-dimensional approximation." 

Time dependence .- A time- depen dent calculation of the buildup and decay of 
0^ perturbations caused by Space Shuttle launches is shown in figure U-T. The 

delay in buildup reflects the slow buildup in the launch schedule. The time 
constant for ozone recovery if operations should cease may be as short as 2 
years and as long as 6 years, depending on the choice of eddy diffusion 
parameters. 
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TABLE U-I.- RECENT MEASUREMENTS OF CRITICAL RATE CONSTANTS^ 


Reaction 

Rate expression. 

Reference 



cm /molecule-sec 


HCl + 0 -*■ 

Cl + OH 

2.5x10"^^ exp(-2900/T^) 

4-1 

HCl + OH 

* Cl + H^O 

2,0x10"^^ exp(-313/T) 

4-2 



4.1x10“^^ exp(-500/T) 

4-3 

CHj^ + Cl 

+ CH^ + HCl 

5.1x10"^^ exp(-1790/T) 
7.44x10“^^ exp(-1226/T) 
( 1.2 ± 0.3)xl0~^^{295 K) 

4-4 

H 2 + Cl ^ 

H + HCl 

4. 7x10"^^ exp(-2340/T) 

— 

Cl + 0^ ^ 

CIO + O 2 

2.7x10“^^ exp(-257/T) 
(1.3 ± 0.3)xl0‘^^(295 K) 

— 

CIO + 0 

Cl + O 2 

1.07x10“^° exp(-224/T) 
3.36x10"^^ exp(+76/T) 

— 

CIO + NO 

+ Cl + NO^ 

1 . 9 x 10 "^^ 

6 . 3 x 10 “^^ exp(+335/T) 

— 


(Cl + 20„)| 

8x10"^^ (300 K) 



2 1 


CIO + 0 

1 

->■ ClOO + O 2 1 
- C 102 + 02 ) 

1.3xl0~^®(300 K) 


0(^D) + H 

3 O products 

(2.1 ± 1.0)x10”^°(298 K) 

4-5 

0 (^D) + CHj^ -*■ products 

(1.3 ± 0.3)x10"^°(298 K) 

4-5 

0(^D) + H, 
< 

5 products 

(1.3 ± 0.05)x10"^°(298 K) 

4-5 

0(^D) + N, 
< 

, products 

(3.0 ± 0.01)x10”^^(298 K) 

4-5 


°This table contains meeisxired rates of important stratospheric reactions 
which have appeared subsequent to the Climate Impact Assessment Program (CIAP) 
table of recommended rates. 

= ambient temperature. 
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TABLE U-I.- (Concluded) 


Reaction 

Rate expression, 

3 

cm /molecule-sec 

Reference 

0(^D) + 0^ products 

(»*.l ± 0<,05)><10"^^(298 K) 

k-5 

0(^D) + N^O -*■ products 

(l.U ± 0.1 )x10"^°(298 K) 

h-5 

Cl + HCl + HO^ 

(6.2 ± 1.5)><10"^^(295 K) 

— 

Cl + HOg -*■ HCl + 0^ 

(3 ^ J]^)xl0"**(295 K) 

— 

Cl + HNO^ products 

( 6.8 + 3.it)xl0"^^(295 K) 

— 

CINO^ + hv ->■ products 

— 

— 
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M TABLE COMPARISON OF 0^ DEPLETION EFFECTS PREDICTED 

BY VARIOUS STRATOSPHERIC MODELS 


|BBd 

ftad 

Modeling 

Space 

Shuttle 

asymptotic 

Cl 

X, ppbv 
(a) 

Percentage 0^ 

Efficiency, 

JfAO^ 

X 

'^OH + HO 2 , 

’'O + HO^, 

Diurnal 

average 

Sources for 

Asymptotic 
odd N, 
ppbv 

CH|^(U0), 

group 

depletion, 
60 launches/ 

yr 

cm /sec 

cm^/sec 

OH{35), 

-3 

cm 

OH(UO), 
cm ^ 

O^D 

rates 

ppbv 

m 

NASA ARC 

0.056 (B) 

0.11 

2 

6 x 10 “^^ 

6 x 10 “^^ 

7 x 10 ^ 

1.6x10^ 

CIAP'^ 

15 

390 


Univ. of 
Michigan 

.15 (C) 

.29 

2 

2x10"^^ 

6 x 10 “^^ 

7 . 6 x 10 ^ 

2.3x10^ 

York/ 

NOAA*^ 

2h 

300 

pa 

LLL 

.1 (D) 

.26 

2.6 

2x10"^^ 

1.5x10"^^ 

6x10^ 

l.Pxlo"^ 

York/ 

NOAA*^ 

16 

360 

pa 

NASA JSC, 
GSFC 

.16 (E) 

.2 

1.3 

2x10“^^ 

.3x10"^^ 

6x10^ 

1 . 6 x 10 ’^ 

York 

noaa'^ 

2 l* 

260 


Harvard 

Univ. 

.08 (F) 

.2 

2.5 

— 

-- 

— 

~ 2 xlo’^ 

__ 

11 

— 




^Vertical eddy diffusion coefficients from figure U-2 are designated by letters in parentheses. 
'^Climate Impact Assessment Program. 

'^The York/National Oceanic and Atmospheric Administration (NOAA) study was withdrawn. 
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TABLE 1+-III.- IMPROVEMENT FACTORS FOR 0^ DEPLETION RESULTING 
FROM USE OF ALTERNATE PROPELLANTS 


Modeling 

group 

Alternate "A" 

Alternate "B" 

3.1 percent 
HCl 

S.8 percent 
HCl 

3.1 percent 
HCl 

5.8 percent 
HCl 

ARC 

1+ 

— 

h 

— 

LLL 

Uo6 

2.8 

k.l 

2.5 

JSC 

5 

3.3 

h 

2.5 
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Figure U-5.- Ozone reduction as a function of 
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Figure U-T.- Time dependence of ozone depletion from Space Shuttle operations. 
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5. AFTE31BURNING EFFECTS 


Frozen chemical equilibrivnn in the rocket nozzle and interaction of the 
hot exhaust plume with the atmosphere result in the formation of NO. This 
conqjovnad also catalyzes the decon 5 )Osition of The effect of NO from after- 

biiming is assessed in the following peuragraphs. 

The NO produced by the Space Shuttle in the stratosphere has been calcu- 
lated by using a model that includes the effect of the shock waves generated 
by pl\ane interactions and afterburning. These results are given in table 5-1 
for the baseline AP propellant. The model is relatively new, and the figures 
quoted are from the initial rvm of the code. The calculation procedure is de- 

2 

scribed in detail in references 5-1 and 5-2. It should be noted that the 
chemistry in the stratosphere plume model is the same as that used for the 
troposphere plume and ground cloud. Aircraft measurements of the NO^ in the 

Titan groimd cloud agree with the calculated values within 30 percent as re- 
ported in reference 5-3. The Titan propellsuit and the Space Shuttle baseline 
propellant are identical in composition. 

The procedure used to evaluate the NO^ deposition rate from the alternate 

propellant was to establish the significant differences between the two pro- 
pellants and to assess their influence on the mass of NO^ leaving the nozzle 

and the enhancement in the plume. The procedure vised to obtain the values 
shown in table 5-II is a simplification of, and more conservative than, that 
used to generate the results of table 5-1. It is also given in references 5-1 

2 

and 5-2 and an unpublished article. 

The differences between the total amount of NO^ deposited by the baseline 
propellant and by the alternate ”A” propellants are due to the following. 

1. The lower chamber temperature for the alternate "A" propellant allows 
the NO mass fraction to reach its equilibrivnn vEd.ue within the chamber. This 
value is approximately a factor of 50 less than the calculated NO mass fraction 
in the combustion chamber for the baseline propellant. 

2. The alternate "A" propellant exhaust products contain more unburned 
fuel (Hg +00), which results in more intense afterburning. 

3. The shock-induced NO^ enhancement in the plvmie shovild be similar for 
both types of propellant. 


D, Thorpe, H. S. Pergament, and B. Hwang: NO^ Deposition in the 

Stratosphere by the Space Shuttle Solid Rocket Motors. Proceedings 9th JANNAF 
Plume Technology Meeting, Feb. 1976, in press. 
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The total rate of NO^ deposited by the baseline and alternate propellants 

at 15 and 30 kilometers (i*T 000 and 95 000 feet) is compared with the HCl de- 
position rate in table 5-IH. It is notable that the NO^ deposited in the 

stratosphere is not significant relative to the HCl for either propellant. 
Assuming approximately equivalent chemical efficiencies for NO^ and 01^, the 

effect of NO deposition on 0^ depletion is of the order of 0.005 times the 

X j 

effect of Cl^ deposition, or less than 0.001 percent 0^ depletion. 


REFERENCES 


5-1. Pergament, H. S. ; and Thorpe, R. D. : NO^ Deposited in the Stratosphere 

by the Space Shuttle. NASA CR-132715> 1975. 

5 - 2 . Pergament, H. S.; Thorpe, R. D. ; and Hwang, B.: NO^ Deposited in the 

Stratosphere by the Space Shuttle Solid Rocket Motors. NASA CR-litU928, 
1975. 

5-3. Stewart, R. B. ; and Gtomberg, R. I.: The Production of Nitric Oxide in 

the Troposphere as a Result of Solid Rocket Motor Afterburning. NASA 
TN D-8137, 1976. 
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TABLE 5-1.- BASELINE PROPELLANT NO DEPOSITION 

X 


[Values in grams per second] 


Parameter 

Altitude 

15 km 

30 km 

Nozzle 

P 500 

2500 

Plume enhancement 

9 500 

1500 

Total 

12 000 

i4000 


TABLE 5-1 I.- BASELINE AND ALTERNATE "A" PROPELLANT 

NO DEPOSITION 

X 

[Values in grams per second] 


Parameter 

Altitude 


15 km 

30 km 

Nozzle 



Baseline 

2 500 

2500 

Alternate "A" 

55 

55 

Plume enhancement 



Baseline 

9 500 

1500 

Alternate "A" 

l6 000 

3000 

Totals 



Baseline 

12 000 

1*000 

Alternate "A" 

16 055 

3055 
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TABLE 5 -III.- COMPARISON BETWEEN NO AND HCl DEPOSITION RATES 

X 


Deposition rate 

Altitude 

15 km^ 

30 km^ 

No^ deposition, kg/km (ton/km) 



Baseline 

20.865 (0.023) 

3.629 ( 0 . 00 i 4 ) 

Alternate "A" 

28.123 (0.031) 

2.722 ( 0 . 003 ) 

NO^/HCl deposition ratio 



Baseline 

0.01 

0.005 

Alternate "A" 

0.11 

0.027 


locity = 561 m/sec. 
^Velocity = lli +0 m/sec. 
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6. EFFECTS OF PARTICULATE EMISSION FROM 


THE SPACE SHUTTLE EXHAUST 


Concern has been expressed that the particulates formed from the Al^O^ 

emission of the Space Shuttle would remain suspended in the stratosphere and 
cause climatic changes in the troposphere. It has also been postulated that 
these particles would increase the 0^ destruction rate. To investigate these 

concerns, studies were carried out which included the analysis of particles 
from the Titan rocket exhaust, the measurement of 0^ surface reaction rates 

on AlgO^i and analysis of these results to predict the environmental conse- 
quences of Al 20 ^ dust in the stratosphere. 


RADIATION SCATTERING BY THE EXHAUST PARTICLES 


Estimates of the Earth mean temperature can be obtained by equating the 
amount of solar energy absorbed by the entire Earth with the amount of thermal 
radiation it emits to space (ref. 6-1). By performing such calculations for 
perturbed and \inperturbed conditions, one can estimate the change in mean sur- 
face temperatiure produced as a resiilt of the injection of Al^O- particles by 
the Space Shuttle vehicle. 

Mie scattering theory was used to relate the optical depth perturbation 
caused by the Space Shuttle to the steady-state mass loading of the added par- 
ticles (ref. 6-1). These added particles are found from the mass of Al^O^ gen- 
erated by a single Space Shuttle vehicle, the projected number of Space Shuttle 
flights, and the mean dwell time of particles in the stratosphere. For a pro- 
jected traffic level of one Space Shuttle flight per week, an optical depth 

perturbation of 3.3 ^ 10 an average albedo change of it x 10~^, and a temper- 

-It 

ature decrease of 1.5 10 K for the Northern Hemisphere were found. The 

corresponding ntmibers for the Southern Hemisphere are 0.it3 as large. These 
temperature changes are far smaller than those associated with significant ob- 
served climatic changes (ref. 6-2). 

3 

Kessler studied the effect of Space Shuttle exhaust particles on the 

Earth's albedo and obtained an upper limit for the effect of 3 ^ lO”^. From 
these results, it can be concluded that the uncertainty in the previous res\ilts, 
are approximately a factor of U due mainly to \ancertainties in particle residence 
times in the stratosphere, particle scattering theory, and particle size 
distribution. 


u. J. Kessler: Current Assessment of the Effect of the Altuninum Compound 

Emission From the Space Shuttle. Rep. JSC-09913, Aug. 1975. 
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The preceding analyses are very sensitive to the number of particles in 
the size range between 0.05 and 1 micrometer. Although particles of 0.01- 
micrometer radii are too small to interact efficiently with light, they could 
serve as nucleation centers for H^SOj^ and grow into larger H^SOj^ droplets. It 

was found (ref. 6-3) that the nimiber of Space-Shuttle-injected particles of 
this radius (and larger) would exceed or equal the number of particles of this 
size naturally fo\md in the stratosphere. The growth of these particles has 
not yet been studied in detail, but any increase in the stratosphere sulfate 
layer could be more significant than the introduction of Al^O^ particles alone. 

Until the nucleation and growth is understood, the climatic effects cannot be 
easily estimated. Further studies are being performed at ARC. 


DESTRUCTION OF 0^ BY SURFACE REACTIONS 


Studies conducted at the NASA Jet Propulsion Laboratory (L. F. Keyser) and 
at ARC (T. Fujiwara) provide strong evidence that dispersal of Al^O^ in the 


stratosphere by the Space Shuttle rocket engines will not result in a signifi- 
cant depletion of 0^. The essential findings summarized in table 6-1 are that 

^ -9 

the collision efficiency for 0^ destruction is 2 x iQ or less at room temper- 


ature and decreases with temperature. These values are approximately three 
orders of magnitude below the threshold of detectable long-range 0^ destruction 

rates (O.l percent of the normal 0^ loss rate). The experimental studies fur- 


ther showed that there was no significant enhancement of catalytic activity by 
exposure to HCl, as might occur in the exhaust pltune. Outgassing at temperatures 
between 600 and 1000 K activates A1^0_ for 0_ destruction. Collision efficien- 

cies following outgassing are estimated to be as high as 10~ to 10~ for A1„0_ 

-7 ^ 

and 10 for samples of alumina from the solid rocket motor exhaust. Even at 


-7 

10 , the destruction rate of stratospheric 0^ is a factor of 10 below the de- 

tectable threshold for 0^ destruction. Furthermore, exposure to air or 0 de- 


activates the alumina in less than 2k hours. The current assessment is that 
the long-term effect of alumina exhaust particles on the 0^ layer is negligible. 


REFERENCES 


6-1. Pollack, J. B. ; Toon, 0. B. ; et al,: Estimates of the Climatic Effect of 

Aerosols Produced by Space Shuttles, SST's, and other High-Flying Air- 
craft. J. Appl. Meteorol. , vol. 15, no. 3, Mar. 1976, pp. 247-258. 

6-2. Pollack, James B.; Toon, Owen B. ; et al.: Volcanic Aerosols and Climatic 
Change: A Theoretical Assessment. J. Geophys. Res., vol. 8l, no. 6, 

Feb. 20, 1976, pp. 1071-1083. 
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tribution of the Space Shuttle Effluent to the Natural Stratospheric 
Aerosol. Geophys. Res. Letters, vol. 2, no. 3, Mar. 1975, pp. 113-116. 
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TABLE 6-1.- OZONE DECAY RATES AND COLLISION EFFICIENCIES ON Al^O^ 


Designation 

Sample 

« a 

Form 

Crystallite 

0 

size, urn (A) 

Specific surface 
2 , 

area, m /g 

Collision 

efficiency 

A 

Laboratory preparation 

Gamma 

0.0050 

(50) 

123 

7.8x10"^° 

B 

Degussa aluminum 
oxide-C 

Gamma 

.0065 

(65) 

116 

1.6x10“^° 

C 

Alcoa F-1 activated 
alumina 

Boehmite, major; 
Gamma, minor 

.0050 

(50) 

208 

-9 

l.lxlO ^ 

D 

BaJier reagent powder 

Alpha, major; 
Theta, trace 

— 


9.1 

2.8x10“^° 

SRM-A^ 

JPL SRM inside tank 

Alpha, major; 
Gamma, Delta 
minor 



k.d 

l+.0xl0"^° 

SRM-B 

JPL open burn inside 
tank 

Theta, major; 
Gamma, Delta, 
Alpha, minor 



5.1+ 

1.9x10"^° 

SRM-C^ 

■ 

JPL SRM inside tank 

Alpha, major; 

Gamma, Delta, 

minor 
— — 

— 


1.8 

l.TxlO"^ 


Analyses performed Ly Technology of Materials, Santa Barbara, California. 

^Gaseous residence time in the solid rocket motor (SRM) was approximately 15 milliseconds 
for SRM-A and approximately 90 milliseconds for SRM-C. Motor SRM-C more closely approaches 
conditions expected for Space Shuttle solid booster rockets. 







7. BIOSPHERIC IMPACT OF uv RADIATION INCREASE 


RESULTING FROM 0^ DEPLETION 

The estimated average percentage of 0^ depletion presented in section 

can be converted by a series of assumptions and calculations into an estimate 
of the increase in BHuv radiation from the S\m, Most estimates indicate that 
a 1-percent 0^ reduction results in a 2-percent increase in BHuv radiation. 

Considering the uncertainty in the estimate of 0^ depletion, the attendees of 

the workshop adopted a 0.2-percent 0^ depletion and a 0.4-percent increase in 

BHuv radiation as a worst-case situation for an assessment of the biospheric 
impact of the projected Space Shuttle operations outlined in the previous 
section. 

Based on the limited available biological data, the impact on the biosphere 
of an assimied 0.4-percent increase in BHuv radiation was considered (appendix 
H), and the following conclusions were drawn. 

1. The effects on the biosphere of a 0.4-percent increase in BHuv radia- 
tion will not be detectable with decades of observation because the natural 
fluctuations of BHuv radiation and the ecosystem are so much larger than this 
increase. These fluctuations include (l) greater variability of BHuv irra- 
diances; (2) statistical imcertainty involved in the responses of organisms 
and ecosystems to given doses of BHuv radiation at given dose rates; and (3) 
large normal variation in the response of organisms to other, more pervasive 
or important environmental factors to which they are exposed such as tempera- 
ture, moisttire, nutrition, competition, and predation. The fact that a cause- 
and-effect relationship of an increase in BHuv radiation cannot be statistically 
detected, however, does not mean that an increase in BHuv radiation can be 
ruled out as a contributor to some deleterious futijire event. Lack of detecta- 
bility is not equatable with a lack of effect. 

2. There is a very low probability that an average 0.4-percent increase 
in BHuv radiation will have any unacceptable effect on agricultural plants or 
natural ecosystems, independently of whether the effects of such an increase 
are detectable or not. The basis for this conclusion is outlined in appendix H. 

3. There may be some increase in the number of melanoma and nonmelanoma 
skin cancer cases among susceptible individuals resulting from 0.2-percent 0^ 

depletion, but it will not be detectable because of (l) the natiiral variations 
of BHuv radiation and biological variability mentioned previously, (2) the long 
latent period (20 to 60 years) for induction of skin tumors, and (3) the many 
other factors that already may be tending to increase or decrease the number 
of reported skin cancer cases. Factors possibly leading to an increased num- 
ber of melanoma cases in the United States include (l) the increased proportion 
of people in the population living long enough to contract skin cancer; (2) the 
increased reporting of skin cancer cases because of Medicare; (3) the changing 
of lifestyles, which, in recent years, involves more leisure time activity in 
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the sunshine; and (U) the net southward migration of the population. (Regarding 
the latter factor, between 19^0 and 1970, the center of population of the United 
States moved west and approximately 35 miles south from approximately 39° N 
latitude; this amount of southward movement corresponds to a 2-percent increase 
in the annual dose of BHuv radiation. ) Factors possibly leading to a decreased 
incidence inc-lude action based on publicity-induced recognition by the popula- 
tion of the dangers of overexposure to solar radiation and on more accurate 
identification of susceptible individuals as a result of research generated by 
the 0^ depletion problem. 

The number of skin cancer cases resulting from a 0.2-percent 0^ reduc- 
tion is not realistically predictable. The extant experimental and epidemiolo- 
gical data, althoiagh clearly suggesting some contribution of solar radiation 
exposure to skin cancer incidence, are inadequate for making quantitative pre- 
dictions with reasonable limits of uncertainty. (See appendix H. ) 
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APPENDIX A 


ALTERNATE BOOSTER PROPELLANT STUDIES 
By J. Q. Miller^ 


The alternate propellant study was initiated in August 197^ to define a 
solid-propulsion supporting development program for eliminating or significant- 
ly reducing the release of gaseous hydrogen chloride (HCl) in the upper atmos- 
phere while minimizing Space Shuttle Program cost, performance, and schedule 
impacts. Activities included propellant formulation and characterization and 
solid rocket motor design development. The major gro\md rules for this study 
were as follows. 

1. Conduct the first operational flight in 1983. 

2. Eliminate or minimize HCl release above 20 kilometers (63 000 feet). 

3. Consider class II propellants only. 

h. Retain 1-lA configuration envelope (baseline). 

Approximately 250 different propellant formulations containing various 
combinations of ammonium perchlorate (AP) and non-HCl-producing oxidizers such 
as ammonium nitrate, cyclotetramethylenetrinitramine (HMX), trimethyletha- 
nitrinitrate, nitrocellixlose, nitroglycerine, and cyclotrimethylenetri- 
nitramine were evaluated theoretically. These formulations were compared with 
the baseline AP-oxidized propellant. From this comparison, the most promising 
alternate propellant formulations were selected for preliminary process, com- 
patibility, hazard, burn rate, pressure exponent, and density impulse charac- 
terization. Based on results from these characterizations, the best HMX- and 
non-HMX-containing alternate propellant formulations were selected for addi- 
tional evaluation. A formulation and ballistic property comparison of these 
two alternate formulations ("A" and ”B") with the baseline formulation is 
presented in table 3-1. To achieve acceptable ballistic properties, both 
alternate formulations require the addition of AP. The theoretical exhaust 
composition for the baseline and alternates "A" and "B" is presented in 
table 3-1 I. 

The most promising design concept developed was a double-web stgir 
propellant grain configuration (dual-grain concept). With this concept, the 
baseline propellant is completely burned out at an altitude of approximately 
20 kilometers (63 000 feet) and the alternate propellant (low HCl) is burned 
at altitudes above 20 kilometers (63 000 feet); thus, minimum cost and perform- 
ance impacts are afforded. Alternate propellant formulations "A" and "B" 
were used in conjimction with the baseline propellant to design two booster 
solid rocket. motors incorporating the dual-grain concept. Figures 3-1, 3-2, 


^ASA George C. Marshall Space Flight Center. 
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and 3-3 contain thrust as a f\inction of time and mass flow characteristics for 
the baseline, alternate "A," and alternate "B" designs, respectively. By 
comparing mass flow characteristics with exhaust compositions below and above, 
the transition altitudes (change from baseline to alternate propellant burn), 
the predicted exhaust constituent weights can bf- establisiied. 

Results from preliminary burn rfoe tests conducted on the alternate "B" 
propellant indicate that AT levels as high as PO percent by weight may be 
required to obtain acceptable burn rates and combustion efficiencies. This 
increase in AP level will increase the ]IC], ])erc'-ntage weight in the exhaust 
from 3.1 to approximately 5.8. Other exhau.st constituent percentage weights 
are not significantly affected by this change. Therefore, for environmental 
assessments at this time, it is recommenneeJ that an HCl percentage weight range 
of .^.1 to 5.8 be used for the alternate propellant. 

The alternate booster propellant studies are planned for continuation 
throiigh December 1976. The data contained here concerning exhaust composi- 
tion and mass flow characteristics are not expected to change significantly. 
Most of the data generated between now and December 1976 will be used to 
update Space Shuttle Program cost, payload, and schedule impacts. 
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APPENDIX B 

STRATOSPHERIC REACTION RATES 
By W. DeMore^ and L. F. Keyser^ 


Current measurements of selected critical rate constants used for Space 
Shuttle exhaust modeling are svimmarized in table ^i-I. Additional data are 
contained in reference U-1. A recent review of rate data for oxides of chlorine 
has been prepared (ref. k-2) . 


Studies conducted at the NASA Jet Propulsion Laboratory and the NASA 
Ames Research Center provide strong evidence that dispersal of aluminum oxide 
in the stratosphere by the Space Shuttle rocket engines will not result 

The essential finding is that the 


in significant ozone (O^) depletion. 


collisional efficiency for 


0^ destruction is in the range of 10 


-9 


to 10 


-10 


at room temperature and decreases with increased temperature. These values 
are approximately three orders of magnitude below the threshold of detectable 
long-range 0^ destruction rates (O.l percent of the normal 0^ loss rate). 


The experimental studies showed further that there was no significant en- 
hancement of catalytic activity by exposure to hydrogen chloride, as might 
occur in the exhaust plume. 


®NASA Jet Propulsion Laboratory. 
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APPENDIX C 


MEASUREMENTS OF STRATOSPHERIC MINOR CONSTITUENTS 
By the Environmental Effects Project Office^ 


CHLORINE COMPOUNDS 


Personnel at the National Center for Atmospheric Research now have 
three sets of hydrogen chloride (HCl) measurements obtained using the filter 
technique. The measurements made by balloons reveal a concentration ranging 
from less than 0.1 part per billion by volume (ppbv) near 15 kilometers 
(Ut 000 feet) to 0.5 to 1.0 ppbv at 25 to 30 kilometers (T9 000 to 98 000 feet). 
These measurements clearly indicate a high-altitude (approximately 30 kilo- 
meters) (98 000 feet)) soiirce and a lower altitude sink (almost certainly 
rainout in the troposphere). Remote infrared spectroscopic measurements 
(refs. C-1 and C-2) show the same general shape and also indicate an apparent 
approach to an asymptotic mixing ratio above 30 kilometers (98 000 feet). 

These measxirements are all shown in figure C-1, together with a calculated 
profile (ref. C-3) based on a specific set of tropospheric halogen sources. 

The Lazrus et al. measurements (ref. C-4) indicate an asymptotic mixing 
ratio in the range of 0.5 to 1.0 ppbv. Latitudinal measurements were also 
made at 15 and I 8 kilometers (^7 000 and 58 000 feet) from an airplane. As 
expected for a tracerlike species, these show poleward increase (fig. C-2). 

Whereas measurements of HCl give a good indication of the total chlorine- 
containing molecules (Cl^) and information about sources, sinks, and transport, 

a measiirement of chlorine oxide (CIO) would tend to indicate the activity of 
Cl toward ozone (0_) in the stratosphere. Resialts are thus far inconclusive. 

X j 

Waters (ref. C-5) has looked for CIO with his microwave spectrometer and can 
set an upper limit of 2 to 3 ppbv, more than an order of magnitude over pre- 
dictions. Carlson (private communication) has a tentative detection of 

lU —2 

column content by way of solar absorption. The value (~9 ^ 10 cm” ) is 
much larger than expected, but the error bars still encompass predicted values. 
Anderson (private communication) conducted his in situ measurement in March 1976. 


HYDROXYL MEASUREMENTS 


An important measure of the odd-hydrogen (H) content of the atmosphere is 
hydroxyl (OH), which can catalytically destroy 0^ in the stratosphere or which 

can alternatively tie up oxides of nitrogen (NO ) as nitric acid (HN0_). 

X j 

®NASA Lyndon B. Johnson Space Center. 
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Previously, only Anderson's measurements at altitudes from U6 to 70 kilometers 
(li+T 000 to 237 000 feet) existed. Recently, there have been measurements at 
ground level (ref. C-6) , at 7 and 12 kilometers (22 000 and 36 000 feet) 

(ref. C-7), and from 30 to h3 kilometers (95 000 to 137 000 feet) (ref. C-8) 
as shown in figure C-3. The lower altitude measurements are of great interest 
in dealing with such problems as the oxidation of sulfur dioxide (SO^) to 

sulfuric acid (H^SOj^) and the destruction of hydrogen-containing chlorocarbons 

in the troposphere. The stratospheric OH measurements are of direct use in 
the evaluation of Space Shuttle Cl injections. They have already been used 
in an attempt to tune a model for use in Cl perturbation predictions (ref. C-8). 

Further information is provided by Burnett's ground-based OH column con- 
tent measurements (ref. C-9). Burnett frequently sees column contents of 
ik -2 

several times 10 cm . However, the OH concentration seems to be patchy; 

13 

and an upper limit can frequently be set to the column content at 2 x 10 to 
13 

3 X 10 , which is smaller than that implied in the in situ measurements. The 

situation is still uncertain but is many times improved by recent measurements. 


MEASUREMENTS FOR OXIDES OF NITROGEN 


The measiirement situation for NO^ constituents is quite confused. There 

are numerous measurements at different times and places of nitric oxide (NO), 
nitrogen peroxide (NO^), and HNO^j which have been summarized in reference 

C-2. (See figs. C-U to C-6.) Odd N (NO + NO^ + HNO^) is produced from 

nitrous oxide (NgO) by reaction with 0^, where ^ represents an energy state. 

It is transported eiround the stratosphere and eventually back to the troposphere 
for loss by rainout. Below an altitude of approximately UO kilometers 
(127 000 feet), the chemistry is such that odd N is split somewhat evenly among 

the three species. Individual measiirements of one or two of these species in 

this region do not yield accurate estimates of the total odd-N content of the 
stratosphere to which a model can be adjusted. Predicted Cl perturbations are 

quite sensitive to the amount of NO present in the model. Measurements of NO 

range over about an order of magnitude and represent a combination of trans- 
port and chemical processes. Estimates of odd-N mixing ratios in the asymp- 
totic region at an altitude above 30 kilometers (95 000 feet) made from the 
existing measiirements range from a few parts per billion by volume to a few 
tens of parts per billion by vol\jme. 

Most models contain a combination of transport coefficients and O^D rate 
rate coefficients to yield asymptotic odd-N mixing ratios of 10 to 20 ppbv. 

The most direct measurement of odd N is the NO measurement by C. J. Mason and 
J. J. Horvath (private communication) at an altitude between kO and 60 kilo- 
meters (127 000 and I85 000 feet), where chemistry predicts odd N to be 
virtually all NO — only about 5 ppbv. More high-altitude NO measurements are 
necessary in addition to in situ simultaneous measurements of all major odd-N 
constituents to test chemistry. 
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NITROUS OXIDE MEASUREMENTS 


R. A. Rasmussen’s ground-level measurements (private communication) yield 
results in the rsinge of 0.3 to 0.35 part per million by volume (ppmv) of 
N^O. These values are higher than the earlier measurements (refs. C-10 and 

C-ll), which are approximately 0.25 to 0.27 ppmv. In contrast to previous 
results, measurements in remote areas (ref. C-12) show almost no random varia- 
tion. This uniformity possibly indicates a longer atmospheric residence time. 
If the stratosphere is the only sink and the longer residence time is real, 
a source/sink balance problem, which is currently being debated, remains. 

New stratospheric profiles of N 2 O by A. Schmeltekopf (private communi- 
cation) show a much more rapid decrease with altitude than previous results 
(ref. C-I 2 ), as shown in figure C-7. These profiles are probably a mani- 
festation of the variability of atmospheric dynamics. On the same flight, 
two fluorocarbons (Freon -11 (CFCl^) and Freon-12 (CF^Cl^)) were measured; both 

showed a sudden drop at an altitude of I 8 kilometers (58 000 feet) as shown in 
figures C -8 and C-9. These are either due to some unknown dynamical effect 
or are a warning of the difficulty involved in these types of measurements. 
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Figure C-1,- Remote infrared spectroscopic measurements vith a calculated 

profile. 
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C-2.- Latitudinal measurements from an aircraft at altitudes of 15 and 
l8 kilometers (4T 000 and 59 000 feet). 
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Figure C-3.- Ground-based measurements of OH concentration, where x is solar 

zenith angle. 
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Figure C-k,- Height profiles of meas\ired nitric oxide number densities in the 

chemosphere. 
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Figure C-5.- Height profiles of NOg number densities measured by ground-based 
aircraft and balloon-borne instruments. Constant volume mixing ratios are 
roughly indicated by the straight lines marked 10~®, 10“^, and 10~^^. 
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Figure C-6,- Envelopes representing the measured height profiles of nitric 
acid n\mber densities. The envelope designated L-G corresponds to the data 
of Lazrus and Gandrud (19T^), M to those of Murcray et al, (l97^), G to those 
of Fontanella et al. (197^+), W to those of Fried and Weinman (l970), and H 
to those of Harries et al. (l97k). The values represented by the circles and 
by the upper limit designated T have been evaluated on the basis of the atmos- 
pheric experiments of Murcray et al. (1969 ) and of the laboratory data 
published by Fontanella et al, (197^). Constant volume mixing ratios are 

roughly illustrated by the straight lines marked 10 , 10 and 10 
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Figiire C-7.- Altitude profiles of nitrous oxide obtained from stratospheric 
bal-loon flights (Wyoming and Saskatchewan) and tropospheric aircraft flights 
(Denver environs) compared to a theoretical photochemical model. 
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Figure C-8.- Altitude profiles of CF^Cl^ obtained from stratospheric balloon 

flights (Wyoming and Saskatchewan) and tropospheric aircraft flights (Denver 
environs) compared to a theoretical photochemical model. 
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Figure C-9»- Altitude profiles of CFCl^ obtained from stratospheric balloon 

flights (Wyoming and Saskatchewan) and tropospheric aircraft flights 
(Denver environs) compared to a theoretical photochemical model. 


APPENDIX D 


ASSESSMENT OF UPPER ATMOSPHERIC EFFECTS 
OF SPACE SHUTTLE OPERATIONS 
By I. G. Poppoff,^ R. C. Whitten,^ J. B. Pollack,^ 
and R. G. Prinn^ 

INTRODUCTION 


Recent aeronomical investigations suggest that trace contaminants 
introduced into the air by anthropogenic activities may affect the stratos- 
pheric ozone (O^) layer. Some of these substances can catalytically destroy 

0^; that is, each molecule of pollutant gas can systematically remove any 

0^ molecules. For example, global 0^ could be reduced by more than 10 percent 

"by nitrogen oxides deposited in the stratosphere in the exhaust from a fleet 
of supersonic transports (refs. D-1 to D-i<), Stratospheric 0^ is important as 

a natural molecular barrier against excessive solar ultraviolet (uv) radiation. 

Concern has been expressed that hydrogen chloride (HCl) present in the 
exhaust plume of Space Shuttle launch vehicles could also deplete the 0- layer. 

U ^ 

The mechanism for 0^ destruction is catalytic (ref. D-5) , just as in the case 

of the nitrogen oxide destruction cycle. Although HCl itself is quite inert 
with respect to off-oxygen allotropes (O and O^), it is efficiently decomposed 

to chlorine (Cl) and water (H^O) by reaction with hydroxyl (OH) . The Cl then 

reacts with the odd 0 so that the overall effect is to accelerate the reaction: 


0 + 0 ^ 20 ^ 


(D-1) 


As will be shown later, the Cl cycle is quite effective in converting odd 0 to 
even 0. 

The deposition of HCl by Space Shuttle launch vehicles commences at launch 
and continues vmtil burnout at U5 kilometers {lk2 000 feet). The rate of HCl 
emission in metric tons per kilometer as recently determined (ref. D-6) is 


^ASA Ames Research Center. 

^Massachusetts Institute of Technology. 

^S. C. Wofsy and M. B. McElroy, lonpublished data, 197^. 
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shown in figure i+-3. At peak operational activity, it is expected that 60 
launches per year would be carried out. Launch operations are expected to 
commence in 1978 and to build up to 60 per year by 1984. In some of the 
studies reported here, the launch schedule shown in figure 4-4 was employed. 

During reentry of the Space Shuttle vehicle into the Earth atmosphere, 
substantial quantities of nitric oxide (NO) are produced in the shock-heated 
wake (ref. D-T). Estimates (refs. D-4 and D-8) of the accumulation of NO pro- 
duced in this manner in the altitude range of 65 to 80 kilometers (211 000 to 
268 000 feet) show that it does lead to a large local buildup of NO, which may 
last for a day. However, it spreads rapidly because of atmospheric motions, 
and in 1 week, the local NO concentrations are reduced nearly to the ambient 
condition. Calculations (ref. D-4) show that even if meridional spreading is 
very slow, the excess NO in the altitude region of 70 to 85 kilometers 237 000 
to 274 000 feet) would be at most a few percent of the ambient concentration 
and would have no discernible effect on stratospheric NO. The NO produced 
during Space Shuttle vehicle reentry will not be considered further in this 
report. 

Tiny particles, made chiefly of aluminum oxide (Al^O^), are contained in 

the exhaust products produced by the currently comtemplated fuel mixture for 
the Space Shuttle. Should future Space Shuttle vehicles use such a mixture, 
they will introduce particles into the stratosphere as they fly through this 
portion of the Earth atmosphere. As many of the added particles are smaller 
than a micrometer in size, they would typically remain in the stratosphere 
for periods on the order of a year; hence, a significant steady-state popu- 
lation of aerosols could be built up in the stratosphere. 

The AI 2 O 2 particles directly injected could serve as nucleation centers 

for sulfuric acid (H^SO, ) vapor and thereby enhance the normal stratospheric 

sulfate layer. Current work suggests that the particles in the normal strat- 
ospheric layer do form on preexisting particles. However, no studies have 
yet been performed that adequately consider the ways in which Al^O^ particles, 

serving as condensation nuclei, might alter the layer. Consequently, this 
discussion of Al^O^ effects will be restricted to the climatic effects of 

AlgO^ particles that have not interacted with vapor. 

One of the most important effects that the added AlgO^ particles might 

have is to cause a change in the Earth's climate by interacting with both 
incident solar radiation and terrestrially generated thermal radiation. An 
attempt has been made to assess the importance of these effects. A series of 
radiative transfer calculations was performed to estimate the sensitivity of 
the mean global temperature to the introduction of these particles. The 
change in the atmospheric opacity and the resultant change in global temperature 


^Hamill, Toon, and Kiang: 
Particles. To be published. 
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was compared with the optical depth changes that accompany particles in- 
troduced into the stratosphere by volcanic activity and with the temperature 
changes characteristic of past epochs of major climatic change. 


CHEMISTRY 


The HCl injected as rocket engine exhaust in the stratosphere can be 
converted into free Cl atoms principally by reacting with OH radicals : 


OH + HCl ^ H^O + Cl (D-2) 

The Cl atoms then react catalytically with odd 0 (O and 0^) 

Cl + 0^ CIO + 0^ (D-3) 

CIO + 0 -»■ Cl + 0^ (D-lt) 

to effectively remove O^. Aeronoraically , reaction (D-4) is the rate-limiting 
process in the 0^ destruction chain. Eventually, odd-0 catalysis is terminated 
when Cl is converted into HCl, mainly by reacting with CHj^: 

Cl + CHj^ -»■ HCl + CH^ (D-5) 

A secondary terminating reaction, similar to reaction (D-5), involves hydro- 
peroxyl radicals. 


Cl + HO^ HCl + 0^ 


(D-6) 


Although the HCl can be recycled into free Cl by way of reaction (D-2), a 
fraction of it is constantly being removed from the stratosphere by downward 
transport into the troposphere. Another important reaction in the Cl photo- 
chemical scheme involves the chemical dissociation of CIO by NO, 


NO + CIO -»■ NO^ + Cl 


(D-7) 


which, followed by rapid nitrogen peroxide (NO^) photolysis, effectively 
short-circuits the Cl catalysis chain. 


D-3 
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For the sensitivity of computed Cl-induced 0^ reductions to variations in 

the rate coefficients and species concentrations for processes (D-2) to (D-T), 
it can be shown that for small perturbations, the fractional columnar 0^ 

decrease AO^ is approximately proportional to the quantity: 

k,k k [0] [OH] 

AO cc — — — — (D-8) 

(kj^[CHj^l + k^[H02l) (kg[N0] + k2[0]) 

where an overbar denotes an average species concentration near the altitude 
of the greatest absolute 0^ decrement and k^ to kg are rate coefficients. 

Some detailed calculations^ support this simple 0^ reduction scaling rule, 

which can be used to obtain quick estimates of the changes in computer 0^ 

reduction resulting from changes in rate constants and ambient species abun- 
dances . 

The methane (CH^) is supplied to the stratosphere by upward transport 
from the troposphere and is decomposed mainly by way of the reaction 

OH + CHj^ -> HgO + CH^ (D-9) 

Because of this simple chemistry, CHj^ serves as a convenient tracer for 

determination above an altitude of 30 kilometers (95 000 feet) of vertical 
eddy diffusivity for use in one-dimensional models. Thus, CH|^ plays an 
important role in Cl chemistry through reaction (D-5). 

The chemistry of the oxides of hydrogen (HO^) system is of great 

significance to the catalytic reduction of 0 by oxides of chlorine (CIO^) 

because of reactions (D-2) and (D-6). The OH concentrations between 30 and 40 
kilometers (95 000 and 12^ 000 feet) have been recently measured (ref. D-9) 
over Texas at a solar zenith angle of 80°. The measured value of approximately 
T 3 

2.2 X 10 moleciiles/cm (±30 percent) at an altitude of lO kilometers 
(126 000 feet) is substantially larger than previously predicted (e.g., 

ref. D-i+),^ The key odd-H (HO^ and OH) reactions in this altitude range are 

the following. 


R. C. Whitten, I. G. Poppoff, and R. P. Turco: Comments on an Assess- 

ment of the Potential Impact on the Stratosphere of Solid-Fueled Rocket Engines. 
To be published. 

S. C. Wofsy and M. B. McElroy, unpublished data, 197^. 
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O(^) + H^o -*-® 20H 
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OH + HO^ -► H^O + Op 


10 

0 + HOp -*■ OH + Op 


11 

0^ + OH HOp + Op 


(D-IO) 


(D-11) 


(D-12) 


(D-13) 
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0 + OH -*■ H + 0 


13 

HOp + hv ^ OH + 0 


(D-lU) 


(D-15) 


H0„ + hv 




H + 0. 


(D-l 6 ) 


is ein energy state, hv is photon energy, and k 0 to 

are rate coefficients. The hydrogen peroxide (HpOp) reactions 

are omitted because they do not result in a significant net production or loss 
of odd-H species. For the atmospheric conditions that prevail between 3^ and 
1+0 kilometers (llO 000 and 126 000 feet), one can easily show that in daylight, 
equation (D-17) is an excellent representation of calculated OH profiles when 

the reaction rate coefficients listed in table D-I and predicted O(^), 0 , 

0^, HpO, and CH|^ concentrations are inserted. 


where 






[OH] 


(k^Q[0] + kg[0(^)] [HpO] 

^9 (^ll^S^ ^12 


1 
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The discrepancy between computed OH abundances and those measured by 
Anderson (ref. D-9) may be due to the rate coefficients and 

neither of which is known to be more accurate than an order of magnitude. 

For example, if these rate constants were equal (which is kinetically 

reasonable) and if rate coefficient were zero, one would obtain OH 

concentrations very close to the lower uncertainty limit of Anderson’s 
(ref, D-9) measured values. Additional uncertainties in the concentrations 

of 0, 0^, 0(^D), and H^O and in the rates of the other photochemical processes 

(D-10) to (D-l6) may account for residual differences in the OH abundance. 
Table D-I contains the important reactions and the most recent, available 
rate coefficient data (refs. D-IO to D-2l). 


DESCRIPTION OF TRACE CONSTITUENT MODELS 


The large number of chemical species and their reactions have been 
approached by use of models. So far, these models range in complexity from 
one dimensional to three dimensional. 


One-Dimensional Models 

One-dimensional models of atmosphere trace constituents include vertical 
transport and photochemistry. Transport is parameterized by an eddy diffu- 
sivity which relates vertical mass flux to the vertical mixing ratio gradient. 
In reality, the motions which distribute the various constituents throughout 
the atmosphere cannot be represented by a single parameter. However, with 
the selection of a diffusivity profile which leads to agreement between 
calculated and observed vertical distributions of convenient tracers such as 
CHi and carbon-li+, globally averaged vertical transport rates can be roughly 

simulated. The diffusivity profile shown in figure D-1 was employed ; this 
profile, a modified form of that suggested in reference D-22, has a 
13-kilometer {k2 000 foot) tropopause level characteristic of middle latitudes. 

One-dimensional models consist of a set of coupled continuity equations: 


9n. 9<{i. 

— + — i- = p 
at 3z i 


l.n. 

X 1 


(D-l8) 
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where t is time, z is the vertical coordinate, is the number density 
of the i-th constituent, and and are the production and loss rates, 
respectively, for the i-th constituent. The flux in equation (D-18) can 
he written as 


e M 


n. 
d — 


9z 


(D-19) 


where is the eddy diffusivity and n^ is the air density. 

All chemistry is represented in terms of and Equation 

(D- 18 ) is typically very "stiff"; that is, the characteristic time constants 
vary over many orders of magnitude. To overcome the problem of stiffness 
while still maintaining solution speed and precision, Turco and Whitten 
(ref. D-23) use the fact that in many aeronomic systems, there are families 
of species in which the members are rapidly intraconverted (or shuffled) but 
are only slowly converted to those outside the family. The family is con- 
served on a time scale much longer than the lifetime of many of its members. 
Accordingly, the equation governing the behavior of the total family con- 
centration is more stable than the equations for its members and may be 
accurately solved over longer time steps. The concentrations of the family 
members themselves are first approximately determined by using a suitable 
stable numerical technique and are then adjusted to correspond to the more 
accurate total family concentration. As in all mathematical models, upper 
and lower boundary conditions must be specified. To model diurnal variations 
accurately without inordinate demands on computer time, it is necessary to 
average diurnally both photolytic rates and concentrations of species that 

imdergo large variations between night and day.^ These modifications have 
been included in the Whitten-Turco model. Wofsy used a similar model for 
Space Shuttle effects computations, euid his work is summarized in the section 
of this appendix entitled "The Harvard Model." 


The Ames Research Center Two-Dimensional Model Of 
Stratospheric Trace Constituents 

In constructing a two-dimensional model (i.e., a model in a meridional 
plane) of trace constituents, one must modify equations (D-l8) and (D-19) 
so that 
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4>., the eddy diffusivity 

is a tensor K , V is the gradient vector operator, and V is the atmospheric 
~e ~ ~ 

advective or bulk velocity. 

The transport computations are time-split into those for vertical and 
horizontal advection and diffusion. Solutions for the transport processes 
are obtained from a mass-conserving, forward-time, space-centered finite- 
difference formulation of the equations. To optimize speed and accuracy, an 
implicit formulation is used for the horizontal diffusion and an explicit 
formulation is used for all other transport processes. Spherical geometry 
is employed in equations (D-20) and (D-2l) with the horizontal coordinates 
extending along a meridian from latitudes 80° S to 80° N in 5° intervals with 
the vertical coordinates extending from the surface to 60 kilometers (l95 000 
feet) in 2.5-kilometer (7000 foot) intervals. End boundaries are taken 
at latitudes 80° N and 80° S because meridional fluxes are expected to be 
small at these latitudes. The end boiandary conditions are taken to be zero 
flux for all constituents across the vertical boundaries. The upper boundary 
conditions are set by putting the constituents in mixing equilibrivim at a 
60-kilometer (195 000 foot) altitude, whereas the lower boundary condition 
for all components except nitric acid (HNO^)* ^2^2 

chemical equilibrium. Because HNO^j HCl, and H^O^ are water soluble, the 

mixing fractions are set equal to zero at the lower boundary. The number 

9 -3 

density of NOp is fixed at 3 10 cm at the lower boundary, whereas that for 

^ 11 —1 
©2 is fixed at 6 x lo cm . 

The mean meridional circulation is obtained by the kinematic method from 
the averaged equation of mass continuity. With the assumption that the den- 
sity field is in a steady state, the approximate form of this equation in 
spherical coordinates is 




where the mass flux <}>^ is now a vector quantity 


\ 


R cos if 


^(pv cos <t> ) 


d<P 


d(pv) 


9z 


= 0 
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where the overbar denotes an average with respect to time and longitude, 
p is the bulk density, v and w are the meridional and vertical velocity 
components, R is the radius of the Earth, and (j) is the latitude. Equation 
(D-22) implies the existence of a "stream function" il) for the total, mass 
flux such that 


2ttRpv cos 4> = - 

dZ 

(D-23) 

27TRPW cos <J. = ^ 


If the distribution of p and v is known, the first of these can be inte- 
grated vertically and w can be obtained from the second. 

To altitudes of 20 kilometers (63 000 feet), the density was obtained from 
data presented in reference D-2i*. Above that altitude, the density was ob- 
tained from vertical integration of the hydrostatic equation by using mean 
(rocket) temperatures. 

The V components were based on measured winds (ref. D-2U) to altitudes 
of 20 kilometers (63 000 feet). These values were then extrapolated to the 
top of the model by imposition of a simple vertical profile that matched the 
20-kilometer (63 000 foot) value. In addition, a meridional variation was 
imposed that resulted in a three-cell circulation structure (latitudes 80® N 
to 80® S) in the stratosphere and lower mesosphere dixring the summer and 
winter seasons and a fo\ir-cell structure in spring and fall. The first 
segment of equation (D-23) was then integrated vertically with the condition 
that ip = 0 on all boundaries. This condition implies that no net mass flux 
exists across a vertical latitude wall; therefore, the specified values of 
V had to be slightly adjusted so that the boundary condition was satisfied 
at the top and bottom of the model. The resulting circulation patterns are 
reasonable in the troposphere and lower stratosphere but need further refine- 
ment in the upper portion of the model. Except in the upper equatorial tro- 
posphere, meridional windspeeds are generally less than 1 m/sec. Vertical 
windspeeds attain maximums in excess of k mm/sec in the equatorial troposphere 
but are usually less than 1 mm/sec except in the vicinity of cell boundaries. 


The eddy fluxes are modeled by diffusion coefficients based on the method 
cited in reference D-25. The data given in reference D-26 were used at 
altitudes as high as 20 kilometers (63 000 feet). These values were extended 
to the higher levels by assuming that coefficient K is proportional to the 


variance of the v component and that the slopes of the mixing surfaces are 
proportional to the slopes of the mean isentropic surfaces. Given the 
other coefficients (ref. D-2U) were obtained from 
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(D-2h) 
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K 
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where a is the mean slope of the mixinf^ surface and is its variance. 

—t2 

The value of a was assumed to he constant. The variances of the v com- 
ponent were taken from reference D-27. 

The chemical rate equations are solved with an implicit technique. 

Omission of the transport term from equation (D-20) yields the finite-difference 
form 


1 

At 


n*? 

1 



(D-25) 


which can be linearized by taking the first term in a Taylor series expansion 

of about Q^, 

1 1 


= Q'l + 

At ^i 


J 


k ^ 



(D-26) 


The members of the set of mass conservation equations are coupled and require 
solution by matrix methods. With the implicit method, large time steps can be 
taken even though the equations are very stiff. 


Massachusetts Institute of Technology Model 

The Massachusetts Institute of Technology (MIT) three-dimensional model 
presently comprises 26 horizontal levels equally spaced in log pressure co- 
ordinates between the ground and 72-kiloraeter (237 000 foot) altitude. It 
solves a set of dynamical equations similaj" to the conventional midlatitude 
beta plane quasi-geostrophic equations except that it is global and allows 
for the variation of the Coriolis parameter with latitude. It also includes 
solution of the full continuity equation for O^. In the horizontal direction, 

predicted quantities are expressed in terms of a series of spherical harmonics 
using planetary waves with maximum wave niamber six and six degrees of freedom 
in latitude (rhomboidal truncation). Northern and Southern Hemisphere orography 
and tropospheric heating rates are also expressed in series form. The chonical 
scheme presently requires specification of the distributions of certain odd-N 


D-10 





I 


and odd-H species, but this limitation will be ultimately removed by inclusion 
of a larger set of chemical reactions. (This set is being implemented by 
recoding for the Illiac IV computer.) Allowance is made for radiative heating 
by 0^ and other gases. Chemical and heating computations are performed at a 

set of grid points chosen to filter aliasing errors during transformation from 
spectral to grid-point representations. Time integration of the model is 
performed using a Lorenz four-cycle scheme with four 1-hour steps in each 
cycle. This spectral model is extremely efficient computationally compared 
to existing stratospheric primitive equation grid-point models; it requires 
about UO seconds of IBM 360/95 time per stratosphere day. A full description 
of the first version of the three-dimensional model appears in reference 
D-28, Further runs, including a simulation of the stratosphere perturbed by 
supersonic aircraft effluents, are described in reference D-29. 

The MIT two-dimensional model extends from pole to pole and over the 8- 
to 68-kilometer (21 000 to 221 000 foot) region. Atmospheric motions are 
simulated by using mean vertical and meridional winds (presently obtained from 
the three-dimensional model), and eddy diffusion coefficients are deduced from 
observations (currently those computed by Luther, ref. D-26). A set of 36 
chemical reactions involving odd-N, odd-H, and Cl species is included (table 
D-l). The full continuity equation for long-lived species is solved on grid 
points spaced 5° in latitude and 2 kilometers (6300 feet) in altitude; a photo- 
chemical steady state is assvimed for very-short-lived species. The lower 
boundary condition for water-soluble species such as HNO^ and HCl is formulated 

in terms of a tropospheric rainout time. A full description of the first 
version of the two-dimensional model is given in reference D-30. 


Radiative Model 

The addition of aerosols to the Earth atmosphere can produce either a 
net warming or cooling depending on their optical constants and their size 
distribution function. Measurements of the transmission and reflection 
characteristics of Al^O^, as given in the scientific literatixre, were used to 

derive their optical constants — the real and imaginary indexes of refrac- 

7 

tion — from the near uv to the middle infrared. Estimates of the particle 
size distribution function were obtained from a combination of measurements 
of their size distribution in the exhaust wake and estimates of the change in 
the size distribution caused by sedimentation process (ref. D-31). 

Estimates of the mean temperature of the Earth can be obtained by equating 
the amount of solar energy absorbed by the entire Earth with the amount of 
thermal radiation it emits to space (ref. D-3l). With such calculations for 
perturbed and unperturbed conditions, it is possible to determine the change 


7 

0. B. Toon, J. B. Pollack, and B. N. Khare: The Optical Constaints of 

Several Atmospheric Aerosol Species: Ammonium Sulphate, Aluminum Oxide, and 

Sodium Chloride. Journal of Geophysical Research, in press. 
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in the mean surface temperature produced as a result of a given perturbation. 
To calculate the solar energy absorbed by the Earth, a computer program that 
permits an accurate solution of the multiple scattering problem was used. 
Allowance was made for scattering and absorption by both gases and aerosols, 
separate calculations were made at a number of wavelengths that spanned the 
near uv to the near infrared, and account was taken of partial cloud cover and 
the vertical inhomogeneity of the atmosphere. Mie scattering theory was used 
to derive the single scattering properties of the natural, background aerosols 
and the perturbing aerosols. 

A separate computer program was used to evaluate the thermal radiation 
at the top of the atmosphere, and both gaseous and aerosol opacity were con- 
sidered. Iteratively, this program also determined the surface temperature 
at which global heat balance was achieved. 


Ambient Atmosphere 

The Cl emitted in the stratosphere by the Space Shuttle will react with 
several ambient trace gases including 0, 0^, OH, NO, and CHj^. Atomic 0 has 

been observed only once in the stratosphere (ref. D-32). Figure D-2 demon- 
strates that the measured data agree quite well with the NASA Ames Research 
Center (ARC) model predictions at altitudes above 35 ■ kilometers (l08 000 feet). 
At lower altitudes, 0^ layering and backscattering of sunlight may affect the 

atomic 0 concentration, which is proportional to the sunlight intensity and 
0^ density. The only measured profile of stratospheric OH is that of Anderson 

(ref. D-9). His data for a solar zenith angle of 80® are shown in figure D-3, 

together with two profiles predicted by a one-dimensional model. ^ As discussed 
previously, the reasons for the discrepancy between prediction and observation 
are suspected but not yet definitely known. 

Numerous measurements of NO have been made in the lower stratosphere, but 
only a few have been made near 30 kilometers (95 000 feet). The one-dimensional 
model predictions agree within the experimental spread with the observational 
data shown in figure D-4 (refs. D-24 and D-33 to D-38). 

g 

Figure D-5 (refs. D-39 to D-42) shows the observed vertical CHj^ distri- 
bution and the one- dimensional model results obtained by using the eddy ^ 
diffusivity profile shown in figure D-1. Figure D-6 (refs. D-43 to D-46)° 
illustrates the measured ambient abundances of HCl and the model calculation 
In both cases, the agreement is satisfactory. 


R, C. Whitten, I. G. Poppoff, and R. P. Turco: Comments on sin Assessment 

of the Potential Impact on the Stratosphere of Solid- Fueled Rocket Engines. 

To be published. 

g 

°E. A. Martell, unpublished NCAR report, 1975. 
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Figures D-7 (refs. D-9 and D-U 7 ), D -8 (ref. D-32), and D-9 (refs. D-3^+, 
D-35 , and D-^ 8 ) show, respectively, some vertical profiles of OH, 0, and NO 
at various latitudes as computed with the ARC two-dimensional model. As with 
one-dimensional model results, the OH concentration profiles are substantially 
lower than Anderson's (ref. D-9) measured values. The computed NO profiles 
are very close to the measured profiles. However, the atomic 0 predicted con- 
centrations are smaller than those measured by Anderson (ref. D-32) because, 
unlike the ARC one-dimensional model, the two-dimensional model does not yet 
include diurnal averaging on concentrations. 

The MIT three-dimensional model predictions of zonal winds, stratospheric 
temperature, meridional circulation, and 0 ^ density compared to altitude and 

latitude, or season and latitude, are shown in figures D-10 (ref. D-i+9), 

D-11 (ref. D-I+ 9 ) , D-12, and D-13 (ref. D-50), respectively. Agreement with 
observations is good in each case. 


RESULTS OF SPACE SHUTTLE IMPACT ASSESSMENT 


The models have also been calculated using the estimated Space Shuttle 
exhaust products. Predictions from those calculations are described in the 
following section. 


Results Of One-Dimensional Model 

The potential effect of Space Shuttle exhaust on the 0^ layer was first 
estimated in reference D-51; an average hemispherical ©2 reduction of approxi- 
mately 0.3 percent for 60 launches per year was calculated. The estimate has 

recently been lowered to approximately 0.1 percent,^ based on the newly 
measured chemical reaction coefficient data presented in table D-I. The time- 
dependent 0^ reduction corresponding to the launch schediile shown in figure 

1+-1+ is shown in figure D-lU. 

The new estimates^ of the 0^ reductions beyond 1978 due to planned Space 

Shuttle la\uiches are shown in figure D-lU. Curve A represents the predicted 
0^ reduction using the new rate coefficients for the reaction of Cl with O^j 

Cl with CHj^, and CIO with NO (table D-l); the other rate coefficients are 

adopted from previous work. The eventual 0^ decrease, approximately 0.08 

percent, is about one-foiirth the estimated reduction (0.3 percent) reported 
previously. 


R. C. Whitten, I. G. Poppoff, eind R. P. Turco: Comments on an Assess- 

ment of the Potential Impact on the Stratosphere of Solid-Fueled Rocket 
Engines. To be published. 
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It should be noted that the effect of diminishing the rate constant for 
reaction (D-3) by a factor of approximately 0.5 and increasing the rate con- 
stants for reactions (D-5) and (D-T) by factors of approximately 1.8 and 1.7, 
respectively, is to decrease the 0^ reduction by a factor of approximately 6. 

This reduction factor is raised to approximately 9 by an increase of approxi- 
mately 50 percent in the predicted NO concentration resulting from the use of 
a new eddy diffusivity profile and to a diurnally averaged model. However, 
these -downward adjustments are partly offset by increases in the computed 0 and 
OH abundances by approximately 50 percent each. 

Curve B shows the 0^ reduction, approximately 0.11 percent at steady state, 

when the rate coefficient of the CIO + NO reaction (D-7) is reduced to its 

-11 3 -1 

older value, 1.7 10 cm sec , with all other rate coefficients remaining 

unchanged. Comparison of curve B with curve A illustrates the large change 
in the predicted 0^ reduction caused by the new data on the important reaction 

(D-7) between NO and CIO. 

Curve C shows the 0^ reduction when the rate coefficient for the OH + HO^ 
reaction (D-ll), which affects the OH concentration, is increased to the pre- 
viously used value of 2 x lO cm^ sec However, in view of Anderson's 
recent OH measurements, this high value seems unrealistic. 

Curve D shows the predicted 0_ reduction of approximately 0.11 percent when 
—11 3 —1 ^ 

^10 ” ^11 “ ^ ^ which leads to a better match with 

the OH distribution observed by Anderson; hence, curve D is probably a more 
realistic estimate of the potential 0^ reduction. 

Curve E illustrates the 0^ reduction when the rate coefficient for 

reaction (D-6) is increased to 10 cm^ sec~^; this value is almost certainly 
an upper limit for this reaction. Adjusting the rate constant (reaction 

(D-5)) to its older value of 5.1 10~^^ ^-1790/1 ^^3 would lead to an 

increase of approximately 70 percent in the predicted 0^ reduction shown in 
figure D-lit. 

The set of 0^ reductions (curves A to E) presented in figure D-l^t 

demonstrates the potential variability in predicted Space Shuttle effects. 
Although the cases treated are quite specific, they are all interesting be- 
cause they illustrate the effects of important aeronomic uncertainties on 
calculated 0^ reductions. Consideration of the uncertainties in reaction 

rates and ambient species concentrations leads to the conclusion that current 
estimates have an uncertainty of approximately 3. 

Figure D-15 shows the steady-state vertical distribution of the incre- 
ments in the concentrations of Cl, CIO, and HCl (labeled [Cl] , [CIO] , and 
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[HCl], respectively). Because these incremental concentrations are approximate- 
ly 0.1 of the computed ambient concentrations of the allotropes of chlorine 
(Cl ) constituents (e.g., the HCl distribution in figure D-6) , the steady- 

state 0^ reduction caused by Space Shuttle operations is approximately 0.1, 

the reduction caused by the ambient Cl^. 

To estimate the stratospheric 0_ recovery time, a case was considered 

■J g 

in which all launches are ended after a steady-state condition is reached. 

The corresponding 0^ reduction shown on figure D-lU (curve F) indicates that 

the recovery is very rapid; the ozone reduction decreases by a factor of 2 
in approximately 1.8 years. The 0^ recovery is aided by the fast diffusive 

redistribution of the injected Cl after the end of launch operations; excess Cl 
above 30 kilometers (95 000 feet) is quickly transported to lower altitudes 
where HCl is not strongly dissociated. It is subsequently transferred across 
the tropopause at a slower rate. The characteristic 0^ recovery time will be 

sensitive to the eddy diffusivity profile and will become longer as the 
diffusivity near 30 kilometers (95 000 feet) is decreased. 


Results of ARC Two-Dimensional Model 

The two-dimensional model described in section was adapted to assess 
the 0^ degradation by Space Shuttle launch vehicle emissions. The laxinch 

schedule shown in figure U-i; was employed, and all launches were assumed to be 
performed at the NASA John F. Kennedy Space Center (KSC) (latitude 30° N). 

The latitude dependence of 0^ column reduction after 10 years of operations is 

shown in figure D-l6 for the four seasons. The most remarkable feature is 
the strong corridor effect at latitude 30° N in sinmner. This seasonal effect 
is due mainly to larger OH and abundance 0 but smaller NO abundance during 
siararaer than during winter in the altitude range of 30 to Uo kilometers 
(95 000 to 126 000 feet). For the same reason, 0^ loss is greater during s\ammer 

than during the winter throughout the Northern Hemisphere. 

The buildup of the globally averaged 0^ distribution is shown in figure 

D-IT. Notably, the globally averaged reduction after 10 years is very close 
to one-half that obtained with the ARC one-dimensional model for the hemis- 
pherically averaged case, as shown on cmrve A, figure D-lU. However, the time 
required for buildup to a steady state is longer than 10 years. Hence, the cal- 
culated steady-state percentage of 0^ destruction will be somewhat larger than 


R. C. Whitten, I. G. Poppoff, and R. P. Turco: Comments on an Assessment 

of the Potential Impact on the Stratosphere of Solid-Fueled Rocket Engines. 
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that computed with the one-dimensional model. After 10 years of simulated 
Space Shuttle operations, further launches were stopped and the 0^ layer was 

allowed to recover. It was found that in approximately 3 years, the globally 
averaged 0 ^ column reduction had decayed to one-half its value when operations 

ceased; the decay time was approximately 50 percent longer than that obtained 
with the one-dimensional model. The investigation also showed that the pre- 
dicted corridor effect centered at latitude 30° N had disappeared within 1 
year after operations ceased. 


Results of MIT Model 

The MIT two-dimensional model has been rvin for a continuous insertion of 
HCL in a 5° latitude belt centered at latitude 30° N. The source magnitude 
was designed to simulate the effect of 50 Space Shuttle launches per year from 
KSC. From this run, an input distribution of CIO for the three-dimensional 
model was obtained. The three-dimensional model has not yet been run with this 
distribution of CIO. 


Radiative Effects 

The results of the radiative calculations (ref. D-3l) are summarized as 
follows. Figure D-l 8 shows the dependence of the global albedo on the optical 
depth enhancement Ax at 0.55 micrometer as a result of the added aerosols. 

The addition of Space-Shuttle-generated AlgO^ aerosols to the stratosphere 

woxild cause an increase in the Earth's albedo or, equivalently, the Earth 
would absorb less solar energy. 

The resulting change in mean surface temperature AT as a function of 
Ax is shown in figure D-19. The four curves of this figure illustrate the 
effect of alternative prescriptions for the change in water-vapor abundance and 
change in stratospheric temperature profile eind the inclusion of the aerosol 
effect on the thermal radiation. Curve A was calculated for a constant re- 
lative hxmiidity and a constant temperature lapse rate and is probably the 
most realistic of the cases shown. Curves B and C show the effects of assuming 
instead a constant absolute humidity and a constant stratospheric temperature 
profile, respectively. Finally, curve D shows what happens when the added 
aerosol opacity in the thermal infrared is neglected. Comparison of curves 
A and D shows that such a neglect is a poor approximation. Quite clearly, 
the aerosol effect at thermal wavelengths largely cancels the effect at visible 
wavelengths. The former causes a surface warming by enhancing the greenhouse 
effect. 

To apply these results to determine the impact of Space Shuttle traffic 
during the next several decades. Ax must first be evaluated. To do this, 

Mie scattering theory was used to relate Ax values to the steady-state mass 
loading of the added particles. The latter parameter was found from the mass 
of AI 2 O 2 generated by a single Space Shuttle vehicle as it traverses the 

stratosphere, from the projected number of Space Shuttle flights, and from 
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the mean dwell time of particles in the stratosphere (ref. D-32). For a 
projected traffic level of one Space Shuttle flight per week, it was found 

that At = 3.3 10 ^ for the Northern Hemisphere; hence, according to curve 

1 -4 

A of figure D-14, AT = 1.5 x lo K. The corresponding number for the 
Southern Hemisphere is three-sevenths as large. 

To assess the significance of these values, the At values were compared 
with those produced as the result of volcanic eruptions and the AT values 
with those typical of epochs of major climatic change. Figure D-20 illustrates 
the At values that have occurred during the last century, primarily as a re- 
sult of the introduction of material into the stratosphere by volcanic ex- 
plosions (ref. D-52). The At values expected from the particles produced 
by Space Shuttle vehicles are several orders of magnitude smaller than those 
generated by volcanic explosions. In fact, even relatively minor volcanic 

_2 

explosions typically generate a At on the order of 10 

The difference in temperature between a major ice age and the current 
epoch is approximatley 5 K; and the difference between 19^0 (the time of a 
recent temperature maximum) and the period from 1U5O to 1915 (the time of 
the "Little Ice Age") is approximately 0.5 K (ref. D-52). During the 
Little Ice Age, glaciers tended to expand and famines were more frequent at 
some higher latitude locations than during the warmer, preceding epoch. Hence, 
a AT of 0.5 K would be a serious change, but a AT value of less than 
0.1 K wovild be less significant. Since the AT expected from the Space 
Shuttle AI2O2 aerosols is about two orders of magnitude smaller than this 

limit, it seems unlikely that such a variation will cause a significant 
climatic change. 

The basic reason for the slight climatic impact of the added Al^O^ 

particles is that there are very few particles large enough to interact 
efficiently with light. Unless the ratio X = 2Tr(r/A) (where r is particle 
radius and A is the wavelength of light) is 1 or greater, particles are not 
efficient in scattering light. For visible light, therefore, r must be 
approximately 0.1 micrometer for X to be approximately 1. The Space Shuttle 
exhaust increases the number of 0.1-micrometer-sized particles in the strato- 
sphere by approximately 1 percent (ref. D-53) . Hence, At is small and no 
significant climatic change results. As pointed out earlier, Al^O^ particles 

could serve as nucleation centers for This process could cause small 

Al^O^ particles with X < 1 to grow into added H^SO^ particles with X > 1. 

The number of ^ 2^2 with r > 0.01 micrometer will exceed or equal 

the nTunber of natural nuclei of this size and will greatly exceed the nxomber 
of natural particles with r > 0,1 micrometer; therefore, such growth could 
be significant (ref. D-53). The growth of these particles has not yet been 
studied in detail; but because any increases in the stratospheric sulfate 
layer could have climatic implications, the problem should be piorsued. Until 
the nucleation and growth is understood, the climatic effects cannot be easily 
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estimated. Further studies to define the problem are being performed at ARC 
by using a model that has been developed for studying nucleation and growth 
of sTilfate aerosol particles. 

Results of Space Shuttle effects obtained by S. C. Wofsy are given in 
the following section. 


THE HARVARD MODEL 


To estimate the effect on stratospheric 0^, a model was adopted that 

provided for 60 Space Shuttle flights per year, steady-state conditions, and 
global spreading of the exhaust products. If one assumes that debris will be 
restricted to the Northern Hemisphere, the adjusted effect would be double the 
calculated global effect. Space Shuttle emission profiles from reference D-6 
were used. 

The global average source strength of CHI from 60 Space Shuttle launches 
is considerably smaller than that of the background sources caused by CH^Cl, 

CClj^, CF^Cl^, and CFC1_. For example, at an altitude of 30 kilometers 

(95 000 feet), 60 launches give an average HCl source of 0.15 cm'^ sec” 

3 —1 

compared to approximately 6 cm sec from the other sources of Cl; that is, 
the Space Shuttle contribution is 2.5 percent. At lower altitudes, the Space 
Shuttle contribution is slightly more important; e.g. , approximately 5.5 percent 
at 20 kilometers (63 000 feet). An 0^ reduction of 2 to 1+ percent from CL 

is predicted because of the background sources; thus, Space Shuttle effects 
lie in the range of a few tenths of 1 percent, or less. 

Meiking a quantitative estimate is hampered by the difficulty of suppress- 
ing numerical fluctuations of the order 0.1 percent for 0^ using a standard 

one-dimensional diffusion model. Nonlinearities of the 0^ equation are respon- 
sible for this problem. Therefore, the estimate is made indirectly by using 
another calculation in which the backgroimd sources were increased by approxi- 
mately a factor of 5. The following table gives the results and a comparison 
with Space Shuttle. 


Model 

Total Ozones, 

-2 

cm 

t total 

mixing ratio 
of HCl + CIO 



percent 

At 30 km, V/V 

^C1 

Standard 
, model 

9.8x10^® 

— 

1 . 05 x 10 "^ 

— 

Perturbation 
of Cl by 
5x 

8.67 

12 

U.75 

3.7x10"^ 
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The from 60 Space Shuttle flights is 2.8 x 10 ^^(V/V) at 30 kilometers 

(95 000 feet); thus, a scaling of results from the table predicts 0.1 percent. 
If the effluent from Space Shuttle flights were restricted to the Northern 
Hemisphere, the Space Shuttle effect would be approximately 0.2 percent. 

The 0.2-percent estimate is on the upper end of the range of values 
predicted by the model for three reasons. 

1. There should be considerable interhemispheric exchange. 

2. Steady-state conditions have been assumed. It is not clear whether 
flights will continue long enough to approach a steady state. 

3. Nonlinear interactions exist between oxides of nitrogen (NO^) and 

CIO^ systems which tend to decrease the CIO^ perturbation over certain limited 

concentration ranges. Neglect of these factors will tend to enhance the cal- 
culated perturbations. 


SUMMARY 


From the foregoing discussion, one can conclude the following about 
stratospheric effects of Space Shuttle vehicles. 

1. The hemi spherically averaged 0^ reduction is of the order 0.1 percent 

subject to a significant uncertainty (approximately a factor of 3) associated 
mainly with uncertainties in pertinent chemical rate coefficients and in 
vertical distributions of important atmospheric constituents (i.e., OH, 0, 

NO, and CHj^). The ARC and Harvard University results agree within the un- 
certainty. 

2. Studies using two-dimensional models revealed a strong corridor effect 
(centered at altitude 30° N) in summer but not in winter. The summertime 0^ 

reduction percentage at latitude 30° N is about twice that at latitude 50° N. 

3. The reentry of Space Shuttle vehicles into the Earth's atmosphere will 
produce NO in the altitude range of 65 to 80 kilometers (211 000 to 268 000 
feet). However, it will be rapidly dispersed (about 1 week) over the Earth, 

so that the long-term addition to the ambient mesospheric NO will be at most 
a few percent. The effect of such NO on stratospheric 0^ will be negligible. 

U. The direct climatic effects of the AI 2 O 2 aerosol are likely to be 

negligible. That is, the decrease of mean surface temperature resulting from 
the scattering of sunlight by Al 20 ^ particles is estimated to be of the order 

of thousandths of a degree. However, it is possible that the small particles 
may act as condensation nuclei and thereby significantly increase the abundance 
of large H^SOj^ droplets. The condensation nuclei aspect of the problem is 

being assessed at ARC. 
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ADDENDUM 


Since preparation of this paper, the importance of chlorine nitrate 
(CINO^) has become apparent. Odd N is responsible for approximately 25 percent 

of the total loss rate of stratospheric O^. Were it not for the nitric acid 

reservoir, the destructive effect of NO^ on 0^ would be substantially greater; 

in the 20- to 30-kilometer (63 000 to 95 000 foot) region, CIWO^ acts as an 

additional NO reservoir and thus decreases the 0_ loss rates caused by the 

presence of NO . Hence, Cl both reduces and increases 0_,; the former effect 

occurs above 30 kilometers as stated in this paper and the latter effect 
below 30 kilometers. Because of the subtractive nature of the two effects, the 
total change in 0^ caused by the presence of stratospheric Cl is very sensi- 

- je to the rates of formation and destruction of CINO^. In a paper to be 

published elsewhere, one of the authors (Whitten) and his associates calcu- 
lated that when CINO^ is included in the computations, 0^ loss caused by Cl 

will be between one-half and one-fourth the values obtained when C1N0- is not 
included. ^ 
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TABLE D-1,- REACTION RATE COEFFICIENTS^ 


Number 

Reaction 

Rate coefficient^ 

Uncertainty 

References 

R1 

OH + HCl -*• Cl + HgO 


±20 percent 

D-10, 





D-11 

R2 

Cl + O3 CIO + O3 


±15 percent . 

(b) 

R3 

CIO + 0 -*■ Cl + Og 

5.3x10“^^ 

±20 percent 

D-12 

ru 

Cl + CHjj -»■ HCl + CH^ 

5.I..10-1" .xp-1125/T 

±20 percent 

(c) 

B5 

Cl + .HO^ -*• KCl + Og 

IxlO"^^ 

± factor of 10 

(d,e) 

R6 

CIO + NO Cl + NOg 

6.3.10-1= 

±50 percent 

(f) 

r8 

OH + CH|^ -*■ HgO + CH^ 

2.1.10-1= .xp-ll=”« 

— 

D-ll* 

R9 

O(^) + HgO 20H 

3.5x10"^° 

±30 percent 

D-1 5 

RIO 

OH + HOg -► HgO + Og 

6x10"^^ 

± factor of 3 

(e) 

Rll 

0 + HOg -► OH + Og 

1x10“^^ 

± factor of 3 

(g) 

R12 

0^ + OH HOg + Og 

1.3.10-1= 

±50 percent 

D-16 

R13 

0 + OH H + O2 

3.8x10“^^ 

±50 percent 

D-17 

R16 

H ♦ Og + M -► HOg + M 

1.9x10"^^ exp^^^^^ 

— 

D-18 

R17 

HOg + HO2 - H2O2 + O2 

1.7x10"^ exp"^°°'^’’ 

— 

D-15 

Jll* 

HO- + hv OH + 0-\ 






1.0x10“-^ 

(e) 

D-19, 

J15 

HOg + hv -*■ H + O2 J 



D-20 

Jl8 

HgOg + hv 2 OH 

l.Uxio"^ 

~ 

D-21 


units per second for unimoleculax processes, cubic centimeters per second 
for bimolecular reactions, and hexadic centimeters per second for termolecular re- 
actions. The symbol T represents ambient temperature. 

^This value is based on a least squares best fit to recent measurements by 
Davis and Watson; Zahniser, Kaufman, and Anderson; end Nip and Cline (19T^). 

Q 

Data frcan R. T. Watson, to be published, 

^eu and De More (ref. D-13) have very recently reported the rate constant for 

nuDober R5 to be 3 ("*^2*8^ ^ 10~^^ cm^ sec”^. Use of the mean value rather than that 

listed vill lead to a decrease in predicted values of 0^ reduction of approximately 
20 percent. 

^Estimate (W. B. De More, private communication). 

f 

Data from M. S, Zahniser and F. Kauftaan, to be published. 

®F, Kaviftnan: Neutral Reactions. DASA Reaction Rate Handbook, M, H. Bortner, ed. 


R. C. Whitten, I. G. Poppoff, and R. P, Turco: Ccsmnents on an Assessment 

of the Potential Impact on the Stratosphere of Solid-Fueled Rocket Engines. To be 
published. 


D-26 


II fl 1 1 


111 


1 £ 1 1 1 


I i £ £ I L 







concentration, cm 


ions for a solar zenith angle of 80° (experi 
points from ref. D-9). 



M 
J* 


Legend 

— 

From 

ref. 

D-33 

o 

From 

ref. 

D-34 

• 

From 

ref. 

D-35 

□ 

From 

ref. 

D-36 

— 

From 

ref. 

D-37 

8#P 

From 

ref. 

D-38 


Model of McConnell 
and McElroy (from ' 

ref. D-39) 


ARC one-dimensional model 


i May 




Number density, cm 

Figure D-U,- Nitric oxide concentrations, 
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Figure D-7.- Hydroxyl concentration profiles obtained using the ARC 
two-dimensional model (day 7823). 
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Figure D-9.- Nitric oxide concentration profiles obtained using the ARC 

two-dimensional model. 
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(b) Modeled (MIT three-dimensional model for season 10), 


Figure D-10,- Vertical cross section of mean zonal windspeeds. Contour 

interval is 20 m/sec. 






(a) Observed (ref, D-U 9 ). 



(b) Modeled (MIT three-dimensional model for season 10). 

Figure D-11,- Vertical cross section of mean zonal temperature distributions. 

Contour interval is 20 K, 
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(a) Determined by Dtitsch (ref. D-50) from observations. 
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(b) Predicted during year 3 in the unperturbed model run. 

Figure D-13.- MIT three-dimensional model columnar 0^ density in the Northern 
Hemisphere as a function of season and latitude. Contour interval is 20 cm' 
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Figure D-lU,- Computed (ARC one-dimensional model) 0^ reductions for the HCL 

deposition rate shown in figure lt-3. Solar photolysis rates recomputed 

whenever the vertical ozone column above 30 or 10 kilometers changed by more 

than 0.01 percent. Curve A represents the reduction if Space Shuttle 

launches are continued indefinitely, and the rate coefficients listed in 

table D-I are used. For curve B, the rate coefficient of reaction (D-T) was 

, -12 351+/T -11 3 -1 

changed from 6,3 ^ 10 exp to 1.7 ^ 10 cm sec . For curve C, 

—11 3 — 1 

the rate coefficient of reaction (D-ll) was changed from 6 x 10 cm sec" 

to 2 X 10 ^ cm^ sec"^. For cvtrve D, the curve A parameters were used, but 

—11 3 —1 

= k^_ = 6 X 10 cm sec was set. For curve E, the rate coefficient 

—10 3 —1 

for reaction (D-6) was increased to 10 cm sec . For curve F, the 
ciirve A parameters were used, but launch operations were considered to be 
terminated after 10 years. 
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Figure D-15,- Incremental steady-state vertical distributions of Cl, CIO, eind 
HCl resulting from Space Shuttle operations. 


D-i*l 


llIMIIIIfillllLIIIIL 






i 



South North 
Latitude, deg 


Figure D-l6,- Computed 0^ column reduction as a function of latitude after 10 
years of Space Shuttle operations, using the ARC two-dimensional model. 
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Figure D-IT.- Time dependence of the globally averaged 0_ reduction using the 

ARC two-dimensional model. 
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O Curve A = constant relative humidity 
□ Curve B = constant absolute humidity 
O Curve C = constant relative humidity and 
stratosphere temperature 
A Curve D = constant relative humidity and 
infrared aerosol opacity 


Figure D- 19 .- Change in the mean surface temperature AT as a function of the 
optical depth perturbation At of AI2O2 aerosols. 
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Figiire D-20,- Pertcirbations of optical depth caused by volcanic explosions 

during the last century. 



APPENDIX E 


PARTICULATE MEASUREMENTS 
By G. Varsi^ 


INTRODUCTION 


After the publication of the July 1972 environmental statement for the 
Space Shuttle Program, questions arose on some of the characteristics of 

9 

the particles emitted by the solid rocket boosters. An experimental program 
was stsirted to resolve some uncertainties and to provide better characteri- 
zation of the particles. Because existing data were obtained mostly from 
firings of motors considerably smaller than the Space Shuttle booster when 
large particles were of principal concern, because two-phase aerodynamic 
losses occurred, and because scaling laws are not established, it was 
decided that new measurements were necessary on the largest solid rocket 
motor possible and that an effort should be made to detect the small particles. 
Current launches of Titan III C/D vehicles appeared the best of available 
choices. A comparison of Space Shuttle and Titan boosters is given in table 
E-I. 


SCOPE OF WORK 


In the Space Shuttle Program, interest was focused on particles emitted 
by the boosters; no consideration was given to particles formed by atmospheric 
reactions of gaseous or liquid compounds of such elements as nitrogen (N), 
hydrogen (H), and chlorine (Cl) (e.g. , NH^ + HCl NHj^Cl). Futhermore, 

because of program constraints (funds and time), no instinment development 
was thought possible. 

PROGRAM PARTICIPANTS 


The following organizations and individuals participated directly in the 
flight measurements: NASA Ames Research Center (ARC) — N. H. Far low, G. V. 

Ferry, and W. A. Page; NASA Jet Propulsion Laboratory (JPL) — R. E. Gaiildin, 
E. G. Laue, L. Strand, and G. Varsi; and Los Alamos Scientific Laboratory 
(LASL) — W. A. Sedlacek. The following organizations and individuals 


^ASA Jet Propulsion Laboratory. 

9 

Letter from Sheldon Meyer, Director of the Office of Federal Activities, 
Environmental Protection Agency, to Ralph E. Cushman, Special Assistant, Office 
of the Administrator, NASA, June 26, 1972. 
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participated in the analysis of the data: the University of Southern 

California — E. Hall; EMS Laboratories; and the University of California 
(Lawrence Berkeley Laboratory) — W. Siekhaus. 


AIRCRAFT AND INSTRUMENTS 


The following additional instrumentation was supplied: (l) the tape 

impactor supplied to JPL by EC&G, Santa Barbara (L. Franks and D. Seaver), 

(2) the photographic equipment supplied to ARC by Lockheed Palo Alto Laborato- 
ries ^ and (3) the photographic equipment supplied to JPL by Edwards Air Force 
Base (Hiroshi Hoshizaki). The experimental platforms used and some of the 
principal characteristics of the instrimientation are shown in table E-II. 


HISTORY AND STATUS OF MEASUREMENTS 


The 
(WTR) in 


following Titan launches have been monitored at the Western Test Range 
California and at the Eastern Test Range (ETR) in Florida: 


Date 

Location 

Measiirement group 

Apr. 197^ 

WTR 

ARC, JPL, LASL 

Oct. 19lh 

WTR 

ARC, JPL 

May 1975 

ETR 

JPL 

Jvine 1975 

WTR 

JPL 


The participants believe that usef\il preliminary data have been obtained which 
significantly correct previous existing knowledge,, The abrupt interruption of 
the activities in August 19T5 prevented the performing of final measiirements . 
All flight experiments have been suspended. Reporting of past activities is 
in progress. 


RESULTS 


Most of the available data pertain to the size of the particles; limited 
information exists on the composition. 


Particle Size 

Data obtained in April 197^ from the coated-wire impactor and from the tape 
impactor by ARC and JPL, respectively, are shown in figure E-1. The agreement 
is quite satisfactory, especially when one considers the difference in tech- 
nique and the fact that the ARC data represent the stim of nearly six transits 

E-2 
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through the pl'jme in the hour following launch of a Titan III-D from WTR, 
whereas the JPL data were obtained during a single transit approximately 10 
minutes after launch. 

From the observation of the filters exposed during the same Titan Launch, 
it is quite apparent that the majority of the particles are smaller than a 
micrometer. The mean diameters computed from both distributions are similar 
and have the following values. 

Dj = 0.1 micrometer (number average) (E-l) 


Q,P micrtimoti-r (volume- average) (E-2) 


Furthermore, the distribution can be well approximated by a power law in the 
range of 0.07 to approximately 10 micrometers. Using an intermediate value 
from the two distribut ions , one obtains 


_k 5 

dU = N D dD; 

o ’ 


0.07 micrometer ^ D < 10 micrometers 


(E-3) 


where N is the number of particles. More recent data obtained by JPL with 
an electrical mobility analyzer are reported in figiores E-2 and E-3. Again, 
satisfactory agreement is observed, but the extension of the distribution to 
smaller sizes by means of this more sensitive instrument seems to show a peak 
in the distribution around 0.07 to 0.08 micrometer, just at the lower limit 
of the impactor data. 

In conclusion, although the experimenters consider these data good pre- 
liminary measurements (rather than definitive determinations), the substantial 
agreement between them suggests that equations (E-l), (E-2), and (E-3) and the 
information presented in figures E-l, E-2, and E-3 can be reasonably used to 
represent the particles found in the plume shortly after launch. 


Particle Concentration 

Less satisfactory information is available on absolute particle concen- 
trations; it is difficult to estimate the effective volume swept by the wire 
impactor because the boiindary of the plume is poorly defined. The tape impac- 
tor is immune to this problem, but it requires absolute calibration of its 
collection efficiency. Because of the program interruption, such a calibration 
has not yet been obtained. Filters have the same difficulty as the wire 
impactors. Under certain flight conditions, the electrical mobility analyzer 
can measure absolute concentrations (when the sampling time is shorter than the 
transit time); however, it was not clear whether such conditions existed during 
the preliminary flight measurements performed before the interruption of the 
program. 

E-3 
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The closvire condition, based on calculated booster emission at an altitude 
of approximately 20 kilometers (63 000 feet) of 930 grams of alumimim oxide 

(AlgO^) per meter altitude, a density for Al^O^ of approximately U g/cm , a 

plume diameter of 2.0 kilometers (6563 feet), and a volume-averaged particle 
diameter of 0.2 micrometer, yields an average concentration C. 


^ 930 g m 

TT 10^ m^m^ TT (10 ^)^ cm^ x g cm ^ 

= 1.76 X 10^*^ particles/m^ (E-4) 


As shown, estimates obtained for the wire impactors and the electrical mobility 
analyzer are in the 10^ m~^ and 10^^ m~^ range, respectively. 

10 3 

The Aitken nuclei counter measured at least 3 x 10 particles/m . Over- 
all, there is agreement within a factor of 2. Of course, a 50-percent increase 
in the diameter of the particles or in the size of the plume can each accommo- 
date about half an order of magnitude. The particle diameter has been discussed 
in the section entitled "Particle Size," whereas the estimate of 2 kilometers 
(6300 feet) for the plume diameter was obtained from a common pilot report 
that approximately 10 seconds are required to cross the plimie at a speed of 
206 m/sec (1+00 knots) approximately 15 minutes after launch and at a 
20-kilometer (63 000 foot) altitude. It appears that a reasonable assumption 

10 3 

for calculations is approximately 10 particles/m in the early plume. 


Particle Density 

Very preliminary measurements were performed at JPL using a gravimetric 
technique to determine the density of the particles collected. Several authors 
reported that many of the spheres collected from rocket firings have 

hollow interiors and therefore an apparent or aerodynamical density inferior 

3 

to the theoretical density of approximately 3.97 g/cm . The preliminary 

3 

measurements have yielded a wide range of densities from 1.5 to 3.5 g/cm ; no 
indication of relative abundances was given. 


Particle Appearance 

Spherical objects were mostly reported, with occasional agglomeration 
producing a covering of, for example, a 10-micrometer sphere with perhaps 
fifty 0.1-micrometer particles. Various textures have been observed. Some 
particles are smooth to a resolution of 20 nanometers; others are finely grained 
and deeply grooved, an indication of a polycrystalline structure (fig. E-U). 


E-i+ 


]| II 1 1 I 1 I 


fi I H I £ 


£ 1 i I i 


\ 



Particle Structiire and Composition 

The crystaU-ine structure of ARC-collected AlgO^ particles was exeunined 

by EIC Laboratories, Pasadena, California. The structures were foxind to be 
mostly hexagonal (a) phase for particles in the 0 . 1 -micrometer size range. 

Both transmission and reflection electron diffraction were used. Significant 
amounts of iron oxide were also fovind; these concentrations presumably were 
caused by a propellant additive for controlling the burn rate. However, 
earlier transmission diffraction measurements on a few 0 . 1 -micrometer particles 
collected by JPL using an impactor showed cubic (y) structure for the Al^O^. 

Other measurements on material collected by JPL from small motors fired in a 
closed tank also indicated a cubic structxire. At present, it would appear 
that the evidence is in favor of hexagonal (a) structure for the material in- 
jected in the stratosphere. 

Elemental comjosition of the particles was obtained by ARC, JPL, and 
LASL, mostly by X-ray fluorescence. Besides aluminum and occasionally potas- 
sixim, sodium, end titanium, significant quantities of iron and silicon were 
found. Very rarely has chlorine (Cl) been observed by X-ray fluorescence. 
However, a series of samples collected by JPL on Nuclepore filters and analyzed 
by electron spectroscopy (ESCA) at the Lawrence Berkeley Laboratory revealed a 
significant, often leirge signal from Cl in varying degrees of oxidation 
ranging from a tight bond (as in a chloride) to a free atom (as for surface 
absorption). It was reported in the literature that such a variability in 
the bond energy of Cl deposited on metal surfaces is normal and dependent on 
the Cl surface concentration. From the ESCA shift of the aluminum peak, it is 
presumed that various compounds of altminiam, chlorine, and oxygen may be formed. 

The apparently contradictory evidence for Cl obtsdned by electron spec- 
troscopy and X-ray fluorescence can probably be best explained by the different 
penetration capabilities of the two techniques. The typical signal for X-ray 
fluorescence originates from a layer a fraction of a micrometer thick, whereas 
the ESCA signal is produced in only a few nanometers. If it is hypothesized 
that Cl is not uniformly distributed in the molten AI 2 O 3 is deposited on 

the surface after the particle has frozen to form an almost complete monolayer, 
one has a situation in xrtiich ESCA may detect 10 to 20 percent Cl in the first 
few layers but the X-ray yield of this small quantity of Cl would be undetect- 
able. Whether the less tightly boxuid portions of such a layer can be released 
to the surroxinding stratosphere under some conditions is not known. 


CONCLUSIONS 


The xinanimous agreement of the various experimenters is that most of the 

particles are smaller than a micrometer and are typically 10 ~^ centimeters 
in diameter, whether a number, surface, or volume average is used. This new 
finding contrasts with widely quoted previous estimates that the mean diameter 

_3 

is of the order of 10 centimeters. The immediate and obvious consequence 
is that the typical residence times in the atmosphere are much longer (about 
two orders of magnitude) . 
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Number concentrations in the early plume of the order of 10 particles/m 
axe probably correct. Other properties of the psirticles are less certain, and 
the results should be viewed as preliminary. 
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TABLE E-I.- COMPARISON OF TITAN AND SPACE SHUTTLE BOOSTERS 


Peirameter 

Titan III 
C/D 

Space 

Shuttle 

Total mass of propellant (2 boosters), 


O 

kg 

38Ux10^ 

922x10 

Nozzle throat diameter, cm 

95.76 

136. ll* 

Propellant formvilation 



Total solids, percent 

8U 

81+ 

Aluminum, percent 

l6 

l6 

Binder 

PBAN^ 

PBAN^ 

Combustion chamber pressure, MPa 



(bar) 

3 to U.5 (30 to 1*5) 

3.5 to 5.5 (35 to 55) 


^olybutadieneacrylonitrile . 


E-T 


]l I 1 1 I 


I I 

\ 




1 


I I 2L & 


i A il I £ 1 I 




TABLE E-II.- EXPERIMENTAL PLATFORMS AND INSTRUMENTATION 


Organization 

Aircraft 

Instruments 

ARC 

NASA U-2 

Wire impactors (coated and uncoated); 
collection efficiency estimated unity 
for particles greater than 0.07 um 


Lear jet 

Photographic cameras 

JPL 

USAF®' U-2 

Continuous tape impactor (supplied by 
EG&G; particle size collection range 
estimated at 0.1 ym to more than 50 ym 

Nonfibrous membrane multiple filter; 
useful collection range estimated at 
0.005 to 10 ym with varying efficiency 

Electrical mobility analyzer; useful 
particle size range from 0.02 to 1 ym 
with automatic size classification 

Panoramic camera (supplied by USAF) 

LASL 

WB-57F 

Stratospheric Aitken nuclei detector 
systems; sensitivity from 0.005 to 
0.1 ym; no classification capability 

Fibrous filters; sensitivity of 0.1 ym 
and larger 


^USAF = U.S. Air Force. 
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Figure E-1.- Impactor data for April 197^, at an altitude of 20 kilometers 
(65 000 feet), for particle diameter D, 






kBH 

tad 

i 

I 

Wi. 

m 

pt 

p9 

pa 

pa 

irw 


t?d 

I 

H 

O 



ETR, May 20, 1975* 7 hours after lift-off, 
IN = 7,1+1 X lO^/m^. 





(c) D, /im (d) D, fim 

WTR, June 8, 1975* 6.5 hours after lift-off, WTR, June 8, 1975, 13.5 hours after lift-off, 

EN = 1.1+7 X 10^°/m^. EN = 2.33 x lO^/m^, 

Figure E-2.- Electrical mobility analyzer data at an altitude of 19 kilometers, 

where N is the number of particles. 



.01 .1 1.0 
D, /i.m 

10 

(a) Niimber, EN = l.i+T ^ 10 /m . 



Figvire E-3.- Electrical mobility analyzer data at an altitude of 19 kilometers 
for WTR, Jime 8, 1975* 6.5 hours after lift-off, where S is surface and 
V is volume. 
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Figiire E-U.- Scanning electron microscope photograph of Titan III-D exhaust 

particle. 
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APPKNDIX F 


ESTIMATES OF TITE OLTMATIO IMT’AO'r’ 

OF AEROSOI.S PRODIJCEP BY SPACE SHUTTLES, 

SST's, AND OT’HER HIOH-FLYINC AIRCRAFT* 

By James B. Pollack,^ Owen B, Toon,^^ 

Q. ft ft 

Audrey Summers, Warren Van Camp, and Betty BaJdwin 


ABSTRACT 


Aircraft and Space Shuttle vehicles flyinp: in the next several decades 
will add hydrogen sulfate and aluminum oxide particles, respectively, to the 
stratosphere. To evaluate the effect of these additional aerosols on the 
global heat balance, solar and terrestrial radiative transfer calculations 
have been performed. The solar calculations employed an accurate numerical 
method for solving the multiple scattering problem for unnolarized light to 
determine the dependence of the global (spherical) albedo on the optical depth 
pertiirbation At, Correct allowance was made for absorption by gases. With 
these results and those obtained from calculations of the terrestrial, thermal 
flux at the top of the atmosphere, the resulting change in the mean surface 
temperature AT was determined as a function of At. In both calculations, 
the measured optical constants of the aerosol species of interest were used. 

To apply these results to the problem of interest, engine exhaust prop- 
erties for the various types of vehicles were used to estimate their optical 
depth perturbation, and the record of past climate changes was examined to 
set a threshold value, 0.1 K, on the mean surface temperature change, below 
which no significant impact is to be expected. Use of this information showed 
that no significant climate change should result from the aerosols produced by 
Space Shuttle, supersonic transport, and other high-altitude aircraft operating 
at traffic levels projected for the next several decades. However, the effect 
of supersonic transports is sufficiently close to the threshold limit to war- 
rant a reevaluation as their characteristics are updated. 


* 

This work was supported in part by the Climatic Impact Assessment 
Program of the U.S. Department of Transportation and by the National Science 
Foundation (grant NSF-DES T^-20l+5B). 

^ASA Ames Research Center. 

^Work initiated while at Cornell University. 
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INTRODHCTJON 


In recent years, a growinp concern that pollutants producer) by this tech- 
nological civilisation may inadvertently change the Earth's climate on a global 
scale has been evident. Much attention has been directed toward material 
introduced into the troposphere. It was pointed out that increases in the 
carbon dioxide (CO^) and aerosol content of the atmosphere could change the 

rve'rall heat balance of the Earth and thereby alter the mean surface tempera- 
ture (refs. F-1 and F-2). Recently, the debate over the desirability of 
fleets of supersonic transports has engendered a careful scrutiny of the in- 
tent ial effects of their exhaust products. Although considerable emphasis 
has been placed on evaluating the change in the ozone content of the 

atmosphere and its resultant biological consequences, much effort has also 
been expended in evaluating the climatic implications of the aerosols that 
will be added to the stratosphere. Similar questions have been raised on the 
effects of the exhaust products expected from Space Shuttle vehicles. 

Aerosols introduced by vehicles flying through the stratosphere might be 
climatologically significant for two reasons. First, since the residence time 
of particles smaller than 1 micrometer in the stratosphere is on the order of 
years, large steady-state perturbations in the number of particles might be 
created. By contrast, similarly sized particles in the troposphere have life- 
times only on the order of weeks. Second, since these vehicles can be ex- 
pected to operate for an indefinitely long period of time, the relatively 
large thermal inertia of the atmosphere surface system offers no effective 
barrier to changes in the mean surface temperature. This paper contains es- 
timates of the change in mean surface temperature caused by aerosols intro- 
duced into the stratosphere by supersonic transports, other aircraft flying 
at high altitudes, and Space Shuttle vehicles. Radiative transfer calcula- 
tions are performed for visible and infrared wavelengths to determine the 
effect of the aerosols on both the solar and thermal energy budgets. Because 
many processes affect the climate of the Earth, changes in the radiative energy 
balance do not completely determine the response of the climate to a perturb- 
ing influence. Although this study neglects atmospheric dynamics and many 
possible climatic feedback loops, it provides a useful estimate of the climatic 
significance of increased aerosol levels. 

Niimerous attempts have been made to calculate the radiative effects of 
aerosols on the Earth's climate in recent years, and the reliability and 
sophistication of the calcvilations have increased rapidly. In reference F-2, 
the effect of an increased particulate burden in the lower troposphere is 
investigated, and it is found that a cooling of the Earth would result. 
Plausible values for the optical constants and for the size distribution 
function of the aerosols were used in performing calculations at one represen- 
tative wavelength in the visible and one in the infrared. For the visible 
albedo calciilations, the two-stream approximation (ref. F-3) to the multiple- 
scattering problem was used. The change in surface temperature was fo\ind by 
balancing the amount of sunlight absorbed with the thermal emission at the top 
of the atmosphere. 
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A number of the subsequent papers dealing with this problem contained 
improvements in the numerical method used to determine the visible albedo. 

Much effort has also gone into investigating the effect of changing various 
parameters that are used in such calculations. In reference F-ii, the principle 
of invariance (ref. F-5) was used to obtain an acciu*ate solution of the radia- 
tive transfer problem at visible wavelengths. It was found that the sign of 
the albedo change caused by an enhanced turbidity depended sensitively on the 
value of the imaginary index of refraction of the aerosols. Large values 
(-^^•05) of the refractive index led to a lower albedo and therefore to warm- 
ing of the surface. 

Improved estimates of the albedo change by including the effects of 
gaseous absorption are cited in reference F-6. However, these calculations 
were inexact in this regard because wavelength-averaged transmission functions, 
which do not obey Beer's law, were used to describe the opacity properties of 
the gases. Because the mmierical technique provided solutions to the multiple- 
scattering problem only for small optical depths, tropospheric water clouds 
could not be included in the model. 

In reference F-7, a modified two-flux approximation was used to study the 
effect of aerosols in both the visible and infrared. It was found that the 
effect of the aerosols on the thermal radiative fluxes cannot be ignored and 
that the sign of the aerosol effect on the surface temperature depends on the 
value of the surface albedo, with net heating occurring when the surface 
reflectivity exceeded 0.3. In the context of the supersonic transport (SST) 
problem, the importance of including the effect of the aerosols in the infrared 
was also emphasized (see ref. F-8). By contributing to the infrared opacity, 
aerosols cause an enhanced greenhouse effect. This acts in the opposite direc- 
tion to the typical effect of aerosols at visible wavelengths, where they 
usually enhance the global albedo. 

In an earlier paper (ref. F-9), initial estimates of the effect of 
aerosols produced by SST emission on the global albedo were presented. This 
paper includes a more detailed description of these results and the infrared 
calculations for estimating the magnitude and sign of the surface temperature 
change. These computations encompass the cumialative effect of all aircraft 
flying in the stratosphere. The term "SST" is used to refer to this case. The 
corresponding calculations for the aerosols produced by Space Shuttle flights 
are also performed. 

These calculations represent an improvement over many of the previous in- 
vestigations because they include more accurate optical properties of the 
aerosols, allow for the effect of the aerosols in the thermal infrared, treat 
accurately gaseous absorption in the visible, and include water clouds in the 
troposphere. The optical properties of the aerosols have been modeled as 
accurately as is presently feasible using measiired optical constants of both 
current and perturbed aerosol components together with a carefully constructed 
model of the other aerosol properties. 

In previous calculations, the optical constants, especially the imaginary 
indices of refraction, either were considered to be free parameters or were 
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set equal to "plausible" values. Since even the sign of the net effect of the 
aerosols on the mean surface temperature depends critically on the values of 
the optical -constants, employment of measured values represents a significant 
advance over earlier effo;M s. 

The particle size distribution functions used are also based on a more 
thorough scrutiny of available empirical data than those used by many earlier 
investigators. The model of the aerosol properties and the sources for their 
optical constants arc described in reference F-10. 

The description of the radiative transfer calculations is followed by a 
discussion of the results. The dependence of various quantities o'f' int.erest 
is plotted as a function of the change in the optical depth of the stratosphere 
at a reference wavelength of 0.55 micrometer. Then, the perturbation optical 
depths that will result from the amount of emission exped.ed from SST’s and 
Space Shuttle vehicles over the next several decades are determined. From a 
comparison of perturbation optical depth values to the results of the calcula- 
tions, the magnitude and importance of the surface temperature change resulting 
from the added aerosols are assessed. 


COMPUTATIONAL TECHNIQTIF 


The radiative transfer calculations consisted of two parts. A numerical, 
multiple-scattering scheme was used to determine the global (spherical) albedo 
and the variation of net solar flux with altitude in the atmosphere. These 
results were subsequently used as boundary conditions in a numerical calcula- 
tion of the infrared thermal radiation that yielded the value of the surface 
temperature. In both calculations, only the global or hemispherical heat 
balance was considered. Thus, the calculations did not allow for latitudinal 
or diurnal variations of the model parameters. In the following paragraphs, 
the numerical technique used to treat the solar radiation is described, the 
sources of the parameters used in this calctilation are summarized, and a 
similar presentation for the thermal computations is given. 


Multiple-Scattering Calculations 

For the visible portion of the spectrum, calculations that considered the 
scattering and absorption of unpolarized light by both aerosols and gases were 
performed. The solution was accomplished with a computer program based on the 
inhomogeneous doubling method described in references F-Q and F-11. According 
to this method, the atmosphere was divided into a series of layers, each one of 
which is homogeneous in its properties but which can differ in its properties 
from any other layer. Thus, allowance was made for the vertical inhomogeneity 
of the atmosphere. An initial step in this computation was determination of 
the single-scattering properties of the aerosols for each layer by using a Mie 
scattering program with allowance being made for the scattering and absorption 
contribution from gases. With this information, the scattering and transmis- 
sion properties of individual layers were determined by means of the ordinary 
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doubling principle (ref. F-12). The surface was modeled as a Lambert reflector 
with a specified reflectivity. Then, a generalization of the doubling method 
was used to combine the separate layers, including the surface. The end resvilt 
of these operations was a specification of the upward, downward, and net fluxes 
at all layer interfaces, including the top of the atmosphere and the surface. 
Integration of the net fl\ixes over solar zenith angle yielded the globally 
averaged value of the solar energy deposition profile, whereas a comparable 
integration of the upward component of the flux at the top of the atmosphere 
resulted in a detemination of the global (spherical) albedo. (These integrals 
are defined in ref. F-U.) The term "global albedo" refers to the total amount 
of solar energy reflected by the Earth to space divided by the amount of solar 
energy incident at the top of the atmosphere. 

This calculation has several attractive features. It employs an accurate 
solution to the multiple-scattering problem that allows for cases in which the 
total optical depth is significantly greater than unity. As a result, opti- 
cally thick water clouds can be included in the lower atmosphere. Allowance 
for gaseous absorption is also made accurately. Ordinarily, inclusion of the 
near-infrared absorption bands of water vapor, CO^, and oxygen presents a 

problem for m\iltiple-scattering calculations. To limit the number of wave- 
length points at which calculations are necessary, one must normally consider 
the absorption that takes place over fairly broad wavelength intervals (~0.1 
to 0.3 micrometer). Such transmission averages of the gas absorption destroy 
its monochromatic nature, and, in particular. Beer's law is no longer obeyed 
by such averages. Thus, it is not correct to employ directly transmission 
averages in mixltiple-scattering calculations. However, in accordance with 
references F-9 and F-11, the usual equivalent width formulas were analyzed 
into a sm of monochromatic transmission factors multiplied by a probability 
coefficient. In essence, this technique is equivalent to transforming inte- 
grals over wavelength into integrals over the absorption coefficient (ref. 

F-13) and allows a small number of monochromatic calculations. Miiltiple- 
scattering calculations were performed for each term of the absorption coef- 
ficient distribution, and the results were weighted by the probability factors. 

A final important characteristic of these calculations was the use of 
measured optical constants for each chemical component of the aerosols. When 
a given layer included several kinds of aerosols, separate Mie scattering 
calculations were performed for each component and the results were added to- 
gether, weighted by their relative cross sections at the wavelength of inter- 
est. The aerosol optical depth of a given layer was found by multiplying the 
aerosol optical depth at a reference wavelength by the ratio of the total aero- 
sol cross section at the wavelength of interest to that at the reference 
wavelength (ref. F-IO). 


Solar Calculation 

The sources of the parameters used in the solar calculation and the nature 
of the model used are described in this section. The solar spectrum from 0.2k 
to h,6 micrometers was divided into 26 wavelength intervals. Separate calcu- 
lations were performed for each wavelength interval; and the resvilts were 
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combined, weighted by the amoxint of solar energy in each interval. In some 
near-infrared intervals that contained gas absorption bands, a number of cal- 
culations were performed, each one corresponding to a given gas absorption 
coefficient. The atmosphere was divided into 15 layers, each 3 kilometers 
(9800 feet) deep. The ground was included as an additional layer and charac- 
terized by a reflectivity of 0.10 (ref. F-?). Water clouds were considered to 
be present over 50 percent of the planet and were located in the layer that 
covered the altitude range from 3 to 6 kilometers (9800 to IQ 600 feet) above 
the ground (ref. F-2). Separate calculations were performed for regions with 
and without water clouds, and the results were combined. 

The properties of the aerosols were specified as follows. The results 
cited in reference F-10 were used to define the aerosol properties in each at- 
mospheric layer for the "base" or unperturbed model. These properties included 
optical depth at the reference wavelength of 0.55 micrometer, chemical composi- 
tion, optical constants, and size distribution function. The properties used 
are summarized in reference F-10. The ma.^or component of the aerosols expected 
to be produced by SST's and by Space Shuttle vehicles are hydrogen sulfate 
(H^SOj^) and aluminum oxide (AI^O^), respectively (ref. F-li* ) . The particle 

size distribution function of the stratospheric H^SOj^ layer under perturbed 

conditions is not known, but some evidence exists that the "base" sulfate dis- 
tribution does not change significantly during epochs of large increases in the 
stratospheric optical depth resulting from volcanic explosions (ref. F-IO). 
Hence, the size distribution of the unperturbed stratosphere was adopted for 
the perturbed component due to SST's. By coincidence, measurements of the 
aerosol size distribution made in the exhaust of prototype Space Shuttle en- 
gines (refs. F-I5 and F-I6), when convoluted with a one-dimensional diffusion- 
sedimentation model (ref. F-17 ) , resulted in a function very similar to the 
"base" stratospheric function. Hence, the "base" distribution was also used 
to describe the pertxirbation caused by Space Shuttle vehicles. In both cases, 
if differences actually exist between the base stratospheric function and the 
perturbed function, the differences will be increased ntunbers of larger par- 
ticles in the perturbed case. This difference would tend to increase the in- 
frared effects of the aerosols and would most probably reduce the temperature 
changes that are found. Thus, the use of the base size distribution function 
is probably a conservative choice in that it may lead to overestimates of the 
net effects. The total optical depth of the two perturbation sources is taken 
as a free parameter in the calculations. In the section entitled "Discussion," 
estimates of this important parajneter are provided from studies of the amount 
of effluent that is expected from SST's and Space Shuttle vehicles during the 
next several decades. When material was added to the stratosphere, the base 
optical depth of each layer of this region was multiplied by the same factor. 
Because the base optical depth of a layer decreases with altitude, aerosols 
were effectively added mainly to the lower portions of the stratosphere. 

Deirmendjian' s (ref. F-I8) cloud model chlorine (Cl) was adopted to 
describe the size distribution function of the water cloud particles; this 
function has a mode radius of It micrometers. The water cloud optical depth at 
a 0.55-micrometer wavelength was chosen to be 6 so that the global albedo of 
the unperttirbed model equaled a value of 0.28, which was determined from 
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satellite measurements (ref. F-19). The optical constants of the water drop- 
lets were obtained from reference F-20. 

The sources of the parameters for the gaseous components of the model are 
described next. Allowance was made for gaseous absorntion caused by 0^, water 

vapor, CO^, and 0^ and for Rayleigh scattering by the air molecules. The ver- 
tical profiles of water vapor and the constant mixing ratios of 0 ^ and CO^ were 
obtained from data given in reference F-21, whereas the 0^ profile was obtained 

from Whitten's results (private communication). The water vapor abundances 
were converted to relative hiunidity profiles since, in many of the calcula- 
tions, the reasonable assumption was made that this variable would remain con- 
stant with changing surface temperature. The equivalent width parameters and 
formulas for the major water vapor and CO^ bands were adopted from reference 

F-22, whereas similar data for the major 0^ bands were obtained from reference 

F-23. These data were analyzed to determine the probability distribution 
functions of monochromatic absorption coefficients for each major band region 
of a given gaseous constituent. Lindzen and Will's (ref. F-2i* ) square wave 
model of the Hartley and Chappuis bands of 0^ and their absorption coefficients 

were used at four gauss-point, wavelength locations for the Huggins band. 
Finally, scattering coefficients given in reference F-25 were used to compute 
the Rayleigh scattering optical depths. The data in reference F-25 include 

_li 

an allowance for slight departures from a X dependence, where X is wave- 
length. 

The objective of infrared thermal calculations was to obtain an estimate 
of the change in the mean s\irface temperature caused by a given stratospheric 
aerosol perturbation. In general, this task was accomplished by performing an 
energy balance calculation at the top of the atmosphere.. However, in a few 
cases, this calculation was supplemented by one in which the entire vertical 
temperature structure was foiind under conditions of radiative-convective 
equilibrium. In the first type of calculation, the thermal infrared flux at 
the top of the atmosphere was determined for a sequence of surface tempera- 
tures, and the results were interpolated to find the surface temperature for 
which the thermal energy radiated to space equaled the solar energy absorbed 
by the planet. The formulas used to calculate the thermal infrared flux are 
given in reference F-26. These equations involve a reformulation of the formal 
solution of the transfer equation in terms of transmission average functions, 
with the source function set equal to the Planck function. The amount of 
solar energy absorbed was obtained by interpolation from calculations made with 
the solar energy program at two values of the surface temperature. Thus, 
allowajice was made for the smaJ.1, but not insignificant, dependence of the 
solar component on temperatiire. This dependence arises from changes in the 
absolute amount of water vapor in the atmosphere under conditions of constant 
relative humidity. 

A few calculations of the vertical temperature structure were performed 
under conditions of radiative-convective equilibrium using the flux iterative 
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method cited in reference F-PJ , The results of the solar calculations were 
used to specify the desired net infrared fluxes for this set of calculations. 
The iteration procedure provided a method to determine the vertical tempera- 
ture structure for which the net infrared flux at each layer in the atmosphere 
was ba.lanced by the net solar fliix at that location. Such a state is one in 
which pure radiative equilibrium holds. However, when the local lapse rate 
exceeded a critical value at which convection commenced, the lanse rate was 
readjusted to equal this limitinp, convective lapse rate. Unfortunately, the 
iteration scheme was not fully convergent, and it led to mean fractional flux 
residuals of approximately 1 to 2 percent after a number of iterations. Thus, 
this type of calculation was applied to cases involving larpie optical depth 
perturbations to determine the general behavior of the stratospheric tempera- 
tures, which may follow quite a different trend from that of the surface and 
tropospheric temperatures. 

Many infrared model parameters, such as aerosol properties, have already 
been described. The infrared was divided into spectral intervals, which 

totally covered the spectral region from wave numbers 25 to 3000 cm ^ (wave- 
lengths i^OO to 3.3 micrometers). The 25- to 50-cm ^ interval was divided into 
a single region, the 50- to 1300-cm ^ region into 50-cm ^ segments, and the 

1300- to 3000-cm ^ interval into 100-cm ^ partitions. Within each spectral 
interval, the transmission behavior of the gases was modeled by using an expo- 
nential of a broad-band opacity variable, which in turn followed a power law 
dependence on gas amount, total effective pressure, and temperature (ref. 

F-26). The value of the power law exponents and the proportionality constant 
for each gas species and wavelength were obtained by fitting this representa- 
tion to a combination of laboratory and theoretical estimates of the trans- 
mission properties of water vapor, CO^, and 0^ at conditions relevant for the 

Earth's troposphere and stratosphere (refs. F-21 and F-28 to F-32). 

The infrared aerosol optical depth was found for each layer of the atmos- 
phere and each wavelength interval by multiplying the aerosol optical depth at 
0.55 micrometer by the ratio of the aerosol absorption cross section of the 
wavelength of interest to the total cross section at the reference wavelength 
(ref. F-IO). This ratio was determined by a Mie scattering program. The 
neglect of scattering by the aerosols at infrared wavelengths was found to be 
reasonable since the single-scattering albedo was usually quite low. A low, 
single-scattering albedo was the result of the large imaginary index of refrac- 
tion and the small ratio of the mean particle size to the wavelength in the 
infrared part of the spectrum. The optical constants for most materials are 
not known as far as itOO micrometers. The values were used at the largest wave- 
length at which they are known to specify the optical constants beyond this 
point . 


The starting point for the temperature/pressure structure of the atmos- 
phere was the U.S. standard atmosphere (ref. F-33) for midlatitude spring and 
fall conditions. In the energy balance calcxxlations, the temperature struc- 
ture of the entire vertical column was generally changed in response to a 
change in sxirface temperature from the U.S. standard atmosphere values in such 
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a way that the lapse rates were kept constant. However, in some calculations, 
the stratospheric temperatiores were kept constant, and only the tropospheric 
temperatures were changed under the constraint of a constant lapse rate. In 
the flux iterative calculations, the convective lapse rate was set equal to 
6.5 K/km, a value characteristic of the present troposphere. 


RADIATIVE TRANSFER CALCULATIONS 


In this section, the radiative transfer calculations are summarized by 
showing the manner in which various quantities of interest, such as global 
albedo and surface temperature, depend on the enhancement of the stratospheric 
aerosol optical depth At at the reference wavelength of 0.55 micrometer. 
Separate calculations were performed for Al^O^ and H^SOj^ particles, which are 

produced by Space Shuttle and SST vehicles, respectively. Because only small 
values of At are of interest for each case, the cumulative effect of both 
vehicles is additive to first approximation. 

The results of the calculations in the visible portion of the spectrum 
are as follows. The variation of the global albedo with At for the two 
sources of added particles is shown in figure F-1. In both cases, the added 
aerosol burden in the stratosphere leads to an increase in the global albedo, 
which, alone, will cause a cooling of the Earth. For an equal At, the 
Al^O^ particles cause a larger increase in the Earth's albedo than do the 

particles. This difference is a result of the Al^O^ aerosols having a 

larger value for the real part of their refractive index (~1.T7 compared to 
1.1+3) and, therefore, backscattering a larger fraction of the incident solar 
radiation. 

Next, the change in the vertical distribution of the absorbed solar 
energy is examined. The At dependence of the fraction of energy absorbed 
in the stratosphere is shown in figure F-?., Here, as elsewhere in the dis- 
cussion, fractional absorption refers to the ratio of the amount of solar 
energy absorbed to that incident at the top of the atmosphere. In contrast 
to the trend for the Earth as a whole, an increase in At leads to an increase 
in the amount of solar energy deposited in the stratosphere. 

Part of the reason for this trend in the stratospheric heating is the 
absorption of solar radiation by the particles. In fig+ire F-3, the fractional 
absorption by aerosols is plotted as a function of At. It can be seen that 
as the number of particles is increased, the amount of energy they 

absorb is increased. However, almost no variation is shown by the Al^O^ 

particles. In the added AlgO^ particles, the total aerosol heating is 

dominated by. the H^SOj^ particles of the unperturbed stratosphere that are 

also present. The particles absorb solar energy principally in their 

near-infrared bands. In contrast, the imaginary index of refraction, which 
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determines the amount of solar energy absorption, remains very small for the 
AlgO^ particles throughout the entire visible and near- infrared spectral re- 
gion. It can be seen that the increase in the total fractional heating of the 
stratosphere with an increase in the burden is due in part to energy 

absorption by the aerosols. In addition, as a result of scattering of the 
incident sunlight by both types of aerosols, an increase occurs in the amount 
of sunlight absorbed by the 0^ in the stratosphere. The 0^ absorntion accounts 

for almost all of the trend of the total stratospheric absorption shown in 
figiire F-2 for the AI 2 O 2 particles and about half of the trend for the H^SOj^ 
particles . 

In figure F-^t, the effect of H^SOj^ aerosol enhancement on various levels 

in the stratosphere is examined. The ordinate is the fractional change in 
the solar heating at a given level and is referenced to the amount of solar 
heating at that level in the unperturbed model. The two curves shown cor- 
respond to two values of the optical depth perturbation At. As expected, 
the greatest changes occur in the lowest portions of the stratosphere, where 
most of the added particles are located. Analogous results are found for the 
Al^O^ particles, although at the highest pressure levels in the stratosphere, 

a decreased heating is found compared to that for the unperturbed model. 

In figure F-5, the results for the fractional absorption below the tro- 
popause are shown. This quantity refers to the sTim of the energy absorbed 
in the troposphere and at the ground. In accordance with the trend found for 
the global albedo in figure F-1, it can be seen that the amount of energy 
absorbed below the tropopause decreases with increasing aerosol burden in the 
stratosphere. Alone, this result implies a cooling trend for the troposphere 
and the surface. However, a judgment on this issue must be deferred until 
the effect of the aerosol opacity in the thermal infrared region is included. 

Besides affecting the global energy budget, aerosols added to the 
stratosphere can also upset the 0^ photochemical cycle by changing the rate 

at which 0^ molecules are dissociated by incident sunlight. In figure F-6, 

the fractional change in the 0^ destruction rate for the entire atmospheric 

column is shown as a function of At. The 0^ destruction rate was calculated 

by summing the number of photons absorbed by 0^ in each of its three major 

absorption bands. The number of photons absorbed by 0^ at a given wavelength 

and in a given atmospheric layer was foiind in a straightforward fashion from 
the amount of solar energy absorbed by that layer and the fractional con- 
tribution of 0^ to the absorption coefficient of the layer. An increase in 

the aerosol content of the stratosphere causes an increase in the rate at 
which 0^ is destroyed and, therefore, a decrease in the steady-state abundance 

of 0^. However, even for rather large values of At, the magnitude of the 
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change is surprisingly small both in absolute value and in comparison with the 
0^ depletions expected as the result of the oxides of nitrogen and Cl compounds 

released by SST's and Space Shuttle vehicles, respectively. The 0^ destruction 

rate is only slightly affected by the added aerosols because much of the 0^ 

photodissociation occurs above the levels at which most of the aerosols are 
found. 


The results of the infrared calculations, which provide estimates of the 
effect of the added aerosols on the mean surface temperature, are described 
next. Most calculations of this type involved balancing the radiative energy 
at the top of the atmosphere. First, such a calculation was performed for the 
base or unperturbed global average model, and a surface temperature of 292.7 K 
was obtained. This value is approximately U.5 K higher than the ground temper- 
ature given in reference F-33 for spring/fall conditions at middle latitudes. 

A global average model is not the same as a midlatitude model. This tem- 
perature difference arises partly because the ground temperature varies with 
latitude for the real Earth and partly because the infrared radiation from 
various latitudes is weighted by a nonlinear, approximate temperature to the 
fourth power, relationship. Subsequent calculations provided values of the 
surface temperatiire for pertvirbed cases, which were compared with the result 
for the base case. The basic parameter derived from this comparison is AT, 
which is the difference in the value of the surface temperature between the 
perturbed and unperturbed cases . 

In figure F-T, the variation of Ax and AT is shown for cases involving 
added H2S0|^ particles resulting from SST's. Calculations were performed for 

four sets of conditions, which are represented by the four curves of figure 
F-T. Curves A and B were determined for conditions of constant relative and 
constant absolute humidity. For these calculations and for those represented 
by curve D, the lapse rates were assumed to be invariant throughout the at- 
mosphere. Curve C represents cases in which both the relative humidity pro- 
file and the stratospheric temperatxires were held constant. Finally, to 
determine the importance of the aerosols opacity in the thermal infrared 
region, curve D was computed. Curve D was calculated under the same set of 
circumstances as curve A, except that the aerosol opacity was held fixed to 
its value in the base case. 

For all four cases shown in figure F-7, an enhanced aerosol burden 
causes a decrease in the mean surface temperature. This result indicates 
that the effect of the aeorsols on the global albedo is more important than 
their effect in the infrared. In turn, the relative importance of the 
aerosols in these two wavelength regions is a reflection of their very small 
mean size. As shown in figure F-5 (ref. F-IO), the stratospheric aerosols 
have a much larger optical depth at visible wavelengths than at longer wave- 
lengths. A comparison of curves A and D shows that the magnitude of the AT 
values is decreased by approximately i+0 percent when allowance is made for the 
infrared opacity of the aerosols. Thus, although this effect is not the 
dominant one, it cannot be neglected. 
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Because curves A and C are fairly close, thf’ resuits are not very sensi- 
tive to the assumed change in the stratospher i c t.emperature profile. Calcu- 
lations performed for conditions of radiative-convective equilibrium for a At 
value of suggest that curve A provider, a more realistic estimate than curve C 
provides. This result indicates that the added aerosols cool the strato- 
sphere to a greater degree by radiating in the thermal infrared than they warm 
this region by causing a slight enhancement in the amoxmt of absorbed solar 
radiation. Finally, a comparison of curves A and B shows that larger AT 
values are achieved under conditions of constant relative humidity than under 
conditions of constant absolute humidity because, for the constant relative 
humidity case, the amount of water in the atmosphere decreases with decreas- 
ing temperature. In turn, less water vapor implies less solar energy absorp- 
tion in its near-infrared bands and a smaller greenhouse effect in the thermal 
infrared. The assiamption of constant relative humidity is probably the more 
realistic one. 

Analogous curves are shown in figure F-8 for the Al^O^ particles produced 

by Space Shuttle vehicles. Again, a net cooling of the Earth is found for 
all fouir cases. However, the infrared opacity cancels a good fraction of the 
effect of the particles in the visible region. The infrared effect is more 
powerful for the Al^O^ particles than for the particles because the 

imaginary index of refraction is larger for the fornner than for the latter. 

An enhanced imaginary index of refraction in the infrared results in an en- 
hanced absorption cross section. 


DISCUSSION 


In the preceeding section, estimates of the decrease in the mean surface 
temperature AT were obtained as a function of the optical depth perturbation 
At caused by H^SOj^ and Al^O^ aerosols introduced into the stratosphere. The 

results of these calculations are summarized in figures F-T and F-8, in which 
AT is plotted as a function of At. Besides applying to the global average 
situation, these results can also be used for optical depth perturbations 
occurring in a single hemisphere. (The angular averages over solar zenith 
angle are identical in these two cases.) In this section, the values of AT 
that have characterized epochs of major climate changes in the past are dis- 
cussed so as to estimate a lower bound on AT, below which no serious climatic 
impact is expected. Then, estimates of the magnitude of Ax expected from 
SST’s and Space Shuttle vehicles operating during the next several decades 
are provided, and the importance of their effect on the climate is determined. 

The significance of a given value of AT can be assessed by comparing it \ 

with values associated with major climatic events of the past. During the 
past several million years, the climate of the Earth has undergone a mamber 
of major fluctuations, in which the climate has alternated between periods of 
major glaciation and relative thaw. At the present time, of course, the 
Earth is within one of the relatively warm epochs, although frequent variations 
with a smaller amplitude have occurred over the current interglacial period. 


F-12 


BilllllilfillKlIIIllL 


The difference between the cxorrent mean surface temperature and that estimated 
for the peak of one of the major ice ages is approximately 5 K (ref. F-3^+). 

Over the last 1000 years, there have been extended periods of time when the 
mean svirface temperatiire was approximately 1 K less than its current vadue 
(refo F-3^). This apparently small difference had major consequences in some 
cases. For example, during some of these colder epochs, the agricultural 
yield was substantially reduced in countries at high geographical latitude, 
and major famines resulted (ref. F-3U). 

This discussion indicates that a temperature change as large as 1 K could 
have serious consequences. However, values that are a factor of 2 or 3 
smaller are probably on the borderline of being important. The numerical 
calculations of AT represent determinations of the change in surface tem- 
peratiire that result from radiative effects. Because many feedback mechanisms 
were not considered, the actual value of AT might be somewhat larger. For 
example, during colder periods, the amount of surface covered by frost and 
ice increases. This acts as a positive feedback effect by increasing the 
global albedo. Calculations (refs, F-35 and F-36) suggest that the net effect 
of a number of feedback mechanisms is to increase the calculated values of 
AT by a factor of approximately 3. For these reasons, 0.1 K was adopted as a 
threshold value for AT and is believed to be a conservative estimate. Smaller 
calculated values of AT are not expected to change the climate seriously. 

To obtain estimates of the optical depth perturbation At that can be 
expected from SST and Space Shuttle vehicles, the mass of material that will 
be emitted into the stratosphere from these vehicles over the next several 
decades was first estimated. The steady-state mass of aerosols produced in 
one hemisphere by the continued operation of a given type of vehicle is re- 
lated to the emission rate of the vehicle according to 


m 


H 


e.f. — t.c.n. 
1 1 M 111 


(F-1) 


where 


m is the steady-state mass of aerosols 

H is the fraction of the aerosols deposited in a given hemisphere 

e. is the total amo\int of mass emitted by a single vehicle per unit 

time 

f . is the fraction of the total mass that contains the aerosol- 

producing component 

M .is the moleciilar weight of the aerosol 
8 . 

is the molecular weight of the aerosol-producing component 
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is the residence tine in the stratosphere 

c. is the efficiency of convertinR the aero sol -producinp component 
into an aerosol 

n^ is the number of vehicles operating in both hemispheres 
refers to the various altitude ranges in the stratosphere 

Estimates were obtained of the steady-state mass of and A.l^Q^ 

particles produced by SST and Space Shuttle vehicles, respectively, by using 
equation (F-l) in combination with the data given in table F-1. These data 
have been assembled from information given in references F-lh and F-IT. The 
second column in table F-I provides estimates of t^^. The third column gives 

values for c^ for H^SOj^ aerosols produced from the sulfur dioxide (SO^) 

emitted by SST and other high-altitude aircraft. The next four cf)liunns of 
the table give estimates of the SO^ emission rates caused by these types of ve- 
hicles for several time frames. Both nominal and upper limit values are given. 

These emission rates are equivalent to the product (e.n.f.) in equation (F-1). 

_3 111 

These values were derived by using a value of 10 for f, the fraction of 

the total mass that is S0„, and M /M =1.6 for an H_S0, solution that is 

2’ a e 2k 

75 percent H^SOj^ by weight (ref. F-3fi). The final column of table F-I provides 

values of the emission rate of Al2*^3 P^i^tides by the Space Shuttle vehicles 

expected to be operating in 1990. This quantity is equivalent to the product 
(e.n.f.) in equation (F-l), In this case, M /M =1 and c. = 1. Finally, 

1 1 X S. 0 X 

although almost all of the emission takes place in the Northern Hemisphere, 
meridional transport carries approximately 30 percent of the material into the 
Southern Hemisphere. Hence, H = 0.7 for the Northern Hemisphere and H = 0.3 
for the Southern Hemisphere. 

Next, the mass of aerosols m found from equation (F-l) and table F-I 
was converted into an optical depth perturbation Ax. The relationship between 
these two quantities is given by 

2 

/„ Q , ira N(a)da 
m ,/0 ext 

^ ■ Ap (F-2) 

Jo ~ a%(a)da 


where A is the area of a hemisphere 

p is the density of the aerosol 

Q . is the efficiency factor for scattering and absorbing incident 
radiation (ref. F-IO) 
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a is particle radius 

N(a) is the size distribution function of the aerosols 

The integrals appearing in equation (F-P) were evaliiated by using the strato- 
spheric size distribution function (ref. F-IO) in a Mie scattering computer 
program. 

If the aerosols are distributed ovei- a scale height of 10 kilometers 
(32 000 feet), the following numerical relationships exist between Ax 
and m. 


Ax = 0.0503 m. 

so^ 

(F-3a) 

Ax = 0.0310 m. 


(F-3b) 

Ax = 0.0228 m. 

AI2O3 

(F-3c) 


where m is the mass in units of micrograms per cubic meter. The compound 
to which m refers is indicated to the right of the comma in these equations. 
For example, in equation (F-3a), m is the mass of SOp that is converted into 

HpS0|^ particles, and Ax is the optical depth of these aerosols. 

Table F-II is a simimary of the values of Ax at a wavelength of 0.55 
micrometer, as determined from equation (F-2). The SST's cause about 70 
percent of the steady-state burden of H^SO^ particles produced by high-altitude 

aircraft. Although planes flying lower in the stratosphere generate more 
emission than SST*s, the effect of these other aircraft is minimized by the 
short residence time of their products. For the Al^O^ particles produced by 

the Space Shuttle vehicle, the amount of mass found from table F-I was doubled 
to allow for particles emitted at altitudes above 21 kilometers (68 000 feet); 
then, the Ax so found was halved to allow for the rapid loss of the larger 
particles. About half the mass of the P^o<iiJced by the Space Shuttle 

exhaust is contained in particles larger than a micrometer; these particles 
are quickly removed from the stratosphere by gravitational sedimentation. In 
all cases, the optical depth perturbations are quite small; in fact, none of 
them is as large as the perturbation produced throughout the mid-1960's by 
small-scale volcanic activity (as discussed in ref. F-IO). 

The impact of the aerosols produced by SST and Space Shuttle vehicles, 
is as follows. According to figures F-T and F-8, optical depth perturbations 
of approximately 0.01 and 0.015 are required for SST and Space Shuttle, 
respectively, to cause a decrease in the mean surface temperature equal to 
0.1 K. This value for AT is the threshold figure discussed previously. All 
the optical depth perturbations given in table F-II are less than these 
critical values. From this comparison, it can be concluded that SST and 
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Space Shuttle vehicles operatinp, at the levels estimated for the next several 
decades will not cause a serious change in the climate as a result of the 
aerosols they add to the stratosphere. 

This conclusion on the impact of SST and Space Shuttle, vehicles can be 
put into some perspective by estimating the number of vehicles operating on 
a regular basis required to produce the critical value of AT. Since the 
nominal 1990 case treated in table F-II was based on 2S0 SST’Vs flying every 
day for I 4 .S hours in the stratosphere and since these vehicles contributed 
approximately 70 percent of the Ax attributable to high-altitude aircraft, 
approximately 7500 SST's are required to generate the critical value of AT. 
The Space Shuttle case treated in table F-II involved about one launch per 
week. To raise AT to its critical figure would require approximately 70 
launches per day. The critical number for the Space Shuttle is far in excess 
of the number of vehicles that are likely to be operating in the near future. 
Thus, the climatic effect of these vehicles is unlikely to be changed by 
reasonable modifications to the parameters used in these calculations. How- 
ever, for high-altitude aircraft, the projections for the rnamber of SST's 
and other aircraft operating in the year 2000 are sufficiently near the criti- 
cal numbers that estimates of AT for such aircraft should be updated as new 
projections become available. 


CONCLUSIONS 


Calculations of the effect on the global heat budget of adding aerosols to 
the stratosphere have been presented. The global albedo increases and the 
mean surface temperature decreases as a result of increases in the number of 
H 2 S 0 j^ and AI 2 O 2 particles that are produced by SST and Space Shuttle vehicles, 

respectively. The effect of the aerosols on thermal radiation was important; 
furthermore, the thermal radiation effect partly counteracted the cooling of 
the surface caused by the effect of the aerosols at visible wavelengths. The 
results of these calculations' on the decrease in surface temperature as a 
function of optical depth perturbation are summarized in figures F-7 and F- 8 . 

These results were used to assess the significance of the climatic impact 
of aerosols produced by SST's, other high-altitude aircraft, and Space Shuttle 
vehicles operating diiring the next several decades. For this task, the AT 
values associated with past major changes in the climate were used to infer 
tha' AT values below 0.1 K would not have major consequences. Next, the 
emission rates estimated for the vehicles of interest over the next several 
decades were used to determine the Ax values summarized in table F-II. 

After the data were combined, it was concluded that the aerosols produced by 
SST's, other high-altitude aircraft, and Space Shuttle vehicles operating during 
the next several decades would not seriously alter the climate. Very large 
numbers of these vehicles are required before threshold conditions are reached: 
approximately 7500 SST's and approximately JO Space Shuttle launches per day. 

It seems highly unlikely that the conclusions about the impact of Space 
Shuttle vehicles will be altered by changes in their projected nimiber or by 
changes in their exhaust characteristics. However, the same statement cannot 
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be made for SST's and other high-altitude aircraft. As new projections of their 
numbers and other characteristics become available, the estimates of Ax 
attributable to such aircraft should be revised and their impact reevaluated. 
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TABLE F-I.- FACTORS USED TO CALCULATE THE STEADY-STATE MASS OF 
AEROSOLS INTRODUCED INTO THE STRATOSPHERE BY VARIOUS VEHICLES 


Factor 

SST and other high-altitude aircraft 

Space Shuttle 
1990 nominal 
AI2O2 emission, 

rate, kg/yr 

Altitude range, km (ft) 

Residence 
time, yr 

Conversion 

efficiency 

SO2 

emission rate, kg/yr. 

in - 

1990 

nominal 

1990 

max. 

2000 

nominal 

2000 

max. 

9 to 12 (26 000 to 36 000) 

0.160 

0.810 

1*. 37x10^ 

8.16x10’^ 

7.50x10"^ 

0 

1.939x10 

- - 

12 to 15 (36 000 to 1*7 000) 

.1*51* 

VO 

0 

CO 

• 

1.53x10’^ 

2.85x10’^ 

2.98x10"^ 

7.61x10'^ 

l.l*xl0® 

15 to 18 (1*7 000 to 58 000) 

.978 

.869 

7.0x10^ 

1.1*X10^ 

1.8x10^ 

5.6x10® 

1.2x10® 

18 to 21 (58 000 to 68 000) 

1.715 

.91*5 

1.1*6x10^ 

3.7l*xl0^ 

1. 10x10® 

1.10x10® 

1.1x10® 


I 

ro 



I 


TABLE F-II.- ESTIMATES OF THE STRATOSPHERIC OPTICAL DEPTH 
PERTURBATION CAUSED BY VARIOUS VEHICLES 
[Northern Hemisphere] 


Vehicle 

Time frame 

At 

SST*s and other high- 

1990 nominal 

i4,9x]0"^ 

altitude aircraft 

1990 max. 

8,9><10“^ 


2000 nominal 

1,1x10"^ 


2000 max. 

S.B’^IO"^ 

Space Shuttle 

1990 

^3.3x10"^ 


Values for the Southern Hemisphere can be obtained by 
multiplying the tabulated numbers by three-sevenths o 

^Adjusted as described in-text. 


F-22 




JI I I 1 £ L 










Figure F-3.- Dependence of fractional aerosol absorption in the stratosphere 
on optical depth perturbation of H^SO. and A1„0^ aerosols. 
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Figure F-U.- Fractional change in solar heating as a function of pressure in 

the stratosphere for SST-produced H2S0|^. 
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Figure F-5.- Fractional absorption below the tropopause as a function of aerosol 
optical depth, perturbation of H_SO,. and Al^O,, aerosols. 
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O Curve A = Constant relative humidity 

□ Curve B - Constant absolute humidity 

O Curve C = Constant relative humidity and ^ 

stratosphere temperature 

A Curve D = Constant relative humidity and infrared aerosol 
opacity 


Figxire F-T.- Change in mean surface temperature AT as a function of optical 
depth perturbation of SST-produced H^SO^ aerosols. 
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Figure F-8.- Change in mean surface temperature as a function of optical 
depth perturbation of Space Shuttle Al^O^ aerosols. 
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APPENDIX G 


NO DEPOSITION IN THE STRATOSPHERE FROM THE 

X 


By H. S. 


SPACE SHUTTLE ROCKET MOTORS 

3. I 3 C 

Pergament, R. I. Gomberg, and I. G. Poppoff 


APPROACH 


The determination of the amount of oxides of nitrogen (NO^) deposited in 

the stratosphere by the Space Shuttle solid rocket motors has been based pri- 
marily on an extensive series of calculations using state-of-the-art computer 
codes to determine both the mass flow of NO^ leaving the nozzle and the mass 

of NO^ produced in the Space Shuttle plume as a resiilt of shocks and afterburn- 
ing. This is the first attempt to predict total NO^ deposition rates in the 
stratosphere using the new calciilation procedures described in references G-1 
and G-2 and elsewhere. Results are reported (refs. G-1 and G-2) for the base- 
line propellant at 30-kilometer (95 000 foot) altitude; calculations were also 
made at 15-kilometer (1+7 000 foot) altitude using a similar procedure. Es- 
timates of the NO deposition rate at 15 kilometers (i*T 000 feet) and 30 kilo- 

X 

meters (95 000 feet) for the alternate "A” propellant have been made using a 
more simplified (and conservative) procedure than the one described in refer- 

2 

ences G-1 and G-2 and elsewhere. It should be stressed that, because both 
NO and hydrogen chloride (HCl) participate in the ozone (0_) destruction cycle 

(NO^ less efficiently than HCl), the primary interest is in the ratio of NO^ to 

HCl deposition rates in the stratosphere. 

Experimental verification of these calculations is limited; only a single 
field measxirement of nitric oxide (NO) in the exhaust cloud of a Titan III 
booster has been teiken at an l8-kilometer (58 000 foot) altitude. The NO^ 

measurements have also been made within Titan III ground clouds. The useful- 
ness of these data in verifying the predicted NO deposition rates is dis- 
cussed in the following sections. 


& 

AeroChem Research Laboratories, Inc. 

^NASA Langley Research Center. 

^NASA Ames Research Center. 
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R. D. Thorpe, H. S. Pergament, and B. Hwang: NO^ Deposition in the 

Stratosphere by the Space Shuttle Solid Rocket Motors. Proceedings of the 9th 
JANNAP Plume Technology Meeting, in press. 
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CALCULATED RESULTS FOR BASELINE PROPELLANT 


The total mass flow of NO^ leaving the Space Shuttle plume and deposited 
in the stratosphere is determined by first calculating the amount of NO^ exit- 
ing from the solid rocket motors and by then calculating the enhancement in 
NO^ caused by plimie shocks and afterburning of the hydrogen (H^) and carbon 

monoxide (CO) in the exhaust products with ambient air. Because of the rela- 
tively slow rate of NO^ depletion in the nozzle, as determined by a nonequi- 
librium chemistry analysis, the mass fraction of NO^ at the nozzle exit plane 

is considerably higher than that resulting from standard thermochemical equi- 
librium calculations. Table 5-1 gives the mass flow of NO^ leaving the two 

nozzles in grams per second and in parts per million by volume; by this means, 

the NO mass flow is related to the total propellant mass flow. The results 
X 2 

reported (refs. G-1 to G-3 and elsewhere ) show that shocks and afterburning will 
enhance the mass flow of NO^ by a factor of 1.6 at 30 kilometers (95 000 feet). 

A similar analysis for 15 kilometers (1*7 000 feet) shows the enhancement factor 
to be 1*.8. This increase with decreasing altitude is due both to a greater 
percentage of mass flow passing through the strong plume shock and to the 
higher temperatiires produced in the afterburning region of the plume. The 
total mass of NO^ deposited per second at 15 kilometers (1+7 000 feet) and 30 

kilometers (95 000 feet) is also shown in table 5-1. 


COMPARISON OF NO DEPOSITION FROM BASELINE AND ALTERNATE PROPELLANTS 

X 


The procedure used to evaluate the NO^ deposition rate from the alternate 

propellant was to evaluate the manner in which differences between the alter- 
nate and baseline propellants could influence both the mass flow of NO^ leaving 

the nozzle and the plume enhancement. Because of the large amount of ammoniimi 
nitrate in the alternate propellant, the possibility that significantly large 
nonequilibrim levels of NO^ could be produced in the combustion chamber was 

investigated. Nonequilibrium chemical calculations were made at nominal cham- 
ber conditions, assumi ng (for conservatism) that 25 percent of the N in the 
propellant was initially converted to NO. (At equilibrium, only 0.1 percent 
of the N in the propellant is converted to NO. ) The results of these nonequi- 
librivim chemical calciilations indicate that at the chamber temperature of 
2690 K, the decay of NO^ to its equilibrium value occurs in approximately 10 


2 

R. D. Thorpe, H. S. Pergament, and B. Hwang: NO^ Deposition in the 

Stratosphere by the Space Shuttle Solid Rocket Motors. Proceedings of the 9th 
JANNAF Plume Technology Meeting, in press. 
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milliseconds, which is considerably less than typical combustion chamber resi- 
dence times. Thus, the alternate propellant NO mass fraction will reach equi- 
librium value in the chamber at a factor of approximately 50 less than the cal- 
culated chamber value for the baseline propellant. 

The other major difference between the two propellants is that the alter- 
nate propellant contains approximately 15 percent (molar) more fuel (Hg + CO) 

in the exhaust products. This increase in fuel will result in more intense 
afterburning. The effects of shocks in enhancing the NO, mass flow in the 
plume should be about the same for each propellant. 

Table 5-H gives the results of NO^ deposition calculations for the al- 
ternate propellant, compared with the corresponding values for the baseline 
propellant. The combined effects of shocks and afterburning result in larger 
plume NO^ production rates (and thus a larger plimie enhancement) for the alter- 
nate propellant. These results demonstrate that the nozzle exit plane mass 
flow is not necessarily a good indicator of the total NO^ deposition rate. 


INTERPRETATION OF FIELD MEASUREMENTS 


Very limited field measurements exist with which to compare the NO^ depo- 
sition rates predicted by the model. A single set of flight data was obtained 
by flying a U-2 aircraft equipped with an NO measurement device through the 
cloud trailing a. Titan III launch vehicle at an altitude of l 8 kilometers 
(58 000 feet). In addition, considerable data are available on the total 
amount of NO^ contained within ground clouds following Titan III launches. 


Stratospheric Measvirements 

The Titan III is the largest solid rocket using the same fuel as the Space 
Shuttle solid rocket boosters and, except for the lack of the parallel hydrogen- 
and oxygen- fueled rocket motor exhaust from the upper stage, can serve as a 
model for Space Shuttle launches. Titan III thrust is Uo percent of the Space 
Shuttle solid booster motor- thrust. Several attempts have been made to deter- 
mine stratospheric exhaust trail properties of Titan III-C rockets launched at 
the Western Test Range. These measurements were made using sampling instru- 
ments carried onboard high-altitude aircraft. In April 197^, a successful 
meas\irement was made in a Titan III exhaust trail at an altitude of I 8 kilo- 
meters (58 000 feet) some 13 minutes after launch by Loewenstein and Savage 
of the NASA Ames Research Center (ARC) using equipment carried onboard the 
ARC U-2. The instriiments measured a reduction of 0^ greater than 4o percent 

below the background level and an average of 12 peirts per billion by voliime 
of NO, approximately 30 times the ambient level. The NO instrument works on 
the principle of chemilminescent reaction with an onboard 0^ reactant gas; 
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the measiirement must be treated as an upper limit since it is not known whether 
other exhaust trail species have residual chemiluminescent reactions. 

The duration of the NO signal during exhaust trail penetration by the air- 
craft suggested that the horizontal dimension was l 600 meters (5200 feet). 
Ground-based photographs, on the other hand, showed that the exhaust trail 
had irregular boundaries, moved in different directions in various altitude 
regions because of winds aloft, and appeared to have a horizontal dimension 
of approximately 3000 meters ( 98 OO feet) at the altitude and time of aircraft 
penetration. 

A major problem in the interpretation of the data in terms of the NO dep- 
osition rate is determining the amo\int of conversion of NO to nitrogen per- 
oxide during the 13 minutes in which the exhaust cloud mixed with the ambient 
air before the measurements were made. This problem, in conjunction with the 
uncertainty of the other chemiliaminescent reactions noted, makes interpreta- 
tion of the measured value in terms of an NO deposition rate excessively 
difficult. 


Ground Cloud Measurements 

A ro\Agh indicator of the adequacy of the chemistry model in the strato- 
spheric calculations can be obtained by a comparison between NO^ data obtained 

in Titan III ground clouds and calculations with a simplified rocket plume 
code which used the identical chemistry model. The results of this comparison 
show that the predicted NO^ deposition is within 30 percent of the measurements, 

an indication that the overall chemistry model used to predict NO^ production 

rates in afterb\irning Space Shuttle pliomes is reasonably accurate. 


DEPOSITION RATES FOR NO and HCl 

X 


A comparison between NO^ and HCl deposition rates for both the baseline 

and alternate propellants is given in table 5-IH. The results show that the 
amovint of NO^ deposited in -the stratosphere is negligibly small compared to 

the amoiint of HCl. In fact, even if the NO^ deposition rate predictions are 

low by as much as an order of magnitude (an unlikely possibility), the levels 
woiild still be negligible for the baseline propellant and the alternate pro- 
pellant at 30 kilometers (95 000 feet). 
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APPENDIX H 


MEDICAL/BIOLOGICAL EFFECTS OF INCREASED uv RADIATION 
RESULTING FROM 0^ REDUCTION 
By D. S. Nachtwey®' 


Space Shuttle operations with the baseline solid rocket motor are 
estimated to reduce the stratospheric ozone (O^) column thickness by 0.2 

percent with an uncertainty of a factor of 3. A percentage reduction in 
the 0^ column can be converted by a series of assumptions and calculations 

into a rough estimate of a percentage increase in biologically harmful ultra- 
violet (BHuv) radiation from the Sun. The rule of thumb that emerges from 
the asstimptions and calculations is that for 0^ reductions of less than 

5 percent, a 1-percent decrease in 0^ colmn thickness represents a 2-percent 

increase in BHuv (refs. H-1 to H-3). The 1-percent /2 -percent rule of thumb 
(also referred to as an optical amplification factor of 2) is based on cal- 
culations involving the increased amount of shorter wavelength uv radiations 
(the band of 285 to 320 nanometers designated uv-B) penetrating through a 
reduced 0^ layer and the increased biological effectiveness of those shorter 

wavelength radiations. 

Not all wavelengths of uv are equally effective in producing biological 
effects; the shorter wavelengths are generally more effective. The relation- 
ship of biological effectiveness to wavelength is termed an action spectrum. 

The calculation of the amoimt of BHuv involves the siimmation of all wavelengths 
weighted by the action spectrimi. The optical amplification factor is calcula- 
ted by dividing the percentage increase in BHuv by the percentage decrease in 
O^. The true optical amplification factor for 0^ reduction varies with 

initial 0^ level, solar zenith angle, and percentage 0^ reduction and on the 

action spectrum for the particular biological effect. The action spectrum 
for most plant effects and for hman skin cancer induction is unknown; however, 
many photobiologists believe that the human sunburn action spectrum is an 
adequate approximation for most purposes. Because of the variation and the 
unknowns, the optical amplification factor of 2 is widely used. It should 
be emphasized, however, that near the Equator where the current levels of 
uv-B are already high, the natural average 0^ level is low, and the solar 

zenith angle at midday is low (Sun nearly overhead) , the true optical ampli- 
fication factor for a smaller 0^ reduction (e.g., 1 percent) may be closer 

to 1.2 to 1.3 for sunburning uv than it is to 2. 


®NASA Lyndon B. Johnson Space Center. 
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Applying the l-percent/2-percent rule of thumb, 0.2 percent 0^ reduction 

obviously converts to a O.^^-percent incr''ase in BHuv. It should be emphasized 
that such an increase in BHuv radiation does not necessarily mean that a 
O.i^-percent increase in some biological harm wilt occur. The biological 
response to BHuv is not necessarily linearly related i.o dose; in some cases, 
the dose of BIIuv has to be greater than some threshold fimount before it causes 
a perceptible effect (e.g., sunburn). 

The impact on the biosphere cf a 0.1*-perrent. inci'cane in BHuv radiation 
can be assessed from extrapolations of the re luti'-ely i.imited 'urperimental data 
obtained with high doses of simulated solar BHuv (ref. H-U). In addition, 
a more intuitively convincing assessment of the eftects of such an increase 
in BHuv can be based on observations made during "nature's experiment,"; the 
effects or lack of effects on organisms and ecosystems of exposures to the 
natural, wide variations in solar BHuv irradiances . The highlights of such 
an assessment follow. 


First, some organisms can be eliminated from concern. Not all organisms 
in the biosphere are vulnerable to an increase in BHuv. Most animals, including 
invertebrates, avoid the Sun by living underground or in the shade during at 
least the heat of the day. The important vulnerable components of the biosphere 
are human beings, unshaded agricultural plants and penned livestock, and natural 
terrestrial vegetation. 


The effects of a O.it-percent increase in BHuv on the vulnerable compon- 
ents will not be detectable. In the first place, the current natural levels 

of BHuv fluctuate widely because of fluctuations in 0^ levels^*^ and because 


of clouds and suspended aerosols. Second, even in controlled experiments, the 
responses or organisms to BHuv are highly variable; biological systems are 
"noisy detectors." Coefficients of variation of 15 to 20 percent are common. 
Third, superimposed on the effects of BHuv are the effects of other factors 
to which plants are exposed such as temperature, moisture, nutrition, com- 
petition, and predation (e.g., refs. H-5 and H-6) . However, the fact that a 
•ause-and-effect relationshix< for an increase in BHuv cannot be established 


^^Measurements from the BHuv instrument on the Nimbus satellite show 
that uhe average 0^ over the entire globe fluctuated during a period from 

April 20 to July 20 with a standard deviation of ±2.k percent. The range of 
values from the mean (31 kilopascal-centimeters ( 306 milliatmosphere- 
centimeters ) ) was ±6.2 percent. On one particular day in April 1969, there was 
a 300-percent difference between the maxiraiun 0^ level (6l kilopascal- 

centimeters (600 milliatmosphere-centimeters) in northern Russia) and the mini- 
miom 0^ level (20 kilopascal-centimeters (200 milliatmosphere-centimeters) near 

New Guinea) . Measurements from ground stations show that variations of 
approximately ±10 percent can occiir in a single day with the movement of a 
weather front. The monthly averages at individual stations (e.g., Arosa, 
Switzerland) vary by ±30 percent throughout a year. The year-to-year range 
for a particular month can be approximately ±5 percent. 
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does not mean that an increase in BHuv can be ruled out as a contributor to 
some deleterious event; lack of detectability of an effect is not equatable 
with a lack of effect. 

Even an undetectable effect can be estimated by extrapolation from ex- 
periments at high dosages, although such estimates have a great uncertainty. 
Such an analysis of experimental data {refs. H-T and H-8) indicates that 
the probability is very low that an average O.U-percent increase in BHuv would 
have any effect at all on plants or ecosystems. It can be demonstrated experi- 
mentally that unnaturally high levels of BHuv can be tolerated by many organ- 
isms with no significant effect. On the other hand, it can be demonstrated 
experimentally that current levels of solar BHuv can be detrimental to some 
organisms that thrive in nature. These findings indicate that organisms have 
evolved strategies for coping with the harmful uv radiation to which they are 
exposed. Many of the mechanisms involved in the strategies for coping have 
been documented. Some examples are 

1. Behavioral avoidance (e.g., nocturnal behavior or complete avoidance 
by living underground, deep in the water, or under trees) 

2. Screening of critical target molecules by hair, feathers, thick skin, 
pigment, or nonessential absorbing molecules in cells 

3. Enzymatic systems that repair solar -uv-radiat ion-induced damage in 
critical molecules or that replace damaged molecules 

U. At the population level, replacement of killed organisms by repro- 
duction of survivors (This mechanism implies that some other factor limits the 
rate of popixlation growth. ) 

It should be remembered that some organisms need tolerate only certain peak 
levels of BHuv for given times of the year because they may be dormant, re- 
sistant, or reduced in ecosystem importance at other times (e.g., grasses 
dry out during summer droiight in some regions and various plants normally 
succeed one another in predominance throughout the year). Other organisms 
(such as trees) may acciamulate damage over the entire year and on a year-by- 
year basis; but, in this regard, many trees, even evergreens, replace their 
more sensitive elements, their leaves, on a routine basis. In many cases, 
the buds for new leaves are in protected regions of the plant. 

The key question in assessing the potential impact on a species of an 
increase in BHuv is whether the strategies for coping will be overwhelmed by 
a given increase or the excessive cost of coping may prevent continued com- 
petition. If nature as it exists today is assumed to represent the unperturbed 
situation, some estimates can be made about whether strategies for coping will 
be overwhelmed by uv increase from a given percentage reduction in the 0^ 

column thickness (e.g., by the 0.2 percent reduction at equilibrium from 60 
Space Shuttle flights per year using solid rocket boosters with aluminum/am- 
moniiam perchlorate ) . 
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If an individual organism can tolerate the maximum amount of BHuv currently 
irradiating its locality at any particular time of year, then by deduction, the 
organism should be able to tolerate anything less than the maximum irradiance. 

By the same token, if an organic species can tolerate the BHuv at a low 
latitude and if all other conditions are equal, other members of the same 
species at higher latitudes (lower amount of BHuv) should be able to tolerate 
an increase in BHuv up to the tolerable amount demonstrated for the species. 
Clearly, all other conditions are never exactly equal, but certainly within 
the range of a few dozen miles, similar organisms under similaj’ conditions 
(e.g., plants in similar soil and moisture) qualify for comparison. In this 
context, it is noteworthy that a 0.2-percent 0^ reduction over a point at 

latitude 35° N will lead to a maximum daily dose of BHuv (in July) that corres- 
ponds to the dosage currently received at a point only 23 kilometers (lU 
statute miles) south (calculated from a table in ref. H-3). 

It should be emphasized that this consideration reflects a maximum 
calculated dose for cloudless skies; the frequent cloudy conditions in stammer 
months are not considered. Moreover, these maximum daily doses are based on 
average 0^ coltomn thicknesses; there will be temporary excursions beyond these 
levels . 

These considerations emphasize the great variability of BHuv doses to 
which organisms are currently exposed without apparent detrimental effects. 

It seems reasonable that organisms in a given ecological niche are adapted 
or adaptable to cope with much more BHuv (up to some maximum tolerance level) 
than they are exposed to on the average. Organisms in nature will only show 
some response to an increase in BHuv when it exceeds their maximum tolerance 
level. Likely, the maximum tolerance level might be exceeded only when the 
increase in BHuv from some anthropogenic 0^ reduction is added on to a natural 

maximum BHuv irradiance resulting from a normal 0^ minimoom occurring for- 
tuitously at the time of maximum insolation and under cloudless skies. 

Conversely, it might be argued that organisms respond to average BHuv 
levels. If so, then an increase in BHuv will add to the average and affect 
all exposed organisms at all latitudes. However, the similarity of physio- 
logical function in members of species living under widely divergent average 
BHuv irradiances suggests that organisms in general do not respond to average 
BHuv levels. Even complex ecosystems can be quite similar under very 
different average BHuv irradiances. Thus, it seems more likely that the max- 
imoom value of BHuv irradiance is the important factor. 

Solar BHuv is probably not a limiting factor in the distribution of living 
organisms. (However, examination of high-altitude plants near the Equator 
may reveal some places where BHuv may be limiting.) Animal life is generally 
limited by food availability, which basically means vegetation. Total vege- 
tative biomass is limited primarily by water; the new dry weight of above- 
ground biomass produced each year by plants is 0.001 of the weight of the annual 
rainfall (ref. H-6). This relationship holds for regions as diverse as the 
grasslands of southwest Africa, the creosote bush deserts of southeni Cali- 
fornia, the beech forests of Central Europe, and the tropical rain forests of 
the Ivory Coast of Africa. The relationship generally holds eveiywhere severe 
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limitations of temperature or moist'ore are absent. Thus, biomass and produc- 
tivity involving natural vegetation and associated animal life are very much 
more limited by wa.ter and temperature (see the following paragraph) than they 
are likely to be by BHuv. 

Of course, the species composition or community structure of the biomass 
may respond to BHuv as one determining factor. Competition usually determines 
the community structure. If one organism in a particular niche has a slightly 
greater competitive edge, it will eventually predominate; but it will not 
necessarily eliminate other organisms. A change in some environmental factor 
may favor a different organism and thus lead to a community structure shift. 
Such shifts already occur normally and repeatedly in response to the niamerous 
changing and interacting factors to which the component organisms of an eco- 
system are exposed. Such changes are usually not detrimental, but they can be. 
For example, nutrient enrichment of lakes may favor algal blooms and micro- 
organism growth that kills game fish. Solar BHuv radiation may be one factor 
in determining community structure; however, because of the n\iraerous inter- 
acting factors, it is not likely that the influence of a small change, such as 
a 0.i+-percent increase, can be estimated. Accordingly, it is probably im- 
possible to establish a cause-and-effect relationship for a BHuv increase in 
any particular docvunented community structure shift that may occur. As dis- 
cussed previously, it is more reasonable to implicate the major essentials, 
temperature and moisture, than it is to implicate BHuv. 

As with natural ecosystems, agroecosystems are currently subjected to a 
much wider range of BHuv irradiances than would occur with an increase in 
BHuv from a 0.2-percent 0^ depletion. Extensive studies have shown that 

given sufficient water by natural means or by irrigation, the factor that 
limits the northern extent of growing regions for partictalar crops is usually 
temperature (including the number of frost-free days). The factor that 
usually limits the southern extent is also temperature. Yields of some plants 
decrease at lower latitudes partly because of water shortage and partly 
because of higher rates of respiration (and, therefore, less storage of 
photosynthate) at the higher summer temperatures. Some plant species are also 
limited by excessive amounts of light. In these species, it is primarily the 
visible and near-uv components that are responsible; the high-intensity 
visible light causes photo-oxidation of plant constituents. Usually, these 
plants are adapted to living under canopies of trees. 

In addition to nature's experiments, preliminary short-term experiments 
performed in the field and in the laboratory show that many agricultural 
plants (15 of 2k tested) exposed to very high doses of simulated solar BHuv, 
under conditions otherwise approximating nature, showed no significant effect 
on dry weight productivity (refs. H-T, H-8, and H-9). If they showed no re- 
sponse to conditions simulating more than a itO-percent reduction in 0^ for 

the particular area, then it is unlikely that a 0.2-percent 0^ reduction will 

have an effect. In the organisms that did show an effect, there was a varying 
and sometimes inconsistent response to the high BHuv dose. Even in these 
cases, the response (percentage decrease in dry weight relative to controls) 
was not proportional to the increase in BHuv. Extrapolation using the 
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conservative assiomption of linearity of effect with dose indicates that any de- 
crease in productivity by the susceptible plants occurring as a result of a 

0.2-percent 0^ reduction and a O.U-percent increase in BHuv would amount to 

approximately 0.1 percent. This analysis leads to the conclusion that there 
is a very low probability that an average O.U-percent increase in BHuv will 
have any effect on agricultural plants or natural ecosystems. 

A 0.2-percent 0^ reduction and a O.U-percent increase in BHuv may lead to 

some increase in the incidence of nonmelanoma skin cancer among susceptible 
individuals. The following factors support an association between nonmelanoma 
skin cancer incidence and BHuv. 

1. Skin cancers are associated with exposed areas of skin (head, neck, 
arms, hands). 

2. Less skin cancer is found among pigmented races than among Caucasians. 

3. Among Caucasians, skin cancer incidence is associated with 

a. Decreased pigmentation 

b. Relative inability to tan 

c. Tendency to sunburn 

d. Increased exposure to the Sun (e.g., outdoor occupation) 

e. Increased light intensity (closer to the Equator) 

U, Genetic diseases (albinism, xeroderma pigmentosum) which predispose 
victims to greater skin sensitivity to solar-uv-radiation damage also pre- 
dispose them to skin cancer induction. 

5. Skin cancer can be produced experimentally in albino and hairless 
mice and albino rats with uv radiation. 

Because of the great variability of BHuv and biological responses mentioned 
previously, such an increase in skin cancer incidence will not be detectable. 
Moreover, the long latent period (20 to 60 years) for induction of skin 
tiimors and the many other factors that already may be tending to change skin 
cancer incidence will also prevent detection. In this latter regard, epi- 
demiological data indicate that certain types of skin cancer are already in- 
creasing independently of any known changes in BHuv. Such increases may con- 
tinue. Factors possibly leading to an increased number of cases in the 
United States include the increased proportion of older people in the popu- 
lation; increased reporting of skin cancer because of Medicare; changing 
lifestyles, which currently involves more leisure time activity in the sun- 
shine; and the net southward migration of the population. (Between 19^0 and 
1970, the center of population in the United States moved westward and approxi- 
mately 55 kilometers (35 statute miles) south from latitude 39° N; this move- 
ment corresponds to approximately a 2-percent increase in annual dose of BHuv. 
Conversely, future decreases are also conceivable. Factors possibly leading 
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to a decreased incidence include action "based on publicity-induced recognition 
of dangers of overexposure to the Sun and more accurate identification of 
susceptible individuals as a result of research generated by the 0^ reduction 
problem. 

Despite the many iinknowns involved in the relationship of BHuv and skin 
cancer incidence, several groups of investigators have made sincere attempts 
to estimate the possible increase in the number of skin cancer cases that might 
result from an 0^ reduction. The most widely publicized estimate is the 

following (from ref. H-IO). 

" Estimated Increases in Skin Cancer for Reductions in Equilibrium Ozone 
Concentrations — Although experimental studies on animals have demonstrated 
the link between UV-B radiation and skin cancer, human epidemiological data 
provide the best basis for quantitative estimates of cancer incidence in humans. 
The current understanding of the skin cancer development from an epidemiologi- 
cal standpoint, however, is somewhat limited. Answers to many important 
questions, including the identification of as yet unidentified parameters 
contributing to skin cancer, are needed to refine estimates. 

"Estimates "based upon changes in the observed cancer incidence with variations 
in latitude for each percent reduction in the equilibrium concentration of 
ozone range from 2,100 to 15,000 (6,000 median) new cases of non-melanoma skin 
cancer per year in light-skinned individuals in the U.S. at steady state. 

This conclusion is based upon the following assumptions (Urbach, personal 
communication, 1975): 

"(l) A l.U to 2.5^ (2^ median) increase in UV-B radiation 

will occur at the Earth's surface for a 1% reduction in 
the equilibrium ozone concentration of the stratosphere. 

"(2) For each Vfo rise in UV-B radiation, a corresponding 

increase in the incidence of non-melanoma skin cancers 
can range from 0.5 to 2% (l^ median) at steady state 
(i.e., after the latency period). 

"(3) Items (l) and (2) combined lead to a O.T to 5% {2% median) 
increase in non-melanoma skin cancers at mid-latitudes for 
each 1^ reduction in the equilibrium ozone level. 

"(U) Based upon recent National Cancer Institute (NCl) estimates, 
there are approximately 300,000 new cases of non-melanoma 
skin cancer in the U.S. annually. 

"(5) Therefore, a 1% reduction in the equilibrium level of 
ozone in the stratosphere might lead to from 2,100 
to 15,000 (6,000 median) new cases of non-melanoma 
skin cancer per year among light-skinned individuals 
in the U.S, at steady state." 
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This approach has been used to estimate skin cancer effects of various 
human activities. For example, the U.S. Department of Transportation (DOT), 
in the addendum (ref. H-ll) to its Environmental Impact Statement (EIS) on the 
Concorde supersonic transport (SST) aircraft, employed the same conceptual 
approach (but a baseline figure of 250 000 cases per year) in calculating the 
increased number of skin cancer cases resulting from various activities in- 
cluding Shuttle. The DOT calculations are based on the Climatic Impact 
Assessment Program (CIAP). The linreasoning application of this approach to 
Shuttle effects is incorrect, as outlined below. 

The mmiber of skin cancer cases resulting from a 0.2-percent 0^ reduction 

lasting for 10 years is not realistically predictable from the existing models 
and data. The extant experimental and epidemiological data, although suggesting 
a relationship between solar radiation exposure and skin cancer incidence, 
are inadequate for making quantitative predictions with reasonable limits of 
uncertainty for this case. 

The data are inadequate for a variety of reasons. Experimental studies 
with mice indicate that skin tumor induction is a complex multistep process 
involving some unresolvable \incertainty (ref. H-12) and having systemic and ^ ^ 
immunological aspects in addition to simple response to BHuv dose (ref. H-13). 
These and other data indicate that not all skin cancer is induced solely by 
BHuv. One-tenth to one-third of the basal cell carcinomas occur on areas of 
the skin receiving little or no solar radiation (refs. H-lU and H-15). The 
non-uv- induced fraction must be ascertained and accounted for to enable the 
relationship of incidence and BHuv dose to be established. The existing 
epidemiological data are inadequate because they were collected in studies 
that were not specifically designed to determine the relationship of skin 
cancer incidence to BHuv dose or to latitude; thus, there are unexplained in- 
consistencies both within and between studies (ref. H-l6) . Many important 
aspects (lifestyles, genetic background, previous history of tumors) were 
omitted from the questionnaire; inclusion might have allowed the resolution 
of some of the inconsistencies. The factor on which the expected increased 
number of cases is based (i.e., the present incidence in the United States) 
is not accurately known; estimates range from 120 000 cases per year (ref. H-17) 
to 300 000 cases per year (ref. H-li*) to itOO 000 (ref. H-l8). The value of 
300 000 cases per year is considered in reference H-l4 to be the most likely 
figure. 

The fact that many of the cases arising in a particular year occur in 
patients who have had a previous history of skin cancer suggests a suscep- 
tibility to skin cancer within the population. This suggestion is borne out 
by the finding that certain ethnic groups (Celtic peoples from Ireland and 
Scotland) are highly susceptible. If only a certain fraction of the total 


M. L. Kripke and M. S. Fisher: Effect of UV Light on the Host Response 

to UV-Induced Tumors. Abstracts of 4th Annual Meeting, American Society of 
Photobiologists, to be published. 
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population is susceptible, the postulated increase in skin cancer cases may 
primarily occiir in the already suffering fraction. This point should be con- 
sidered in estimating skin cancer incidence increases. Data needed for such 
an estimate are not yet available. The estimated harm to the total population 
is thus inflated. Although counting cases, as distinguished from counting 
people, may inflate the estimated nijmber, factors in the analysis of the 
epidemiological data also deflate the number. The 300 000-case/yr figure is 
reputedly based on the factor 165 /IOO 000 cases/yr, a demographic ally weighted 
average incidence for the entire United States. This factor times the 1970 
white population figure equals a number that, rounded off, equals 300 000 
cases/yr. Other factors tending to obscure the true incidence of skin cancer 
is that skin cancer incidence is frequently underreported because of outpatient 
treatment by general practitioners. In some studies for which biopsies were 
not taken for histological confirmation, the incidence was overreported be- 
cause of inclusion of noncancerous actinic keratoses of self-healing kerato- 
acanthomas in the statistics (refs. H-I 6 and H-19). 

Despite the uncertainty in the data, a variety of mathematical models 
have been proposed that represent attempts to fit the epidemiological data. 

Prom these models, amplification factors (percentage increase in skin cancer 
incidence per percentage decrease in 0^) have been derived. Variation in 

amplification factors from 1 to 6 have been reported. These amplification 
factors are very sensitive to underlying assumptions of the models. 

The foregoing considerations lead to the conclusion that current data 
and models are inadequate to enable accrrrate estimations, within reasonable 
limits of uncertainty, of the degree of skin cancer increase that might result 
from a small, temporary 0^ reduction. Of particular importance is the fact 

that the Shuttle perturbation lasts about 10 years, which is a fraction of 
the average induction period for skin cancer. Lack of knowledge therefore pre- 
cludes an acciirate estimate of the nimiber of additional skin cancer patients 
that might result from a 10-year-duration, 0.2-percent 0^ reduction. Given 

that the 0.2-percent 0^ value will be reached by 1985 if space transportation 

systems continue as scheduled, given that the 0 ^ will rapidly return to previous 

levels, and given that the induction time for skin cancer is 20 to 60 years, it 
is concluded that any increase in skin cancer induced by proceeding as scheduled 
for another decade will be extremely small compared to normal incidence. 
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